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Perturbation theory is applied to nuclei in the cen- 
trally symmetric representation in which the neutrons and 
protons have the wave functions of three-dimensional har- 
monic oscillators. In first order the particles are indepen- 
dent, but the second-order calculation removes this over- 
simplification. In order to calculate the kinetic energy 
correctly, a transformation to the coordinate system in 
which the center of gravity is at rest is introduced. The 
question of convergence has two aspects: whether the 
successive contributions, first, of more and more highly 
excited states, and second, of the successive higher orders, 
diminish rapidly. One of the many interaction assumptions 
which are equivalent for He* is used in calculating the 
binding energy of He‘, with a result only slightly less than 
that given by the equivalent two-body method, a satis- 
factory proof of both methods. The calculation of the Li® 
binding energy with one form of interaction is carried far 
enough to include the second-order contribution of the 


sextuply excited states and the third-order contribution 
of the doubly excited states, the convergence being ap- 
parently sufficiently rapid that further contributions would 
be negligible. The change in some of the smallest of these 
contributions effected by altering the interaction assump- 
tion is also neglected, in calculating the Li® binding energy 
with other forms of interaction. All forms of interaction 
considered have a radial dependence resembling the error 
curve, and all but one treat like-particle and unlike- 
particle interactions symmetrically. Of these, the only 
forms which satisfy the demands of scattering and of the 
H? and He‘ energies, and which also give enough binding 
energy for Li®, involve combinations of all types of permu- 
tation operators with rather large positive and negative 
coefficients. The influence of the second order on the calcu- 
lation of nuclear mechanical and magnetic moments, in 
particular those of Li®, is also discussed. 


ETHODS for the calculation of the energies 
of the two-, three- and four-particle nuclei 
from assumed interactions of the particles have 
been developed.! With them and from scattering 
data? it has been possible to get valuable indica- 
tions of the probable approximate nature of the 
interactions. But the number of accurately 
known data depending critically on the assump- 
tions is hardly greater than the number of 
arbitrary constants involved, so that adequate 
verification of the assumed interactions by pre- 
diction of independent results has not been 
1See especially Feenberg and Knipp, Phys. Rev. 48, 906 
(1935); Feenberg and Share, Phys. Rev. 50, 253 (1936). 


*As analyzed most recently by Breit, Condon, and 
Present, Phys. Rev. 50, 824 (1936). 


possible. The statistical method (Thomas-Fermi) 
has been applied to heavy nuclei to test certain 
assumptions, but has not furnished very definite 
criteria.’ It will therefore be important to develop 
a method for treating the binding and other 
properties of nuclei consisting of more than 
four particles. 

The success of the perturbation theory in 
atomic problems‘ suggests a similar attempt in 
nuclei. The convergence in the atomic case 
permits no certain prediction of convergence in 
the nuclear case, for the atom is favored by the 

°C. v. Weizsacker, Zeits. f. Physik 96, 431 (1935); Breit 
and Feenberg, Phys. Rev. 50, 850 (1936); et ai. 


*Cf. Condon and Shortley, Theory of Atomic Spectra 
(Cambridge, 1935). 
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stationary (although, comparing with celestial 
mechanics, relatively very weak) field of a 
central particle. In atomic calculations, however, 
one uses a central field which is, in a Hartree 
type of approximation, at least directly at- 
tributable about as much to the average of the 
perturbations as to the central particle, and a 
similar average field enters nuclear calculations 
also. The manner of carrying out the average, 
in the atomic case, is strongly influenced by the 
stabilizing effect of the central particle, and 
there is no question of disintegration, such as is 
known in nuclei and might manifest itself as 
divergence in a calculation with localized wave 
functions. Formally, in investigating conver- 
gence, one has to ask whether integrals //,», of 
an interaction multiplied by products of different 
wave functions of the same variable, are con- 
siderably smaller than integrals (or differences 
of integrals) Ha. of a sum of interactions multi- 
plied by squares of wave functions. H/» is re- 
duced by cancellation of positive and negative 
parts, somewhat more effectively if the inter- 
actions vary only slightly within one wave- 
length, but considerably in any case. //gq is kept 
within bounds in the atomic cage by the screening 
of the nuclear charge by other electrons, and in 
nuclei presumably by the exchange nature of 
the interactions which reduces the effectiveness 
of the interaction of one particle with many 
others. The two cases have, then, considerable 
similarity in the question of expectation of 
convergence. There are indications that the wave 
function of a nucleus is not with any accuracy 
separable into simple, single-particle, wave func- 
tions, a sort of complexity which has been dis- 
cussed in connection with scattering by Bohr,' 
and which one would expect from a naive picture 
of the problem. There is the alternative picture of 
a heavy nucleus as a group of alpha-particles 
intact, but this is probably an over-simplification 
in the other direction.6 In cases where an 
assumption of extreme simplicity has proved 
satisfactory, as in the Gurney-Condon-Gamow 
correlation of alpha-emission data, all the 
simplicity assumed is, of course, not essential to 
the result. Thus, although one may reasonably 
hope for convergence, one may also expect that 


5 N. Bohr, Nature, 137, 344 (1936). 
®'W. Heisenberg, Zeits. f. Physik 96, 473 (1935). 
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at least the second order of the perturbation 
theory, which represents a mixing of the simple 
wave functions of the initially almost inde- 
pendent particles, will be important for some 
problems. 

In selecting a type of initial wave functions 
with which to calculate, it is essential not only 
that they bear some resemblance to the true 
wave functions, but that they lead to integrable 
expressions when combined with the assumed 
interactions. At least for light nuclei, with 
interactions assumed to resemble the Gauss 
error function, e~*"*”, both requirements are 
met in the representation in which each particle 
is given the wave function of a three-dimensional 
harmonic oscillator. Since the perturbation 
theory assumes that the initial wave functions 
satisfy a zero-order wave equation, we introduce 
a fictitious zero-order potential 


V°= }(ca)?Er?, 


(1) 


summed over all particles. Due to the separa- 
bility of r? into x, y?, and 2*, we might use either 
a Cartesian representation with wave functions 
separated in Cartesian coordinates, or a spherical 
representation with wave functions expressed in 
spherical harmonics. The former seems to be 
simpler for computing binding energies. 

The initial wave functions are still not com- 
pletely specified, the value of the “‘inverse- 
square width” parameter o being arbitrary. The 
result of a complete perturbation calculation is 
expected to converge to the correct energy as 
rapidly as possible when the initial wave func- 
tions are chosen to make the first-order energy 
as good an approximation as possible. The 
second-order contribution being negative, the 
true energy is surely considerably below the first- 
order energy. We therefore select our initial wave 
functions by variation of ¢ for minimum first- 
order energy. It will be seen, in the cases treated 
below, that this leads to practically the same 
second-order energy as does variation of o for 
minimum second-order energy. 

The center of gravity of all the particles is 
not at rest in the coordinate system of the wave 
functions (the oscillations of the individual par- 
ticles about the origin being independent in zero 
order, for example). This coordinate system has, 
then, an irregular ‘‘shuttling’’ motion. The 
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correct Hamiltonian includes the kinetic energy 
>p"/2M in a coordinate system x’, y’, 2’, in 
which the laboratory or center of gravity is at 
rest. For purposes of calculation, this kinetic 
energy may be expressed in terms of the same 
coordinates x, y, 2, as are the wave functions. 
We calculate, in a ‘‘shuttling’’ coordinate system 
for convenience, the average value of the energy 
in a stationary system (in nonrelativistic approxi- 
mation). The slight difference between the 
kinetic energy which we calculate and the 
kinetic energy of the zeroth order wave equation 
enters the perturbation procedure just as does the 
corresponding difference in potential energies, the 
familiar perturbation term. From p’=p—p/JN, 
N being the number of particles of equal mass 
M, we have the wave equation 


1 1 1 


2 a>b 
+U- V+E|y=0. (2) 


Here —U is the total binding-type interaction 
between all particles and V is the Coulomb 
potential of the protons. The units used (as in 
Feenberg’s papers, from which we shall take 
parameters) are mc* for energy (0.51 Mev) and, 
to eliminate this and the factor h?/M from the 
kinetic energy terms, fic'(mM)-! for length 
(8.97 X10-" cm). In these units the Coulomb 
energy is 

V1840_ 1 4/271 


137 r 8 
summed over proton pairs. The perturbation 
term in the Hamiltonian is 


lV 


(3) 


Although the perturbation potential is not small 
(V° being zero, and — U, as we shall specify it, 
having its deepest negative value when all the 
coordinates are zero) the effective ‘‘perturbing 
force” is rather small: Considered as functions 
of one r for example, averaged over the others, 
—U and V° have roughly a constant difference 


in the region where the wave functions are 
large, and a constant in //’ does not affect the 
calculations, because of erthogonality. (The 
rapidly increasing perturbation at large distances, 
the term V°, does introduce a tendency for the 
calculation to diverge in nuclei, but it may be 
more than compensated, in the comparison with 
atoms, by the fact that V° also keeps the de- 
generate sets of highly excited levels relatively 
farther apart than they are in atoms.) 
The zero-order wave equation is 


= {3 Ey =0. (4) 


Using the variables £=(ca)'x, »=(ca)!y, 
¢=(cea)'z, it has the “Cartesian” solutions 


N 
y= II Il Hyp 2 (5) 


a=1 Ent 
simply products of harmonic oscillator wave 
functions, wherein we normalize the Hermite 
polynomials thus: 
Hy(§) 
Ho(&) = —1), 
H;3(§) — 38), 
H,(&) = 2-'6- (4 128? +3), 
(4 — 208° +158), 
(8° — 90E? — 15). 


(6) 


The zero-order energies are 


Dna. 


It is convenient to introduce antisymmetry 
immediately by forming 


N 
p= (N,N,)7! > ( ) (Prt 
Py Py a= 


(7) 
where P, is an operator permuting coordinates 
(including spin) of the N, neutrons, and P,, of 
the N, protons, and P in the exponent is the 
order of the permutation P. 6,,,(¢a) is of course 
a Pauli spin wave function. The usual high degree 
of degeneracy is apparent. It is only in such 
cases as the ground states of the alpha-particle, 
wherein two protons and two neutrons each have 
three m’s equal to zero, or of O'*, having as many 
particles as possible with no m greater than 1, 
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that there is no degeneracy. The perturbation 
may be diagonalized within a degenerate set, 
yielding wave functions ¢ which are linear 
combinations of the ¥°’s within a degenerate set 
and which approximate the proper functions of 
the perturbed system. The first-order energies 
are then EY = fo//¢. In deriving an expression 
for the second order energies E® it is usual to 
seek y expanded in terms of all the ¢’s, which 
gives E® in terms of integrals /¢./1'¢,. But 
the degree of degeneracy of the excited states of 
nuclei is sufficiently high that the selection of 
the ¢’s may become quite involved. This may be 
avoided, however, in calculating properties of a 
state arising from the lowest degenerate set, for 
example, by expanding y in terms of the ¢’s of 
that set and the y's of all other sets. The second- 
order energy is then expressed in terms of 
integrals y,° thus: 


=D Ep’). 
b 


The evaluation of the integrals is, of course, 
unaltered and considerably simplified if one 
omits the antisymmetry from either ¢. or y°, 
modifying the normalization accordingly. The 
third-order term in the energy is 


av’ — bavdT oo’) 
(8) 
a 6 (Eo— E.)(Eo— Es) 

The quadratic array of matrix elements involved 
makes the computation of E® impracticable 
when many excited states must be taken into 
account. But an estimate using a limited number 
of the lower excited states will serve as an 
indication of the rate of convergence of the 
method in individual cases. 


TYPES OF INTERACTION 


For a preliminary explanation of the trend of 
the mass defect curve within the limits of 
accuracy of the statistical model’ and also for 
reconciling the mass defects of the deuteron and 
alpha-particle’ (with very narrow wave func- 
tions), the assumption of binding-type inter- 
actions between unlike particles only was suffi- 
cient. But the greater stability of even numbers 


7 W. Heisenberg, Zeits. f. Physik 77, 1; 78, 156; 80, 587 
(1932); E. Majorana, Zeits. f. Physik 82, 137 (1933). 
*E. Wigner, Phys. Rev. 48, 252 (1933). 
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of protons and of neutrons’ and, recently and 
more directly, the scattering of protons on 
protons,'° indicate the existence of like-particle 
binding interactions. The proton-proton and 
neutron-neutron binding interactions are pre- 
sumably the same, but the former has super- 
posed on it the Coulomb repulsion. Since the 
like-particle interaction is to be taken of the same 
order of magnitude as the proton-neutron inter- 
action, it must also be mostly of an exchange 
nature to avoid collapse of a heavy nucleus. 
One interaction assumption that satisfies these 
demands is that used by Feenberg and Knipp:! 
tron 
2P +1 
— -. (9,) 


like-particle 
pairs, ab 


The permutation P‘® operates’ only on space 
coordinates; P is the more familiar permutation 
which operates on spin coordinates also.’ It 
happens that the values of the parameters B, C 
and g chosen to fit the neutron-proton and 
proton-proton scattering data and the binding 
of the lightest nuclei by the equivalent two-body 
method" are of such a magnitude as to suggest 
that a simpler interaction is equivalent to (9,) 
in problems concerning the ground states of 
three- and four-particle nuclei.' It is 


U=B P+ gP ar}. (9,) 


all pairs 
This interaction is formally the same between 
like and unlike particles, so we may call it a 
symmetric form of interaction (in contrast with 
(9), which is unsymmetric). It is in effect weaker 
between like particles than between unlike par- 
ticles because of the antisymmetry of the wave 
function in like particles; in the nuclei mentioned, 
the pairs of like particles have opposite spin, and 
between them the term in +gP enters only in 
the “exchange integral,” which has a negative 
sign and is otherwise the same as the direct 
integral when the ground state of a three- or 
four-particle nucleus is involved. This with the 
factor (1—g) in the direct integral makes 
the effective depth of like-particle interaction 


*L. A. Young, Phys. Rev. 48, 913 (1935); Guggenheimer, 
J. de phys. et rad. 5, 475 (1934). 

’° Tuve, Heydenberg and Hafstad, Phys. Rev. 50, 806 
(1936). 

" E. Feenberg, Phys. Rev. 47, 850 (1935); Feenberg and 
Knipp, reference 1. 
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(1—2g)B. If this is equal to C, (9,) and (9,) are 
equivalent in these nuclei. It is apparent that 
the same equivalence does not exist for heavier 
nuclei, so we shall be able to decide which type 
of interaction is more nearly satisfactory. Find- 
ing neither adequate, we shall also consider other 
types (Eq. (17)). 


INTEGRATION 


In the calculation of matrices of these inter- 
actions, we encounter integrations of the type 
Sve Ubi --dry. Because of the separa- 
bility of such integrals in Cartesian coordinates 
and the absence of more-than-two-particle inter- 
actions, they reduce to products of integrals 
which we designate thus: 


favea= (1, | 


x | dvll,(& (o+1)+r) 
XID (E/ (a +1) + 


where r=o+2. We note that fasca=faacs =feraa- 
Instead of writing a closed form for these, which 
would, it seems, be rather intricate, we find it 
convenient to list the integrals which we shall 
need in subsequent calculations : 


foovo = 
fiwo=1, T, 
fiono= 1 1 T, 
Soom = =] (24r), 


fin =1—2/7r+3/r? 

fone = (2/r—3/r?) /23, 
= (—1/7r+3/7?) /23, 
fe200=3 (2r?), 

= 1 —2/7r?+ (3/2) /r?, 
fs100= — (3 

Fano = (3 2)4(—1 r+1/r?), 
= (3 /2)*/(2r?), 

foo = 2/7 —6/7?+ (15/2) 
fain = (3 2)i(—1 r+4 r?—5§ r*), 
= (34/2) (3 r?—5/r3), 
= 34( — 1/72+ (5/2) /73), 
Ss300 = / (27), 

Sano = (3/2)4( 
= 
= —345/(4r%), 

fsioo = (15/8) 4/73, 

= (4r3), 


(10) 


THE BINDING ENERGY OF THE ALPHA-PARTICLE 


An impression of the extent and convergence 
of the perturbation calculation can be had by 
examination of the details of this example, 
using, explicitly, the unsymmetrical form of 
interaction.” A zero-order state of the alpha- 
particle can be described in all detail by giving 
sixteen quantum numbers, describing the x, y, 
and z oscillations and the spin of each particle. 
For a certain excited state, we have for example 
the list 100+, 000—; 000+, 100— where the 
numbers give n’s and the signs give the sign of m,, 
those before the semicolon (;) referring to 
neutrons and the others to protons. This can be 
abbreviated 1, 0; 0, 1 if omitted m,’s and m;’s are 
understood to be zero and + is understood to 
come before —. In calculating the diagonal 
energy for a state, the symmetry between the 
particles and the separability into x, y, and s 
integrals, reduces the integration to simple forms. 
For example, for the ground state, 0, 0; 0, 0, 
we have 


(0, 0; 0, 0) 2¥.2/0, 0; 0, 0) 
=12(0 d*/dx? 0)= —60a. 


The terms in ¥,-V, here give nothing, being 
products of odd functions of different variables. 
The average binding potential is 


(0,0;0,0 0; 0, 0) 
= {(4—2g)B+2C} (0, 0 | 0, 0) 
= 


The Coulomb term is reduced by a factor 1/2 be- 
cause of exchange, leaving’ (ac)!/4. We have, 
then, for the first-order energy of the alpha- 
particle, 


E®©+E%=(9 
(11) 


The second-order energy consists of contribu- 
tions from many excited states. The integrals 
Hq’ are zero between the ground state and any 
excited states for which the sum of the m's is an 
odd number, because /7’ is odd in an even number 
of coordinates and the wave function is odd in 
an odd number of coordinates. The states of the 

® Preliminary results for the still simpler case of unlike- 
particle interactions only were given in an abstract: Inglis, 
Phys. Rev. 50, 399. More accurate calculation shows those 


results to be about two percent low. The first order only 
was given by Heisenberg, reference 6. 
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TABLE I. Doubly excited states of He’. 
Hoa 
TYPE NR Hoa’ (o =2.6, etc.) 
2, 0; 0, O 12 {(2—2) B+C)u/r—3ae 8}, y2 —0.19 
1,0; 1,0 6 —Bu/r+aa/8 —1.45 
1, 0; 0, 1 6 | —(1—g) Bu/r+aa/8 —0.09 
1, 1; 0, 0 6 —Cu/r+aa/8 1.27 
++-- 
at. 6 3 —gBu/r — 1.36 


system which do contribute to the second order 
may then be grouped as doubly excited states 
(degenerate with E,°=E,°+2ac), quadruply ex- 
cited states, sextuply excited states, etc. Of 
these, only those with and Yn; sepa- 
rately equal to even integers have nonvanishing 
nondiagonal elements //,,’. The number of con- 
tributing states is further limited by the fact 
that no term in the Hamiltonian involves more 
than two particles, so only one or two particles 
may be excited in a contributing state, and by 
the selection rule A.1Js=0 which confines our 
attention to one value of total spin projection, 
because our operators in this approximation 
introduce no spin-orbit coupling. The Coulomb 
term is for light nuclei very much smaller than 
the rest of the Hamiltonian, so we neglect it 
beyond the first order, which leaves H’ sym- 
metrical with respect to protons and neutrons. 
Because of the equivalence of the particles and 
of directions in space, there are several doubly 
excited states of the same type, having the same 
value of H,,’. The number of states of one type 
is the product of the number of permutations of 
its set of n’s in x, y, and z, multiplied by two if 
protons and neutrons are not excited sym- 
metrically, and again multiplied by two if the 
excitation is not symmetrical relative to + 
and —. The types of states, with the number of 
states of each type and the corresponding ex- 
pressions for //7,, in terms of the parameters, 
are given in Table I. The sequence of m, is the 
same as for the ground state, except for the last 
line, as indicated. For brevity, we introduce 
u=(e/r)'=(¢/(¢+2))'. Terms in ao do not 
occur in H,,’ beyond the doubly excited states 
(cf. the selection rule for the harmonic oscillator) 
The third, fourth, and fifth types listed in Table I 
have matrix elements J/,,.’ quite similar to that 
of the second type, differing, aside from the 
terms in ag, in the replacement of B by (1—g)B, 
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C, and gB, respectively. The same similarity 
recurs through the higher states, and we shall 
shorten the tables by listing only one of the 
four similar types explicitly and indicating the 
omission by an asterisk (*). With this under- 
standing, the types of quadruply and sextuply 
excited states are indicated in Table II. In order 
to be able to compare the contributions of the 
various groups of states, we shall split the second- 
order term in the energy of the ground state 
into a part due to all doubly excited states, a 
part due to all quadruply excited states, etc., 
thus: Adding 
the contributions from the types of states listed 
in Tables I and II, we have 


— (F+G)u?/ (rac) 
—[(2—g)B+C ]u/r+3ac 16, 


— Eg =3(31F +5G)u2/(8r4ac), 


where 


and G=[((2—g)B+C Ff. 


Now we are in a position to obtain a numerical 
second-order result for the binding energy by 
assuming values of the parameters a, g, B, and C. 
It has been stated that the unsymmetrical form 
of interaction (9,) and the symmetrical form (9,) 


excited states of He'. 


TYPE a Nr. (72/u) Hoa’ 
+ - 
4, 0; 0, O 12 (3/8)8{(2—g)B+C} 
22, 0; 0, O 12 —(1/2)}(2—g)B+C} 
3, @; 12* (3/2))B 
12, 0; 1, (1/2)4B 
2, ©; 2, 0 6* —(3/2)B 
2, 0; 02, 0 12* —(1/2)B 
11, 0; 11, 0 6* —-B 
6, 0; 0, O 12 (53/4) {(2—g) B+ /r 
42, 0; 0, O 24 (33/4) }(2—g¢g)B+C{/r 
azz, ©; 0, O 4 (1/2)3}(2—g)B+C}/r 
5, 0; (15/8)1B/r 
32, 0; . 24* — (34/2) B/r 
122, 0; 1, 0 12* —(1/2)B/r 
4, 0; 12* —(5/4)35B/r 
4, 02, 0 24* (33/4) B/r 
22, 0; 2, 0 24* 3(1/2)§B/r 
022, 0; 42° (1/2)§B/r 
3, 0; 3, 0 6* (5/2)B/r 
12, 0 24* — (34/2) B/r 
21, 21, O 32° —(3/2)B/r 
012, O; a, © 12* —(1/2)B/r 


—B/r 


| 


al 
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are approximately equivalent, with the values of 
the parameters found by Feenberg and Knipp. 
We shall alter those values very slightly, so as 
to make the two forms of interaction exactly 
equivalent in this calculation, while altering the 
results very little, by taking a=16, B=72, 
C=42.5, g=0.205. (With these values, the 
equivalent two-body method gives the binding 
energy of the alpha-particle —E=56mc*. The 
values of a and B were kindly suggested by Dr. 
Feenberg, as better fitting newer mass determina- 
tions.) In varying o for minimum (first-order) 
energy, numerically, it is more instructive, and 
somewhat easier, to calculate energies directly 
than to use the equation for the derivative. In 
Table III are given values of the separate terms 
in the energy, calculated in units mc? with these 
parameters, for values of ¢ about the minimum. 
From these values it is apparent that it makes 
practically no difference in the second-order 
energy, whether we vary o for minimum first- 
order energy and calculate the second-order 
term with the value of o so determined, or vary 
o for minimum second-order energy directly. 
The most striking feature of the results given is 
that the quadruply excited states contribute 
much more than the doubly excited states. This 
does not mean, however, that the calculation 
will not converge, since there is a good reason 
for this exception to the expected rule that each 
group should contribute much less than the 
preceding group. The reason is that the terms in 
ag, which are found for the doubly excited states 
only (see Table I), enter the matrix elements 
with opposite sign to that of B or C, and greatly 
reduce their values. (The value of —Ep® 
listed for o = 2.6 is calculated as 44.8 —44.5 =0.3, 
for example.) So the convergence does not seem 
to commence until after the quadruply excited 
states. We take the small value of Es® as 
sufficient indication that we have quite a good 
value (— — E®@ = 53.0mc*) of the second-order 
binding energy without proceeding further. 

The third-order term in the energy (8) is 
laborious to calculate, so we shall not make as 
careful an estimate of it as we have made of the 
second-order, and shall do it only for our present 
choice of parameters and for o=2.6. First we 
confine our attention to the doubly excited 
states, and evaluate the elements //,,’ for them, 


un 


TABLE IIL. Terms of the second-order binding energy of Het 
in units mc*) for values of o about the minimum. 


o 2.4 2.5 2.6 2.7 2.8 
—E 50.4 50.7 50.8 50.7 50.5 
— 0.. 0.6 
—Eg® | 0 1.7 1.5 
—Es 2 0.2, ( 1 0 1, 


as entered in the last column of Table I. Ex- 
amining these numbers, we see that we may, as 
an approximation, neglect the first and third 
types. The second, fourth, and fifth types have, 
respectively, the following values of the ex- 
pression /7,.'—H1,.’ which enters (8): 


1,0; 1,0: 
1,1;0,0: {2(2—g)(1—1/7r)B 
1+, 0+; 1-, O-: {(2(2—g)(1-—2/r) —3/77]B 
+2[1+1/(37) 


The nondiagonal element J/,,’ between the first 
two of the states here listed is (¢+1)Bu/r 
= 24mc?; all other nondiagonal elements between 
these types of states, or between states of the 
same type, are smaller at least by a factor g or 
1/7 (and are not all of the same sign), so we 
may neglect them. Substituting these values 
and those of Table I in (&), enumerating the 
states of each type, and remembering that the 
“nondiagonal”’ terms enter twice because the 
two summations in (8) are independent, we have 
for the part of the third-order term in the energy 
due to doubly excited states only: 


=O0.5¢mc*. 


Of this about —0.04mc? is due to the “‘non- 
diagonal” terms. A rough estimate of the total 
contribution to (8) of all terms involving both 
doubly and quadruply excited states is —0.01mc*. 
The ‘‘nondiagonal” terms being apparently neg- 
ligible, we estimate further only the “diagonal” 
terms of the quadruply excited states, and find 
roughly 
Eg® = —0.0;mc’, 


the principal contributions (—0.1mc*) coming 
from the types 4, 0; 0, 0 and 22, 0; 0, 0. The 
result is at least accurate enough to show that, 
contrary to the situation in the second order, 
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the third-order contribution of the quadruply 
excited states is considerably smaller than that 
of the doubly excited states. Eg is also much 
less than Eg®, which seems to indicate that the 
convergence is rapid. 

Combining the various terms, we have for 
the binding energy of the alpha-particle, as 
calculated with our choice of parameters, 
—E=52.5mc*. The same parameters give 56mc? 
by the equivalent two-body method, which 
probably indicates an upper limit.' Considering 
the binding energy as the difference between an 
average potential and an average kinetic energy, 
the agreement is really quite close, and the slight 
difference between the calculated binding energy 
52.5mc? and the observed value 55mc? indicates 
that the parameters probably should be re- 
adjusted slightly. 


THE BINDING ENERGY OF LI® 


A calculation of the binding energy of Li® is 
complicated, compared to the above, by the 
degeneracy of the lowest states in zeroth order. 
The first-order calculation consists of selecting a 
linear combination of those states in such a way 
as to diagonalize the submatrix of the energy 
between them. There are no matrix elements 
between states with different total projection of 
spin angular momentum, because the operators 
in (9) do no more than interchange spin co- 
ordinates, without changing their sum. We 
therefore have small submatrices for .W/s5=1 and 
for Ws=0. The neglect of spin-orbit coupling 
further means that the total spin and total 
orbital angular momentum separately commute 
with the Hamiltonian, as in atomic spectra, so 
that the states with diagonalized energy are also 
characterized by the quantum numbers L and S 
of Russell-Saunders coupling. The same first- 
order result has been obtained independently by 
Feenberg and Wigner'* and by Bethe and Rose,'* 
by first diagonalizing the angular momenta and 
then calculating the energy in terms of certain 
parameters, for Li® and other light nuclei. The 
diagonalization of the energy may be carried out 
directly for Li’. The energy may also be had 
from the secular problem with our wave func- 
tions, as we shall now show. 


 Feenberg and Wigner, Phys. Rev. 51, 95 (1937). 
4 Bethe and Rose, Phys. Rev. 51, 283 (1937). 
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The nine states with 1Js=1 may be grouped, 
according to their four small submatrices of the 
interaction, thus: (1; 1), (O01; 01), (001; 001) 
with a three-row matrix, and the pair (1; 01), 
(01; 1) having the same submatrix as has the 
pair (1; 001), (001; 1) and the pair (01; 001), 
(001; 01), all of which may be seen from sym- 
metry. After subtracting (0, 0, 1; 0, 0, 1 27,0, 
0, 1; 0, 0, 1) from the diagonal elements, the 
submatrices are 


0 —-a 
and b —d (13) 
—@ 0 i—d b 


where a = Bf 100, and d= Bf", 
using either of the forms of interaction (9). The 
three consequent two-row secular determinants 
each have the roots 6+d, and b—d=a, using 
(10). The three-row secular determinant is easily 
factored by adding rows and subtracting col- 
umns, yielding the roots a, a, —2a. The fivefold 
root @ corresponds of course to the energy of 
the *D, b+d to *P, and the lowest root —2a to 
the *S. The latter root in the three-term secular 
equations leads to the wave function 


(1; 1)+(01; 01)4+(001;001)}. (14) 


This is the ground state, since the interactions 
have been chosen to make triplets lie below 
singlets in cases like this and the deuteron, where 
essentially only one proton-neutron interaction 
is involved. 

The energies of the singlets may be had explicitly by 
considering the states with Ms=0, which fall into groups 
twice as large as those just considered. The 'D and 'P are 
had as {fin — |B from one of the four- 
row submatrices (factoring by adding rows and subtracting 
columns in blocks of two) and the 4S as { —2fiu+(1—2g) 
(Pino—fiom)}|B from the trace (=0) of the six-row sub- 
matrix, knowing the other terms, and using the sum rule 


The ground state has, then, the first-order 
energy (0, 0, 1; 0, 0, 1/ 27/0, 0, 1; 0, 0, 1)—2a. 
Using the symmetrical form of interaction (95), 
this is® 


% The Coulomb energy is calculated here to first order, 
and roughly for other nuclei from He to O having 2+A 
protons, by taking (reference 6) p(r,7’)=2(ao/7)** 
to be (2+3A+43A2/72) (aa)!/8. 


(15.) 


@ 
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Using the unsymmetrical form (9,), we have 
instead 


(15,) 


In calculating the second-order contribution to 
the energy of the ground state, we may again 
group the doubly excited states, for example, in 
types such that all states of the same type have 
the same value of J/,,’. This follows from the 
isotropic nature of the wave function of the 
ground state (14), which means that exchange 
of directions, as well as exchange of all protons 
with all neutrons, in the description of the 
excited state, does not alter the integral. As an 
illustration of the difference between this calcu- 
lation and that of the alpha-particle, which has 
been given in detail, we list in Table IV the types 
of doubly-excited states which we take into 
account. The matrix elements //,,’ are listed in 
a general form for future reference, and they are 
reduced to the matrix elements of our present 
calculation by putting g:=g,=0. 

Adding the squares of the matrix elements of 
these states, and dividing by 2ac, we get their 
contribution to the second-order term in the 


~ 
~ 


energy : 


Ep® = — {52 \/r 
+ 
+214/7*}(Bu }64—126¢ 
+ (154 —60g) /r-+203/72} (Bu/r)/12 
—(323/288)ac. (16,) 


The corresponding computation using the un- 
symmetrical form of interaction (9,) (for states 
having both a neutron and a proton excited, 
the calculation is of course unaltered) gives 


Ep® = — {((63 
+ 329*)o+237 
X (7/27) +35(7/27)? +[2(12 —20¢/3)0? 
+[120?+48¢+63 ]C?} u?/(2aor*) 
+ |B 
—323ac/288. (16,) 


Using the parameters a=16, B=72, C=42.5 
and g=0.205, as for the alpha-particle hive, 
we evaluate (15) and (16) forric e symmetrical 
form in Table V and for the unsymmetrical form 
in Table VI, for values of ¢ about the minimum. 
We see that there is here considerable difference 
in the behavior of the two forms of interaction. 


TABLE IV. Doubly excited states of Li*. 


Type NR 
+ - + - + 
0, O, 3: 0, 1 6 10g; —9g,) 12} 
0, O, 12; 0, 0 1 12 2-4} 10g, —9g,) 
2, O 1; 0, O 1 6 4+9(1—2¢,— 7? |Bu—S5aa/12 
0, 2 0, O, 1 6 2-4} ]Bu—5S5aa/12} 
02, 0, 1 0, 0, 1 12 2-3} JBu—Sao/12} 
0, 02, 1 0 0, 1 12 2-4} [(3 —4g+2g, —3g,)/r+(2—g—4¢: —3g,)/r? ]Bu—5a0/12 
1. 1: 0, O1 12 [(—2+3¢:+3g¢) /r+3(1 — 2g; —2¢,)/7? ]Bu 
0, 11 i; 0, O1 12 —3g,)/r? ]Bu 
0, 2; 0, 2 3 ]|Bu+aa/6} 
0, 0, 11; 0, 0, 11 3 
0, 0, 2: 0, 0, 02 6 2/r2?—7/(27°) Bu 
0, I, 4; 0, oO, 2 6 |Bu+aa/6} 
1; 0, O, 02 12 2-4} Bu 
Ol, 0, 1; 0, 0, 11 6 0 
0, O1, 1; 0, 12 [—g:/ et |Bu+aa/12 
Ol, O, 1; O, 1 12 —2Bu/r+ac/6 
01, O, 1; 0, 1 12 
1; 0, 1, 1 9 — Bu/r+ac/12 
2 0, oO, 1 6 2-4} [(1+g—3g: ]But+aa/6 
0, 1, 02; 0, Oo, 1 12 2-4[(+1—g—g:—2ge)/r+(—14+2g:+2g,)/r?]Bu 
0, O1, 11; 12 Bu+ac/12 
01, O1, 1; 0, O 1 12 (—1+2¢g+2g,) Bu/r+aa/12 
+ + + + - - 
0, 2, 1; 0, Oo, 1 6 —(3g—3g,)/r?} Bu 
0, 02, 1; 12 2-4} 
0, oO, O1 12 Bu 
0, 1, O1; 0, 0, 11 12 | —ge/7+(—¢+2,)/r?} Bu 


| 
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a 0.8 1.0 1.2 1.4 
—-EO—E" 12 13.3 12.7 8.7 
— Ep® 19.3 25.6 
TABLE VI. Unsymmetrical form (Eq. 9.). 
o 1.2 1.4 1.6 1.8 
—EO— Ei 26.9 27.4 6.9 25.1 
—Ep®™ 18.1 17.9 17.6 17.5 


At least to this order of accuracy, the nucleus 
seems to be more unstable and larger (smaller ¢) 
with the symmetric form of interaction than 
with the unsymmetric. This “repulsive’’ effect 
of the symmetrical form may be traced to the 
interaction between like particles with like spins : 
the entire operator |(1—g)P‘+gP} annuls the 
exchange nature of the “exchange integral” 
introduced by antisymmetry, thus augmenting 
its magnitude »‘‘hout changing its sign. 
Comparing these results with those for the 
alpha-particle. we see that the second-order 
contribution to the energy is of a larger order of 
magnitude for Li®. On the other hand, they are 
both calculated as differences of quantities of 
the same order of magnitude (the value —Ep™ 
=17.9 in Table VI is found as 75.6—57.7 for 
example; cf. 44.8—44.5=0.3 for the alpha- 
particle, in Table III above). This may be 
interpreted as an illustration of fact that the 
®ero-order potential V° can be chosen to fit all 
the particles in the alpha-particle very well, 
since they are equivalent particles (namely, 
s-particles), whereas in Li® the size of V® must 
be adjusted as a compromise between the 
tendency of the four s particles to congregate 
closely and the tendency of the more loosely- 
bound # particles to spread out. Thus V° is not 
as good an approximation to the true average 
potential in Li®, and the correction introduced 
by the second-order contribution is larger. In 
more detail, Ep® is easily seen to depend 
essentially upon a mean-square difference be- 
tween the “U part’’ and the ‘“V® part” of the 
elements //,,’ (the kinetic energy correction also 
enters here, but in a secondary réle). This is not 
true of the further contributions Eg, etc., a 
selection rule excluding the “J” part,” as we 
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have seen. Since the positive part of —Ep® 
alone is of the same order of magnitude in Li® 
as in the alpha-particle, we may expect contribu- 
tions of the quadruply excited states, etc., to be 
of the same order of magnitude in the two nuclei. 
In further investigation of the convergence, 
the types of quadruply excited states have been 
listed (but the reader is spared the sight of 
them !), with the number of states of each type, 
and a numerical evaluation of the matrix ele- 
ments J/,,’, for each form of interaction, and for 
approximately the value of o« which gives a 
minimum value to E+E. Computation with 
these elements leads to the following values: 


Symmetrical form of U,c7=1, —Eg® =6.9me’*. 


Unsymmetrical form of U, o=1.4, 

—Eg® =5.;me’. 
To give an idea of the slope of —Eg®(c), the 
value for the symmetrical form and ¢=1.2 was 
found to be about 6.ymc*. The second-order 
contribution of the sextuply states we only 
estimate very roughly, by analogy with the last 
calculation, by noting that half of the contribu- 
tions to —Eg® (symmetrical form of U) are 
made by states in which the s particles alone 
are excited, and a quarter by the types 04, 0, 1; 
0, 0, 1 and 0, 04, 1; 0, 0, 1 and 2, 0, 1; 2, 0, 1 and 
2, 0, 1; 0, 2, 1. The six analogous types of 
sextuply excited states contribute 0.05, and this 
is apt to be roughly the same fraction of the total 
(a survey of the details that make one type of 
state contribute less than another may be 
necessary to convince one of this), which leads 
to the estimate —Es®=0... Experience with 
the alpha-particle has shown that the ‘‘non- 
diagonal” terms of the summation may be 
neglected in reckoning the third-order contribu- 
tion to the energy. Doing this, and proceeding 
otherwise as for the alpha-particle, we find 
—Ep®=-—1.;, as calculated with the sym- 
metrical form of interaction. This we take as 
sufficiently definite indication of rapid conver- 
gence. Since the neglected nondiagonal terms 
are apt to oppose in sign the neglected fourth- 
order contribution to the energy, we take as our 
(most probable!) result the sum of the terms 
we have cited, or 

—E=40..mc?, for the symmetrical form of U 
— E=49.;mc*, for the unsymmetrical form. 
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The approximations made in estimating the 
contributions of some excited states and neg- 
lecting others lead to errors in the result of 
probably not more than Ime’, in the judgment 
of the writer. 

The parameters here used were essentially 
those previously determined by the equivalent 
two-body method. By our method they gave 
52.5mc? instead of 55mc* for the binding energy 
of the alpha-particle. If, in order not to alter 
the deuteron results, we simply change C to make 
the alpha binding 55mc? by our method, it 
becomes C=46. The consequent change in the 
first-order energy of Li® is —9.1mc? and in Ep 
is 1.0mc*?. With the parameters B=72, C=46, 
a=16, g=0.205, the binding energy of Li® is 
then —E=57..mc*, for the unsymmetrical form. 
The result is thus quite sensitive to this small 
readjustment of a parameter to fit the data on 
the lighter nuclei. This particular readjustment 
seems to be incompatible with the newest 
scattering results.'® Other readjustments are 
possible. Fitting the He‘ energy by a change of 
g does not alone have much effect on the Li® 
energy; (0E/dg) being about 60mc* for He* and 
70mc for Li®. Change of B, a and g simultane- 
ously (approximately scattering requirements) 
is a possible means of fitting the binding of He‘ 
(keeping C=(1—2g)B if desired), but it also 
increases the Li® binding not much more than 
the He‘ binding. Values of the parameters 
selected to bring about agreement with the 
binding energies 4.3mc*? of 55mc? of He* and 
62mc of Li® are a= 16, B=72, C=47.7, g=0.238, 
but here C is considerably greater than (1—2g)B, 
which definitely disagrees with the scattering 
results.'° This unsymmetrical form of inter- 
action, although comparatively flexible, seems 
then to be not quite adequate to explain both 
the binding energies and the scattering data. 
The binding energy of Li® calculated by the 
symmetrical form (9,) with our original choice of 
parameters was in more marked disagreement 
with experiment than was the result we have 
just discussed. Furthermore, it is a less flexible 
form of interaction, containing fewer parameters, 
and less can be altered in the Li® result by 
adjusting the parameters to fit He* by this 


% Breit, Condon, and Present, reference 2, especially 
pages 844-45. 


method: we alter a from 16 to 17, and then 
change B from 72 to 75.3 to fit the deuteron 
binding” and g from 0.205 to 0.198 to keep the 
excited level of the deuteron at zero.’ With the 
new parameters we find the He* energy to be 
—55.1, while for Li® the value of —E“—E® 
has been lowered by only 0.3mc? and —E,y™ by 
0.9mc? below the values given in Table V. The 
final value of the binding energy of Li® is thus 
45mc*. This seems to exclude the interaction 
used as inadequate (see discussion below). 


OTHER SYMMETRICAL FORMS OF INTERACTION 


A symmetric form of interaction would add 
to the elegance of the theory of nuclear binding. 
Such an interaction is compatible with the 
scattering results, which seem to require practi- 
cally equal singlet interactions between like and 
unlike particles. Although the special and com- 
paratively simple form (9,) seems to be inade- 
quate, a more general symmet . form may still 
be successful. The general symmetrical form of 
interaction 


U= with 


all pairs 


(17) 


ge) Pa 


(where P..‘” is a permutation operator ex- 
changing spin coordinates only"’) is equivalent”? 
to the forms already used for the two-, three-, 
and four-body problems if g+g,=g)~0.2, the 
constant determined by scattering. We continue 
to specialize J(r.,)=Be-*"«. The first-order 
energy of Li® with this form of interaction is 


E+E 
+(8—2¢g—20g, Bu+0.7(ac)!. 


The matrix elements H,,’ calculated by use of 
this form of U in (3) (neglecting V) have already 
been listed in Table IV, for the doubly excited 
states. The contribution of these states to the 


1 Feenberg and Knipp, reference 1, Table III. 

As an approximation. Cf. Feenberg and Share, 
reference 1 Table IV. According to reference 16, the values 
of go should be larger, thus demanding an increase of B 
and a, to fit Het, which does, however, not improve the 
calculated Li® binding appreciably. 

19 J. A. Bartlett, Phys. Rev. 48, 102 (1936). 

2° Breit and Feenberg, Phys. Rev. 50, 850 (1936). The 
equivalence does not extend to our third order. 
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second-order energy is 


= — 
— 50g: — 102g, — 
(70+ 10g?-+700g,2+ 3942.2 +40¢ — 3502, 
— 240g, — 70gg, — 124g¢,+1050g:¢.)/7 
+ (180-+27¢?-+ 700g, +5872,2—85¢—325¢, 
— 345g,+ 
X (Bu /r)?/ao+ | 64 
+(154—60g — 540g, — 480g.) 203/72} 
(Bu/r)/12—(323/288)ae. (18) 


We wish by use of these expressions to compute 
the binding energy of Li®, after determining the 
parameters, as far as possible, from other 
considerations. The recent paper on scattering 
by Breit, Condon and Present!® indicates that go 
should exceed (by about 0.04) the value selected 
to make the excited level of the deuteron zero. 
This increase of go effects a decrease of the 
alpha-particle binding, which may be compen- 
sated by shortening the range of the interaction 
and increasing its strength.’ Parameters thus 
selected to satisfy the new scattering require- 
ments and to give the observed H? and He‘* 
energies (the latter by this method) are 


a=22, go=0.22 B=92. (19) 


The limits placed on the selection of the g’s 
by Breit and Feenberg's consideration of extreme 
cases of heavy nuclei®® may, for our purposes be 
summarized 

1.25=14+5g,+4¢, (20) 


(which is the second of their Eqs. (7.3), allowing 
0.03 for the Coulomb term). It appears from the 
foregoing that our present task is to try to find 
values of the parameters which will give a low 
enough energy for Li®. Since (20) is a condition 
that certain extremely heavy nuclei will not 
have too much binding, we wil! obtain as much 
binding as possible by selecting the g’s to satisfy 
(20) as an equality. This limits us to a one- 
parameter family of choices of the g’s, the rela- 
tions between them being 


g=0.22-—g.,, 
£:=0.25 —0.8g.. 


Using (19) and (21) we proceed to compute 
the binding energy of Li® for various values of g.. 
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The first-order energy is independent of g, and 
has its minimum E+E® = —23.9mc? at ¢= 1.2. 
Introducing (21) in (18) yields 


Ep® = — {35.8 —5.3g,+20g,? — (31.7 
+80g,”) 
X (Bu/7)?/ao+ 103.8 —6g.+ (5.8 
+12¢,)/7+203/7?} (Bu/7r)/12 
— (323/288)ae. (22) 


This is evaluated for o=1.2 and for several 
values of g, in Table VII. The contribution of 
the quadruply excited states is more tedious to 
compute, but was found, using (9,), to be rather 
small. We notice that Ep® does not differ very 
much, for rather small values of g,, from the 
value computed using (9,), and we shall here 
assume that this is also true of the contribution 
of the quadruply excited states, and of the other 
smaller contributions to the energy. We thus 
subtract about 5mc? from the sum of the values 
just computed, or about 29mc? from the values 
of Ey® in Table VII, to obtain the energy of 
Li® (experimentally —62mc*) for these para- 
meters. It thus appears to be necessary to go to 
rather large values of the g’s (g-~+2 and the 
others from (21)) in order to find a symmetrical 
form of interaction, of a single error-curve radial 
dependence, which satisfies the demands of the 
scattering data, of the upper limit of nuclear 
masses, and of the H?, He* and Li® energies. It 
is still possible (but does not seem very likely) 
that this conclusion would be altered by more 
complete evaluation of the higher order contri- 
butions. 


Note added in proof: The g’s here required are incom- 
patible with further limitations, implied by an extension 
of Wigner’s considerations of the stability of isobars.* 
This seems to exclude (17) as inadequate for a detailed 
theory of nuclei. A symmetrical form of interaction with 
different ranges a~4 for various P’s, or with another radial 
dependence, is, of course, still a possibility. 


DIsCUSSION OF THE CENTRAL APPROXIMATION 


The question arises, whether it would not be 
better to attack the problem of Li® from the 
other extreme of the binding of a deuteron to an 


_* Cf. Inglis and Young, Phys. Rev. 5i, 525 (1937). The 
limitations also remain approximately valid as one departs 
from the high-density limit. 
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TABLE VIL. Ep® in me? for o=1.2 


g 


—1 0 1 2 
Ep” 39.1 24.7 19.7 24.0 37. 


alpha-particle, assuming that they do not distort 
one another very much. One might expect this 
approach to be valid if the binding energy of 
the deuteron to the alpha-particle were much 
less than the “‘internal’’ binding energy of the 
deuteron itself. In the light of this criterion, let 
us compare the situation in Li® and Li’. The 
experimental ratio of the “external” to the 
“internal” binding energy of the deuteron in 
Li® is approximately 3mc*?/4mc*, and the corre- 
sponding ratio for the triton in Li? is 5mc?/16me*. 
It seems from this rough criterion, then, that we 
are more apt to have a merged nucleus in Li® 
than in Li’. In Li’ there are indications™ that 
the lowest state is a *P, the “orbital” angular 
momentum being essentially that of one proton. 
We may conclude that our central approach is 
probably the more nearly correct in first order. 

The problem carried out from that other 
extreme would presumably yield the lowest 
energy of the system as a function of a parameter 
determining the average distance between the 
particles. The interactions are so chosen as to 
make this energy approach —(55+4)mc? asymp- 
totically (from above, the Coulomb range being 
longest) as the average distance approaches 
infinity. If, by our central approach, which con- 
verges rapidly enough to be practicable only for 
smaller average distances, we compute with 
assumed interactions a value of the energy above 
the asymptotic value, we may conclude that 
these interactions would not lead to a merged 
nucleus (if to a stable nucleus at all). The 
physical argument of the preceding paragraph 
makes it probable that such interactions are not 
the correct interactions. 


CALCULATION OF THE NUCLEAR MOMENTS 


Until now we have assumed that the higher 
order corrections to the energies of the low con- 
figuration do not put some other state below 
the *S. Since the first-order separations of the 
states are of the order of magnitude 3mc? and 
the second-order corrections are about 20mec’, 


1D. R. Inglis, Phys. Rev. 50, 743 (1936). 


this is not entirely obvious without further 
inspection. That the singlets should not descend 
below the corresponding triplets seems fairly 
apparent because the interactions have been 
chosen to make matrix elements smaller between 
singlets than between triplets, a characteristic 
which would extend to most of the J/,,’ also. 
Of the other two triplets, let us consider for 
example the *D, which is next to the *S in first 
order. A wave function with which we might 
calculate the energy of the *D is {(1, 1) 
— (001, 001)}/v2 (an arbitrary choice of one of 
the degenerate functions allowed by the three- 
row secular problem above). Using this instead 
of WY, we would construct a table similar to 
Table IV. The type of state in Table IV which 
makes the largest contribution to the second- 
order energy of the 4S is (01, 0, 1: 01, 0, 1) for 
which =(2fio+2A)/y 3. The correspond- 
ing matrix element in the *D calculation is 
+(fioot+A)/v2 for eight of the states of this 
type and 0 for the other four. The same tendency 
for the second-order correction to be largest for 
the S state appears to pervade the rest of the 
table. Noting this tendency, we have without 
further calculation the result that the *S, which 
is the ground state in first order, is the ground 
state in second order also. This is essentially a 
calculation of the nuclear angular momentum, 
I=1. 

In calculating the higher order corrections to 
the zero-order'* magnetic moment, we_ take 
S summed over protons 
only, where 


= got (Eo- E,) + TT 
— 6.0} Hoo’ / | (Eo— Ea)(Eo— 


The first-order correction includes the terms in 
the first power of //,,’, products of the form 
Hoa’ S The first factor 
H,.' has nonvanishing values only for states a 
with an even sum of the quantum numbers » for 
each direction. The second factor has non- 
vanishing values only for states a having =n; 
odd and <n, odd. All such terms therefore 
vanish. In considering the second-order correc- 
tion, the argument must be extended to the 
case where a nonvanishing term would require a 
to differ from zero by even Yn;, and } from a by 
even YAn;, but from zero by odd Thus 
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we see that the higher orders contribute no 
correction to the orbital part of the magnetic 
moment. This rule is quite general, for all 
nuclei, and is of course independent of the error- 
curve shape of the interaction function here 
employed. It follows also from conservation of 
orbital angular momentum, if one thinks of 
expanding in the @’s. 

We have therefore to consider only excited 
states with Ms,=3/2, Ms,=-—}3 or with 
M s,= —3, Ms,=3/2. The former are exactly as 
numerous as the latter, in the special case of 
Li®, and have exactly corresponding elements 
H,,', so their contributions to the projected spin 
magnetic moment, cancel 
one another. In Li® the zero-order result, which 
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makes the magnetic moment of Li® equal to 
that of the deuteron (as observed),” is exact 
(insofar as V is negligible in J/,,.’, cf. reference 14: 
note added in proof). 

In other nuclei correction terms appear due to 
states analogous to the last four types of Table 
IV. These are small, of order g?, for the forms of 
interaction (9), but may be quite large for (17) 
with large g,. There is a remote possibility that 
they might furnish an additional means of 
testing the interaction assumptions. 

I am especially grateful to Doctors H. Bethe, 
E. Feenberg, and L. A. Young for discussion of 
this and related problems, and for communica- 
tion of certain of their results before publication. 


* Manley and Millman, Phys. Rev. 50, 380 (1936). 
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The scattering of slow neutrons by atoms is considered, 
assuming that, in addition to the ordinary nuclear forces, 
there is a magnetic interaction between the neutron and 
the atomic electrons. It is found that the neutrons scattered 
from an unpolarized beam will be partially polarized in 
virtue of this magnetic interaction. Since the scattered 
intensity depends not only upon the intensity, but also 
upon the spin density of the incident beam, the polarization 
thus produced can manifest itself by a second scattering. 
An expression is derived for the neutron intensity after 
double scattering from magnetized iron plates. Under 


INTRODUCTION 


HE magnetic moment of the neutron has not 
been measured directly, but has been ob- 
tained from the magnetic moments of the proton 
and the deuteron.! The assumption of simple 
additivity of magnetic moments, involved in this 
indirect deduction, is, however, open to some ob- 
jection from the point of view of the 6-ray theory 
of heavy particle interactions and magnetic 
moments.” Since the neutron and proton are sym- 


1J. M. B. Kellogg, I. I. Rabi and J. R. Zacharias, Phys. 
Rev. 50, 472 (1936). 

2G. C. Wick, Lincei Rend. 22, 170 (1935); H. A. Bethe 
and R. F. Bacher, Rev. Mod. Phys. 8, 82 (1936). 


optimum conditions, it is found that the scattered in- 
tensity with parallel orientation of magnetizations is 15 
times that with antiparallel orientation. The partial 
polarization of the scattered neutrons indicates that the 
undeviated neutron beam will also have a nonvanishing 
spin density. Expressions are derived for the intensity 
and spin density of a neutron beam after traversing a 
certain thickness of magnetized iron. These results are 
used in the discussion of three types of experiments for 
producing and detecting a polarized beam of neutrons. 


metrical with respect to interaction with the 
electron-neutrino field, the magnetic moment of 
the deuteron should be equal to the “elementary 
moment” of the proton, eh 2Mc. The observed 
value is 0.85 eh, 2Mc, which is probably to be ex- 
plained by the additional moment arising from 
the process of neutron-proton interaction, and by 
the fact that the proton is decomposed and does 
not possess its “elementary moment” during a 
large fraction of the time. 

Recently, Bloch* has suggested a direct method 
of measuring the magnetic moment of the neu- 


°F. Bloch, Phys. Rev. 50, 259 (1936). 
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tron which depends upon the fact that an atom 
may scatter a neutron either in virtue of the 
nuclear interaction with the neutron, or the mag- 
netic coupling between the atomic electrons and 
the neutron spin. It is this magnetic scattering of 
neutrons which we shall investigate. It will ap- 
pear as a result of this calculation that the ex- 
pression for the scattering cross section given by 
Bloch is in error. The difference between the 
results to be presented and his results arise from 
the use of the correct Dirac value of the current 
density and the corresponding magnetic field 
rather than the “‘classical interaction’’ between 
two magnetic dipoles. 


I 


Van Vleck* has shown that, despite the ex- 
change nature of the forces between nuclear 
particles, the interaction energy of a neutron and 
a nucleus may be approximately described by an 
ordinary potential V(r,).6 Consequently, the 
Hamiltonian of a neutron and an atomic system 
of Z electrons may be written: 


1 z 
2M (1) 


where 5Cy is the Hamiltonian of the unperturbed 
atom, ua¢, is the magnetic moment operator of 
the neutron, and @,, @; are, respectively, the 
Pauli spin-matrix vector of the neutron and the 
Dirac matrix vector of the 7‘ electron. 

If we denote the unperturbed energy of the 
atom and the energy of the incident neutron by 
E, and E, respectively, the wave equation for the 
system becomes: 


=(Eo+ E)¥, (2) 


which can be solved approximately by neglecting 
inelastically scattered waves. We therefore write 


Tz), (3) 


where + a matrix in the spin variables 
of the electrons, is the wave function of the un- 


‘J. H. Van Vleck, Phys. Rev. 48, 367 (1935); see also C. 
H. Fay, Phys. Rev. 50, 560 (1936). 

5 In general, it would be necessary to take into account 
the effect on scattering of the virtual levels of the system: 
neutron+nucleus. However, in what follows we shall be 
primarily concerned with the interaction of iron nuclei and 
neutrons of thermal energy, for which no such level is 
known. 


perturbed state of the atom, and y(r,) is to be 
regarded as a matrix in the spin variable of the 
neutron. 

Inserting (3) into the wave equation (2), we 
obtain the differential equation : 


= 


f +++ 7.) 


With the introduction of the abbreviation: 


z 
i=l r-—r; 3 
X awWolhi, (35) 
Eq. (4) becomes 


(h?/2My?+E— V(r))¥(r) = (6) 


which is the wave equation of a neutron interact- 
ing with the nuclear field and the static magnetic 
field of an atom. 

In order to solve this equation we shall treat 
the right side as a small perturbation, that is, 
we write 

V(r) (7) 
where 
(?/2MV°+E— Vir))y(r) =0, (8) 
and 


(9) 


Since V(r) is independent of the neutron spin, 
a solution of (8) can be found in the form of a 
spatial wave function yo(r), times a spin wave 
function which may be expanded in terms of the 
eigenstates x,, of the s component of the spin. 
Hence, 
and 


(h?/2MV?+E— Vir))po(r) =0. (11) 


If we choose the direction of motion of the 
incident neutron as the axis of our polar co- 
ordinate system, the wave function yYo(r) can be 
expanded in a series of Legendre polynomials, 


() | 
t 
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e 
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= > (2/+-1) t'e*——-P (cos 8), (12) 
i=0 kr 
with 
k=(1/h)(2ME)!}. (13) 


The functions u,(7) satisfy the equation 


h? 
( +(e u(r) =0 (14) 
2M\dr* 


and have the asymptotic form 


ui(r)~ sin (15) 


We shall confine our attention to slow neutrons 
since it is only in this case that the long duration 
of the collision compensates for the small magni- 
tude of the magnetic forces. Under these circum- 
stances the neutron wave-length \=27/k is very 
large compared with the range of the nuclear 
forces, and it is well known that all the “‘phase 
shifts,” 6,, will be small except that of the s 
partial wave, do. 

If r be greater than the range of the nuclear 
forces, ro, Eq. (14) reduces to: 


d* I(l+1) 


dr’ 
whose solution, having the desired asymptotic 
form (15), is 


rkr\ 
u(r) = (= ) (cos 


+(—1)! sin (17) 
Hence, if 140, u.(r)(r>ro) is approximately 
Aakr)'Ji.;(kr), ie.; the wave function of a free 
particle with angular momentum /. 69, however, 
is not small and wo(r) becomes sin (kr+6)). The 
wave function which describes the nuclear scat- 
tering is, therefore, for r>7o: 


= F(r, #) = 


l=0 Dkr 
eikr 
Jia Pi(cos v)+- (e250 —1) 
2ikr 
etkr 
=eitr 4 (18) 


2tkr 
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Since ¥“(r) is a matrix in the spin variable of 
the neutron, it can be written 


Substituting this and the similar expansion of 
y(r) into Eq. (9), we find 


? 


M 
8 
h? ; 


2° ms’ 


(20) 


These equations must be solved subject to the 
boundary condition that they contain only scat- 
tered waves. The asymptotic form of the desired 
solution is :° 


Mu,e'*" 
Vm. — | Flr’, DH (r’) 
2rh? r ms" 
)dr’, (21) 
where 


cos © =cos J cos +sin sin cos (g—¢’). 
Here, (r, 3, ¢) and (r’, 8’, ¢’) are the polar co- 
ordinates corresponding to the vectors r and r’. 

The integral occurring in (21) can be approxi- 
mately evaluated by using the expression (18) for 
F(r, 3) since the region in which this formula fails 
gives a negligible contribution to the integral. 
Furthermore, we can, by the same reasoning, 
replace F(r, with since, for slow neu- 
trons, the scattered wave in (18) is comparable in 
magnitude with the plane wave only if 7 is of the 
order of nuclear dimensions. Replacing 
by exp (ixo-r’) and © by exp —(ix-r’), 
where xo and x are, respectively, the propagation 
vectors of the incident wave and the wave scat- 
tered in the direction of r, we obtain 


mi 
(22) 


From the expression (5) for H(r) we find, after 


a simple calculation, 


* Mott and Massey, The Theory of Atomic Collisions 
(Oxford Press, 1933). 
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XE f exp (ie (23) 
=1 


We shall restrict our considerations to the slow 
neutrons which are strongly absorbed in cad- 
mium, the so-called C neutrons. These neutrons 
have energies in the thermal region,’ so that the 
corresponding de Broglie wave-lengths are of the 
order of 10-8 cm. If the neutron temperature is 
sufficiently low, the C neutron wave-lengths will 
be large compared with atomic dimensions.* The 
term exp (7(%»—«)-r;) will, therefore, not vary 
appreciably over the atom and can be replaced 
by the first term in its expansion which gives a 
nonvanishing contribution to the integrals in (23). 
With this approximation, Eq. (23) becomes 
Kok 


Wot(Ko—x) --dr-. (24) 


i=1 


This expression for exp may 
be put in a more convenient form by using the 
approximate relation : 


(xo —%) — X X 
1 - 
2the 
The diagonal matrix element of this equation 
referring to the ground state gives 


Therefore, 


fox (2(%» -r)H(r)dr = 


7E. Amaldi and E. Fermi, Phys. Rev. 50, 899 (1936). 

* The ferromagnetic elements occupy a singular position 
in this respect, since their magnetic moment arises from an 
incomplete inner shell. 
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Since the average magnetic moment of the atom, 
u, is given by the expression 


e 
‘r.)(>oriXa;) 


X (28) 


we obtain, finally: 


fesp (2(%o —%)-r)H(r)d7= 


——-ptirp. (29) 


This expression is valid when the neutron wave- 
length is large compared with atomic dimensions. 
The value of the integral decreases rapidly with 
decreasing wave-length. 

Collecting our formulae, we may write for the 
asymptotic solution of Eq. (4) 


V(r) ~ exp Xm, 


ms 
r ik mes h? Ky—k 
| Ko—k | 
| 
ms ms" 
nel 
Ms, ms" 


Il 


The intensity of the neutrons scattered in the 
direction of the vector r is given by 


(31) 


¥'(r) denotes the scattered wave in Eq. (30) and 
v is the neutron velocity. It is easily shown that: 


2 | Ko 
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where | Cr, |*, 
(a= Cn,’ (33) 
Me, ms" 
are the intensity and the spin density, respec- 


tively, of the incident beam. According to the 
above formula, the total cross section for scatter- 
ing from an unpolarized beam ((@9),=0) is greater 
than (47, k*) sin? 59. All experimental scattering 
cross sections are very small compared with 47 /k? 
at C neutron energies, which indicates that 
If we write 69 =a and neglect squares and 
. (32) simplifies to: 


higher powers of 49, Eq 


To wo 
[== 
hi Ky — 


1 4u,ka 
oh 
Ky 
u— —— *(Go)a— (Go) (34) 
Bloch’s formula, however, has 
1 tu, ka 


for the term in the scattered intensity which de- 
pends upon the spin density of the incident beam. 

The spin density of the scattered beam at the 
point r is: 


1 | 
x~ 
Ko—k 


(Go) a 


Kg 
Ky —k 


= ) 
hs 


Evidently the neutrons scattered from a beam 
whose average spin density is zero will be par- 
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tially polarized. Furthermore, Eq. (34) indicates 
that this polarization is detectable by a second 
scattering. Consider, therefore, the following 
double scattering experiment: A beam of neu- 
trons, whose propagation vector is xo, is incident 
upon a magnetized plate of iron. The neutrons 
scattered in the direction.of x; fall upon a second 
magnetized plate of iron, placed at a distance of 
r, cm from the first, and the intensity of the 
neutrons which are rescattered in the direction 
of x, say, is measured at a distance of rz cm 
from the second scatterer. 

If there are NV, atoms in the fi 
considered for simplicity of negligible dimensions, 
the intensity of the singly-scattered neutrons 


will be 2 


first piece of iron, 


2 
(36) 
ht 


while their spin density will be: 


Io 
(¢:)2= —Ni———— 
ry h?v 
x(- (ne): (37) 


We have denoted averages over the atoms by the 
subscript a, so that (u:)., for example, is the 
average magnetic moment per atom. 

The expression for the intensity of the doubly- 
scattered neutrons is: 


I; 4.107 2 
Te? hi 


4.1/7 No 24unka 
38) 
h 


Substituting the expressions (36) and (37) for J; 


® The following formulae, although derived for amor- 
phous solids, are approximately applicable to crystalline 
substances. 


— 
Io 
h* 
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and (@;),, we obtain 


To | Ko 2 
| ht |Ko— 
417? 4.1/7 
4My,a\7/ 
+( ) (. 
h? | Ko— 


| 
{| ——— (39) 
& 


A convenient set of experimental conditions is 
described in part by the equations: 


Ko =k, 
(Ko—k1)* (Mia =O, 
—K1)* (Meo =O. 


Under these circumstances, the resultant inten- 
sity depends upon the angle between the direc- 
tions of magnetization of the two scatterers. The 
asymmetry ¢, defined as the difference in inten- 
sity between parallel and antiparallel orientation 
divided by the average intensity, is then given by 


4Mu,a\? 
2( ) (wide, (Badu 
h? 
( 4M? %o— )| ) 
ht ida ht (40) 


If the magnetizing fields are sufficiently strong to 
produce saturation, it is permissible to neglect 
[ (xo —v1/ -w)? and replace with 
u. *. The expression for the asymmetry then 
becomes 


amc?*)? 
——, (41) 


(1+ ame?)?) 
X (1+ 
where and 


The saturation value of the intensity of mag- 
netization is about 1700 gauss for iron at ordinary 
temperatures, which corresponds to a value of 
vi,2~2.2. Adopting provisionally the atomic 
beam value of y,= —2, and utilizing the experi- 
mental result we obtain e=1.75. 
From the definition of « we see that the intensity 
of the double scattering with parallel orientatioa 
of magnetizations is 15 times that with anti- 
parallel orientation. However, despite the large 
magnitude of the asymmetry, this effect will be 


Mitchell, C. J. Murphy and M. D, Whitaker, 
Phys. Rev. 50, 133 (1936). 


difficult to detect with present methods because 
of the small intensity of the doubly-scattered 


neutrons if conditions of single scattering obtain 
in both iron plates. 


Ill 


We have seen that the neutrons scattered from 
an unpolarized beam will be partially polarized. 
It is apparent that, as a result of this effect, the 
undeviated beam will also become polarized, the 
amount of polarization per neutron increasing 
with the penetration of the beam. In order to 
investigate this effect, suppose that an unpolar- 
ized neutron beam of intensity Jo is incident upon 
a magnetized sheet of iron. The neutron beam 
will be described by the solution of the wave 
equation: 

| 
=) 
loam 


r; is the position vector of the i‘® iron nucleus, 
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and the symbols V(r), H(r) represent the same 
quantities as in Eq. (6). The resultant conserva- 
tion theorem, applied to a section of the iron of 
area A and of thickness dx, such that dx is large 
compared with atomic dimensions but yet so 
small that J and «, do not vary appreciably 


over it, is 
th 
[ vivdr -as. (43) 
ot 2M 


We may at this point take into account the 
neutron absorption in the iron caused by nuclear 
capture by replacing the left side of Eq. (43) 
with: 


0 
[vivar+ne TAdx, 
ar 


where is the number of atoms per cm and ¢, is 
the capture cross section. When the stationary 
state is reached, 


th 

2M 
If we neglect the waves scattered by other 
atoms in the calculation of y¥/(r—r,) (this is 
equivalent to neglecting multiple scattering), we 
»btain an approximate solution of Eq. (42), viz.: 


¥(r) =exp (ix Cm,Xm,+ LY (45) 


The Cn, are to be regarded as slowly varying 
functions of x in order to describe the decreasing 
intensity of the undeviated beam with increasing 
penetration. Inserting (45) into Eq. (44) and 
evaluating the surface integrals, we obtain: 


dx 


4ryu,ka Kp Ko 


(46) 


If we denote the quantity in brackets multiplying 
I by o and neglect the third term therein, Eq. (46) 
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becomes 


Ky Ko 


su.) “Ca. (47) 
kk 


dI ( 
dx h 


Returning to Eq. (42), we easily derive the 
relation : 


0 
f 
ot 


| 


f (48) 
i 


which, in the stationary state, gives approxi- 
mately 
de. Ko Ko 
= noe, ——— — —— 
dx h?v kk 
1677 
(49) 
Shv 


We shall solve Eqs. (47) and (49) subject to the 
arbitrary initial values 7° and (¢°),. The ab- 
breviations : 


Ua Ko 4h 
=e;, —=n, e3:n=cos 0, ——=6, 
| 3Mva 
4an Maun Wao 
he 
simplify these equations to 
d I’ 
——=(cos On—3e;)-(@’)., 
(51) 
d 
—(e’),=(cos On—3e;)— + 
dé v 
where 
(52) 


These equations can be solved exactly. However, 
the complexity of the resultant solution is such 
that it is convenient to make certain approxima- 
tions which involve the fact that 6~3400 at 


ie 


)- 
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thermal energies. We thus obtain: 


1 
sin © cos O(e"), e:) 
0 
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1 
Xe @)E -(- + 


1 
+— sin © cos (3—cos* 
6 


1 


—-—sin © cos O(e@°),-es cos 


+— sin © cos O(e"),-e; sin 
6 


1 
+-—sin © cos O(e°), 
6 


1 
+- sin cos Of ——(e°), :)es)e cos? @)E 
6 v 


1 f° 1 

+( —+(e°), sin (-) cos O(e"), e:)e: 

v 6 
1 0 

sin (-) cos (-) (o°), e:)es)¢ (3—eos2 ee 


1 0 
-e;+— sin © cos ©- e:) cos 6é 


v 


1 0 
sin © cos sin 6& (53) 
6 v 


where 


1 1 
——(n—cos @e;), e:= (54) 
sin © sin 0 


If we apply these results to the case under 
consideration, i.e.; J°=J, (@°),=0, we find that 
the intensity and spin density of the beam after 
traversing a thickness of x cm of iron are given by 

I= Iye~"** cosh 


4an | 
x x), 


h? 
(55) 
uo 
sinh 
}Ua| 
Mau, | we! 
x ( (3—cos -x). 


The polarization thus produced can manifest 
itself either by an absorption or a scattering ex- 
periment. The first procedure, which we shall call 
the double transmission method, gives the follow- 
ing expression for the intensity of the undeviated 
beam : 


T= cosh cosh Bs 


(Mie (Me) a 
—+— tanh 6, tanh (56) 
| (Made, 


where 


4rn Mau, (wi, 2). 
Bi, 2=(3—cos? ©), 2) ». (57) 


The subscripts 1, 2 refer, of course, to the first 
and second sheets of iron, respectively. The 
asymmetry e, defined as the difference in inten- 
sity between parallel and antiparallel orienta- 
tions divided by the average intensity, is then 
given by 

e=2 tanh 8, tanh Bo. (58) 


The magnitude of the asymmetry which may 
be obtained experimentally is limited by the un- 
desirability of reducing the intensity to such a 
large extent that the effect is hidden by the high 
speed neutron background. Thus, if we do not 
wish to diminish the intensity by more than 75 
percent, we can take 0, = = 30° and use thick- 
nesses of 0.35 cm (x; =x2=0.7 cm), since ¢=12 
10-*."! If both plates are saturated, the asym- 
metry, as calculated from Eq. (58), is 37 percent. 

The second procedure, which we shall term the 
transmission scattering method, leads to the fol- 
lowing expression for the intensity : 


«| -w) | 
| oda he 


Io 4u,Ma Ko 
h? | (wi)a| Ko—k 
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"J. R. Dunning, G. B. Pegram, G. A. Fink and D. P. 


Mitchell, Phys. Rev. 48, 265 (1935). 
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The subscripts 1, 2 refer to the polarizing plate 
and scatterer, respectively; xo and « are the 
propagation vectors of the incident and scattered 
wave; and ¢ is the distance from the scatterer to 
the point of observation. 

The best experimental conditions are obtained 
when 


Ko — 


Ko— Kk) 


Under these circumstances, the intensity is given 


by 
Io 
r hs 
———-y) | 1.) 
24a ht 
To 4u,Ma 
— sinh (ue) a- (60) 
r (wi)a! 


The asymmetry, defined in this case as the differ- 
ence in intensity between antiparallel and parallel 
orientation of magnetizations divided by the 
average intensity, is: 


KATHARINE 


WAY 


SunMa (ue). tanh 


+——[ w (61) 
he 


It is interesting to note that the maximum inten- 
sity occurs with antiparallel orientation of mag- 
netizations, in agreement with what one would 
expect by elementary considerations. 

If both polarizer and scatterer are saturated, 
x,=0.7 cm, and 0,=30°, the asymmetry is 81 
percent. With given values of | | (ue)a| and 
x;, the maximum asymmetry is obtained at 


-0,=90°. For example, under the above condi- 


tions, but with 0,=90°, the asymmetry, as cal- 
culated from Eq. (61), is 92 percent. 

There is still a fourth possible type of experi- 
ment in which a neutron beam is polarized by 
scattering, and then allowed to pass through a 
magnetized iron plate. If the iron plate is of such 
dimensions that it is permissible to neglect the 
fact that the scattered waves are spherical and 
not plane waves, the intensity is given by a for- 
mula identical with Eq. (59). 

In conclusion, the author wishes to express his 
indebtedness to Professor I. I. Rabi and Professor 
E. Fermi for helpful discussions and suggestions, 
and to Professor F. Bloch for an interesting con- 
versation on this subject. 
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Photoelectric cross section curves for a Majorana-Heisenberg potential of the type 
and a velocity dependent potential determined by Jo= —(2B are 
compared with a cross section curve for a square hole Majorana force calculated by Breit, 
Condon and Stehn. In each case the values of the constants used are those which have been 
determined as the best for accounting for the binding energies of H?, H’, and He‘. Results show 
that the cross section values for the first two potentials differ considerably from the third but 
very little from each other. A general formula for the area under the cross section curve, which 
holds for exchange as well as for ordinary forces is derived. For exchange forces fa(v)d(hv) 
2(me*h /2Mc)(1+aa) and this depends only on a, the range of interaction, a being defined by 
a*h? /M =e, the binding of energy the deuteron. The addition of a long range repulsive force to 
the velocity dependent interaction is found to decrease the cross section for this potential type 
considerably. The classical equivalent of the velocity dependent potential operator is 


determined. 
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Cross SECTION CURVES TY] 
HREE hypotheses about the law of interac- 


tion between neutrons and protons have 
been shown to be consistent with the observed 
binding energies of H*, H®, and He*. The require- 
ment that the interaction account for these bind- 
ing energies establishes a best value for its depth 
and range. It is interesting therefore to compare 
the photoelectric cross sections predicted by the 
best values for these three types of interactions 
to see if they are sufficiently different so that an 
experimental determination could furnish a cri- 
terion for deciding in favor of one of the interac- 
tions, and to consider further what general in- 
formation about the neutron-proton interaction 
the cross section can give. 

The three types of interaction considered are 
(1) the square hole ordinary force, (2) the bell- 
shaped Majorana-Heisenberg force, and (3) a 
velocity dependent exchange force. 

(1) Mohr and Massey' show that for the square 
hole ordinary force the masses of H*, H*, and Het 
are reasonably consistent if the range of interac- 
tion is between 1.7 and 2.2X10-" cm, and the 
attraction between like particles between 0.2 
and 0.3 of that between unlike particles. In what 
follows the width 2.0 x 10-" cm has been used. 

(2) Bethe and Bacher? find that for a potential 
V=— the values a = 2.32 X 10-" cm, and 
Vo=34.1 Mev are consistent with the observed 
binding energies of H?, H*, and He*. The depth of 
the like-particle interaction is 21.0 Mev and that 
of the unlike-particle force in singlet states 
20.5 Mev. 

(3) Way and Wheeler*® show that a velocity 
dependent force* whose potential is defined by 
Vy= SJ (x, e)(e)de, where Jo = 
will account for the same binding energies if 
a=1.17X10-" cm and the depth B=48.3 Mev 
for *S states and 26.2 Mev for 'S states. They as- 
sume equivalence for 'S states of interaction 
between like and unlike particles. 

Morse, Fisk, and Schiff* consider a potential of 
the form — with essen- 
tially four parameters 7;, 7%, D, (singlet), D; 

1 Mohr and Massey, Proc. Roy. Soc. A156, 634 (1936). 

* Bethe and Bacher, Rev. Mod. Phys. 8, 145 (1936). 

§ Way and Wheeler, Phys. Rev. 50, 675 (1936). 

‘John A. Wheeler, Phys. Rev. 50, 643 (1936) for 


notation and discussion of velocity dependent forces. 
5 Morse, Fisk, and Schiff, Phys. Rev. 50, 748 (1936). 


Fic. 1. Photoelectric cross section of the deuteron in units 
of 10°” cm? plotted against E, the total kinetic energies of 
the neutron and proton after dissociation, for three different 
neutron-proton interactions. 


(triplet). D, and D, are determined in terms of 
ro and 7; so that they will give a fit for the singlet 
and triplet levels of the deuteron. Since the value 
of the other parameters necessary to account for 
the binding energies of H* and He* have not been 
determined, this type of interaction is not 
included. 

The photoelectric cross sections of the deuteron 
for the best values of the other different interac- 
tion types are plotted against E, the total kinetic 
energies of the neutron and proton after dissocia- 
tion, in Fig. 1. The graph shows immediately that 
there is very little difference for the bell-shaped 
Majorana-Heisenberg (2) and the velocity de- 
pendent (3) interaction. The square hole cross 
section curve (1), however, is considerably lower 
than the other two and experiment might make 
it possible to decide between it and either one of 
the others. The best energy to work at would be 
about E=2 Mev or about hvy=4 Mev. 

The ratios of the different cross sections at 
E=15.3 Mev and 0.48 Mev (hv=2.62 and 17.5 
Mev) are given below. Here again it turns out 
that experimental determination of this ratio 
would help only to decide in favor of (1) or 
either (2) or (3). 

SQUARE (1) BELL (2) VELociTy (3) 
1.55 1.95 2.07 


DETAILS OF CALCULATIONS OF CURVES 


Curve ! 

Although Mohr and Massey consider a square 
hole interaction of an ordinary type, the cross 
section curve given (curve 1, Fig. 1) is for a 
square hole potential of Majorana type since 
this seems the more logical type and its use 
would probably not have greatly altered the 
binding energy calculations. Moreover, for the 
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width a=2X10-" cm, the cross section curves for 
the two types differ very little. For energies up to 
E=8.8 Mev the points on this curve have been 
calculated by Breit, Condon, and Stehn® from a 
formula given by Breit and Condon.’ It should 
be noted that the corrected values of the Majo- 
rana curves are given in the first paper referred to. 
The point at E= 15.3 Mev has been added. Breit 
et al. take the value of ¢, the binding energy of 
the deuteron as 2.2 Mev. For the calculation of 
the other curves the value used was 2.14 Mev. 
The lower value of € would slightly increase the 
ordinates of Breit’s curves. 


Curve 2 

Four points on this curve were calculated by 
numerical integration of the wave equations and 
subsequent numerical integration of the matrix 
element. The point at E = 15.3 Mev is from extra- 
polation with the help of the theorem about the 
areas under the curves given below. In the course 
of the work it was noticed that the wave func- 
tions obtained were almost identical with those 
for a Majorana square hole of width 1.5 times 
that of the 1/e-width of the bell-shaped hole, 
1.€., 3.48X10-" cm. It would be interesting to 
know if this were true in general. 
Curve 3 

The points on this curve can be obtained by 
direct integration from the well-known equation 
for the cross section 


4 
o=— 
c hv 
where M= 


The notation used is the same as that of Breit et a/ 
and is as follows: af’, 1J = «, the binding energy 
of the deuteron ; .1J = mass of proton ; hy =energy 
of photon; No=normalization factor for ground 
state; k°h?/\J=E, sum of kinetic energies of 
neutron or proton after dissociation ; a=‘‘width” 
of interaction, either 1/e-width or end of square 
hole; v=relative velocity of proton and neutron 
after dissociation. 

When /=0 the radial wave equation for the 


® Breit, Condon and Stehn, Phys. Rev. 51, 56 (1937). 
? Breit and Condon, Phys. Rev. 49, 904 (1936). 
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velocity dependent force considered by Wheeler 


1s 


h? —2B 
VW a 


and fo(r) = 


For other values of /, /:(r,p) is assumed equal 
to zero, so that the wave functions for the upper 
states have the simple form 


fe(r) =sin kr kr—cos kr. 


Integration of the matrix element then gives 
for the cross section 


82 e* aa(1+aa) 
- 


3 he (1—aa)? 


1 1 
(1+ 

Increasing the value of a increases this cross 

section considerably. For a=1.17X10-" cm, 

the best value, and hv=2.62 or E=0.48 Mev, 

o=1.57 X10-*’ cm’. At this same energy, with the 

same orders of magnitude, o=2.66 for a=1.87 
and ¢=5.31 for a=2.80. 


AREA UNDER Cross SECTION CURVES 


For ordinary forces it is well known that the 
area under the cross section curve, {o(v)d(hv), 
is independent of the range and type of interac- 
tion and is given by the f sum rule as re*h/2 Me. 
Several writers have pointed out that this rela- 
tion does not hold for exchange forces. It is 
interesting to examine, therefore, the general 
expression for fo(v)d(hv) for any type of force.® 
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* Lam indebted to Dr. John A. Wheeler for pointing out 
the possibility of deriving the following relation. 


d(hv) 


al 


PHOTOELECTRIC CROSS 


Applying Parseval’s theorem 
9 
e 
= f —32 


and on substitution of sy=(i 


M 
2Mc h? | 


2 
=——{1+T7,]} 
2Mc 


when V represents the potential function. If V 
and zs are commutative operators the second term 
is seen to vanish leaving the familiar expression 
of the f sum rule. 

In the general case, expressing Vy as //(x,e) 
¥(e)de the second term becomes 


M 
6h? 


For Majorana forces J(x, ©) = 1’(x)é(x+e). For 
a square hole potential the integration gives for 
the second term 


aa 


T;=—————_ +078?) 
9a*8?(1+aa) 


where h?6?/ 1 = D—«, D being the constant depth 
of the hole. For the bell-shaped Majorana a 
numerical integration is necessary. 

For the velocity dependent force for which J is 
given above 


Table I gives values of 7, for these three types 
of potentials for different values of a. 

It turns out that in all cases 7, is very nearly 
equal to aa so that fo(v)d(hv)=reh(1+aa), 


TABLE I. 72100 percent is the percentage increase in 
So(v)d(hv) for exchange forces as compared with ordinary 
forces. 


da Tz (SQUARE) | 72 (BELL) T2 (VEL) 
1.17 | 0.267| 0.242 | 0.286 

2.32 | 529 495 0.526 | 542 

3.00 | 654 691 
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2 Mc for any of these three exchange forces. A 
determination of the area under the cross section 
curve would thus establish the width of the inter- 
action without giving any information as to its 
type, provided, of course, that the true interac- 
tion is any of the three considered here. An 
ordinary force could be detected by the absence of 
the extra term, aa. Curves (2) and (3) can thus be 
distinguished by the different areas under them, 
the ratio being A;/A2=1.526 1.28621.2. For 
the curves shown, extrapolations can be made 
beyond E=15 Mev with fair accuracy by assum- 


ing that they fall off either as 1, E* or 1, E**. 


LOWERING OF VELOCITY DEPENDENT CROsS 
SECTION 


The only experimental value of the cross sec- 
tion determined so far is that of Chadwick and 
Goldhaber* who find ¢=0.5X10-*? cm? with a 
possible error of a factor of 2 for hy = 2.62. It must 
be remembered that this value includes the 
magnetic dipole cross section which for low 
energy y-rays (2.62 Mev) may be } or } the 
photoelectric cross section depending upon 
whether the 'S level of the deuteron is stable or 
unstable, according to calculations of Breit and 
Condon in the paper already mentioned. For 
y-ray energies above 4 Mev the ratio is 0.05 at 
most. All the calculated cross sections are thus 
above this experimental limit, the square hole one 
being nearest to it. However, more accurate ex- 
perimental determinations are needed to show 
whether there is really a serious contradiction 
between theory and experiment. It seems inter- 
esting in the meantime to calculate the effect on 
the velocity dependent cross section of adding a 
long range repulsive force to the assumed inter- 
action, and it is found that such an addition does 
lower the theoretical cross section for this inter- 
action type considerably. 

The interaction is assumed to be 


Jo=—(2B, “+ (2B, b 


Three of the constants B;, Be, a, and b can be 
chosen arbitrarily. The fourth is then determined 
by a relation with the binding energy of the 
deuteron. The expression for the cross section is a 


® Chadwick and Goldhaber, Proc. Roy. Soc. Al51, 479 
(1935). 
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TABLE II. Velocity dependent cross sections for hv=2.62 
for attractive force alone and attractive force plus two different 
additional repulsive forces. 


| Bi(Mev) | | Bs (Mev) o 
1.17 48.3 0.0 | 0.0 1.57 
1.17 49.3 12.0 0.268 | 0.386 
1.17 49.3 16.6 0.307 0.187 


simple extension of the one already given. In the 
following B,, a and } were chosen arbitrarily, a 
being taken as the best value width, B,; a depth 
very near to the best value depth 48.3 Mev and } 
a width large in comparison with a. o for hy =2.62 
is given for two different values of } in Table IT. 


CLASSICAL EQUIVALENT OF VELOCITY 
DEPENDENT POTENTIAL 


The velocity dependent cross section curve is 
very similar to those of the other interactions and 
the area beneath it obeys the 1+aa rule where a 
is the 1/e-width of the operator 


Jo=- 


although this width does not seem to be exactly 
comparable to that of more familiar interactions. 
It is interesting, therefore, to determine the 
classical equivalent of the velocity dependent 
potential operator. 

Dirac has shown that the best quantum-me- 
chanical representation of a classical interaction 
potential V(y, p) dependent on position and 
momentum is given by an integral operator con- 
nected with V(y, p) by the following equation 


where y= }(x+e). In equivalent form we have 


Vy, p) = f I(x, (xe), 


where y is to be kept constant in the integration. 
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For 

J(X, ©) = 
the integration can be carried out for the case 
sin (y, p)=0 giving 

—8Bh ap 


—¢ tele tan—! — 


V(y, 
ap h 

2y 1 /2y\? ap ap 

+: +(-) 

a 2!X\a h oh 
1 s2y\* 1 2y\4 ap ap 
3'\a 4!\a h oh 


If sin (y, p)=1, the integration is more difficult 
but the resulting expression reduces to a similar 
series of which the first two terms are identical 
with the first two of the above terms. 

For ap/h small, say 0.1 (this would require a 
relative velocity of about c, 30) 


V(y, p)= 


This approximation to the interaction has 
1/e-width of a/2 or 0.585 X10-" cm and a depth 
8 times the “‘effective’’ depth B or 386.4 Mev. 
Bethe and Bacher!’ show, however, that for a 
strictly exponential hole of 1/e-width 0.5 10-" 
cm, the depth would have to be 310 Mev to fit 
the binding energy of the deuteron. Wider holes 
require smaller depths. The above approximation 
is therefore not justified. The velocity dependent 
interaction of Wheeler can be thought of as a 
deep and narrow exponential type of interaction 
modified in a complicated way by the relative 
momentum of the particles. 

The writer is very grateful to Dr. John A. 
Wheeler for much helpful advice during the whole 
course of the work as well as for the suggestions 
about the area under the cross section curves. 


1° Bethe and Bacher, Rev. Mod. Phys. 8, 111 (1936). 
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Unidirectional Measurements of the Cosmic-Ray Latitude Effect* 
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Cosmic-ray intensities in the vertical direction, and in tude effect have been found, one about 12 percent and the 
two inclined directions, 45° eastwards and 45° westwards other about 20 percent, and it is uncertain in the present 
from the zenith, have been recorded as a function of lati- measurements which one is correct. The intensity in the 
tude. The apparatus, embodying several coincidence vertical and western directions is strongly influenced by 
counter trains, was arranged for unattended operation and local anomalies in the earth’s magnetic field, whereas 
automatic recording. It has made several voyages between eastern intensities agree more closely with the distribution 
New York and southern Chile through the Panama Canal _ expected if the earth's field were that of a dipole. Shower 
on ships of the Grace Line. A latitude effect of 12 percent intensities, recorded on one voyage, showed less of a lati- 
was found in the eastern direction compared with 4 percent _ tude effect than vertical rays. 
in the western direction. Two values for the vertical lati- 


the theory of the geomagnetic effects the magnetic effects as recently developed by , 
intensity of cosmic radiation at the surface Lemaitre and Vallartat would enable one to 
of the atmosphere depends upon the latitude and_ calculate the intensity from each direction at 
upon the two angular coordinates required to every point on the surface of an earth whose 
specify the direction. Absorption in the atmos- field were that of an eccentric dipole, if the 
phere introduces an additional dependence upon following data were given: (1) the functions 
elevation and zenith angle. Extensive measure- giving the intensities in external space (outside 
ments have been made of the variation with the magnetic field) and their distributions with 
latitude and elevation of the total intensity, respect to magnetic rigidity, mv,/e, of each type 
integrated with respect to direction, and of the of radiation, positive and negative electrons, 
difference between eastern and western intensi- photons, protons, etc., and (2) the attenuation 
ties at a variety of zenith angles, elevations and factors, functions of energy and atmospheric 
latitudes, but little attention has been given to path length, corresponding to each type of 
variations of intensities from particular directions radiation. The real problem is to reverse this 
with respect to latitude. Some results of this procedure and to make such measurements of 
nature have been reported by Auger and _ the Variation of intensity with direction, eleva- 
Leprince-Ringuet! by Pickering,? and by Clay, tion, and latitude as will lead to the derivation 
Bruins and Wiersma,’ who have measured of the distribution functions and the attenuation 
vertical intensities with coincidence counters factors. Since ionization chambers measure the 
mounted on shipboard for voyages across the intensities integrated with respect to the di- 
equator, but no systematic survey has previously _ rectional coordinates, they are not well suited 
been attempted. to the requirements of this problem and uni- 
The results of such a survey will be of value’ directional measurements which can be made 
in making an analysis of the primary cosmic with coincidence counters are preferable. The 
radiation. Considering the matter from the _ first of these to have been accomplished were the 
opposite point of view, the theory of the geo- measurements of the east-west and north-south 
~ * Preliminary accounts of this work were given before aSymmetries and the variations of intensity with 
the American Physical Society in February and in. De- zenith angle which have now been made in a 
cember, 1936, and before the American Geophysical Union 
in April, 1936. Phys. Rev. 49, 639 (1937); Atlantic City Variety of latitudes and elevations.* As a second 


Meeting of the Am. Phys. Soc., Dec. 1936. Trans. Am. step in the program we are now engaged in a 
Geophys. Union p. 176 (1936); J. Frank. Inst. 222, 647 ae Prog Bas 


(1936). *G. Lemaitre and M. S. Vallarta, Phys. Rev. 50, 493 
'P. Auger and L. Leprince-Ringuet, Comptes rendus (1936). The present state of the theory would enable one 
197, 1242 (1933); Nature 133, 138 (1934). to make this calculation for the zone extending to thirty 
?W. H. Pickering, Phys. Rev. 50, 495 (1936). degrees on either side of the geomagnetic equator. 
’J. Clay, E. M. Bruins, and J. T. Wiersma, Physica 3, ®* These results have been summarized by Johnson, 
746 (1936). Phys. Rev. 48, 287 (1935). 
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Fic. 1. (Top) Photograph of the ‘‘counter box"’ containing the counters and 
circuits for selecting the coincident discharges. The support is from two mutu- 
ally perpendicular horizontal axes acting as a leveling gimbal. (Bottom) The 
‘recording box,”’ containing counter dials and other instruments, flood lights, 


camera, and power converters. 


survey of the variations with latitude of the 
sea level intensities measured at several definite 
directions. 


Tue AvuToMATIC UNIDIRECTIONAL Cosmic-RAY 
RECORDER 


Our intensity meter consists of a box con- 
taining six directional trains of triple coincidence 
counters, mounted on a leveling gimbal, and 
provided with an automatic recording device so 
that it can make extensive unattended voyages 
on shipboard and bring back photographic 
records, taken at hourly intervals, of the dial 
readings of various instruments. These include 
the counter dials, giving the numbers of coinci- 
dences recorded by each counter train, a barom- 
eter, a thermometer giving the temperature 


inside the counter box, a compass for coordi- 
nating readings with course directions, ulti- 
mately taken from the ship’s log, a chronometer, 
and various voltmeters useful as trouble indi- 
cators. In its present set-up two of the counter 
trains are mounted vertically and four are 
inclined ‘at 45° from the vertical in two opposite 
azimuths. No provision has been made for 
automatically adjusting the azimuthal orienta- 
tion of the counter box since on the South 
American run, with which we have been engaged, 
a large part of the voyage adheres closely to the 
meridian. On future voyages this adjustment can 
be made manually by one of the ship’s officers 
and recorded automatically. Power for the 
electrical circuits is taken from the ship's dynamo 
and control devices regulate the voltage. The 
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counter box is shown exposed in Fig. 1 (top), and 
Fig. 1 (bottom) represents the box containing 
the recording equipment, the power supplies, as | o-dmnenen 
and auxiliary apparatus. The arrangement of the « +, 
counters is indicated in Fig. 2. The three ele- * 3 
ments of each triple coincidence train consist of ¢ 
bundles of six counters arranged to cover an he : + ++ oll 
area 20X11 cm’. The outside units of each train 8 
are 50 cm apart so that the angular aperture oo .— 
within which half or more of the area is effective ° 
measures 12° in zenith angle by 23°. A detailed 


account of the apparatus is reserved for another 
publication. 

During the summer of 1935 the equipment 
was mounted on the top deck of the Grace 
Liner Santa Barbara where it completed three 
voyages between New York 41°N and Val- 
paraiso, Chile, 33°S geographic latitude. Because 
of some trouble with the recording apparatus, 
only one of these voyages resulted in a satis- 
factory record. After some improvements in the 
recording mechanism the apparatus was again 
installed in the spring of 1936 on the top deck 
of the Grace Liner Santa Lucia for five voyages 
between New York and Talcahuano, Chile, 
36.5°S and on all of these voyages the records 
were good, although for a part of the time some 
of the individual counter trains have misbehaved. 


Fic. 2. The arrangement of counters. Coincidence trains 
are indicated by rectangles. The ‘‘out-of-line’”’ arrangement 
below is for recording shower intensities. 


GEOGRAPHIC LATITUDE 


Fic. 3. Vertical intensity recorded on one voyage. Counts 
per hour and percentages of the high latitude value are 
plotted against geographic latitude. 


TREATMENT OF DATA 


In reducing the results the average counting 
rates for the twenty-four-hour period from 
midnight to midnight have been corrected for 
barometric fluctuations on the basis of the work 
of Johnson and Stevenson® who found a variation 
of 0.36 percent per mm Hg for the vertical 
intensity. The same correction has also been 
applied to the measurements at 45°, though 
eventually we will have more accurate determi- 
nations of the barometer effect and its depend- 
ence on direction and latitude. Since these 
corrections are always small no appreciable 
error is introduced by the present procedure. 

In the case of the vertical intensities all of the 
readings are used, but a selection has to be made 
of the readings of the inclined counter trains 
since for a part of the time, in port and on 
certain parts of the voyage, the course deviates 
too far from the meridian. In the present report 
we have discarded readings for courses more 
than 45° from the meridian, and for the remain- 
ing readings on a typical voyage the average 
deviations were 18° southbound and 20° north- 
bound. 

As an illustration of the results, the intensities 
recorded by one of the vertical trains on one 
voyage are plotted in Fig. 3. The points represent 
the averages of the twenty-four readings from 
midnight to midnight plotted against average 
latitude. These are expressed as counts per hour 
on one scale and as percents of the high latitude 
value on another. The probable errors, indicated 


® T. H. Johnson and E. C. Stevenson, Phys. Rev. 47, 578 
(1935). 
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560 T. H. JOHNSON 
by the length of the vertical lines drawn through 
the points, were calculated from the actual 
dispersion of the hourly readings and here, as 
with practically all of the probable errors 
calculated in this way, the agreement with the 
estimates based upon the total number of 
coincidences recorded is very close. In spite of 
this agreement the consistency of results recorded 
on different days at the same place was not 
always as good as might have been expected and 
the scattering of points seems to depend in some 
of the records upon other factors than the finite 
number of rays counted. The results have been 
examined for temperature effects by comparing 
day and night readings when the actual inside 
temperatures have not been satisfactorily re- 
corded and, for the results represented in Fig. 3 
and much of the other data, no such effect was 
found. On the other hand in some of the records 
temperature effects and other erratic variations, 
possibly caused by humidity changes, have made 
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their appearance. Such variations have now been 
traced to the high resistors (310° ohms) used 
for stabilizing the counter discharges and we 
expect in the future to be able to overcome 
completely this type of trouble. In the results 
reported herein we believe we are free from 
effects of this nature except for some of the 
results represented in Fig. 4 recorded in and 
near New York harbor where severe fluctuations 
of temperature may have produced excessive 
humidity inside of the counter box. 

The five best sets of measurements of the 
vertical intensity, adjusted to agree at the 
minimum, are plotted against geographic latitude 
in Fig. 4. Here the consistency is good except 
for points above 30°N where the ship leaves the 
Gulf Stream and severe fluctuations of tempera- 
ture are liable to set in. Although no systematic 
dependence upon temperature to the extent of 
the fluctuations in this region was evident, the 
humidity effect noted above would be of a more 
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Fic. 4. Composite results of vertical counters. Different symbols indicate separate voyages. 
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subtle nature and we cannot place too much 
reliance upon the results recorded in the higher 
northern latitudes. Below about 30°N however 
the results seem dependable. The averages of 
the points lying within each five-degree range of 
latitude are represented by the black dots and 
by the smooth curve drawn through them. A 
similar representation of the results of the 45° 
inclined counters in the western and eastern 
azimuths is shown in Figs. 5 and 6. This is 
believed to be the first time that the strong 
dependence of the latitude effect upon direction 
has been given direct experimental demonstra- 
tion although differences between the eastern 
and western latitude effects had to be expected 
in view of the east-west asymmetry found in the 
equatorial belt. The present results indicate a 
western latitude effect at this zenith angle of 
about four percent compared with an eastern 
latitude effect of about 12 percent expressed in 
terms of the high latitude value. The difference 
of 8 percent is in fairly good agreement with the 


value of the asymmetry measured directly at 
the equator.* 

In the case of the vertical intensities the value 
of the latitude effect remains uncertain in the 
present measurements. The experimental points 
seem to divide into two rather distinct branches 
at about 30°N one of which seems to level off at 
this latitude with a total latitude effect of about 
12 percent and the other continues to rise to 
New York harbor corresponding to a total 
latitude effect of 18 or 20 percent. We are not 
able to say at the time of writing which branch 
is correct or what the cause of difference is. 
Although we are looking for some such instru- 
mental effect as that suggested above to clear 
up the discrepancy, at the same time it is perhaps 
not out of place to note that a somewhat similar 
discord exists in the ionization chamber results. 
In his summary of the results of Clay, Hoerlin, 
Prins and of his own world wide survey, Comp- 
ton’ gives a latitude effect of 8 percent between 

7A. H. Compton, Rev. Sci. Inst. 7, 71 (1936). 
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Fic. 5. Composite results of counter trains inclined 45° toward the west. 
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Fic. 6. Composite results of counter trains inclined 45° toward the east. 


the equator and 30°N on the 75th meridian and 
an additional 5.5 percent increase from 30°N to 
41°N. Millikan and Neher® on the other hand 
find an 8 percent increase over the same range 
from the equator to 30°N with less than a 
percent increase from there on up to higher 
latitudes. Allowing for a greater intensity in the 
vertical direction of the softer, field-sensitive 
rays, a greater latitude effect than that measured 
by the ionization chamber is to be expected. 
If, then, our 12 percent latitude effect is asso- 
ciated with Millikan and Neher’s value of 8 
percent, it turns out that our 20 percent value 
when reduced on the same scale just agrees with 
Compton's value of 13.5 percent. 

In the curve of Fig. 4 representing the vertical 
intensity it is noted that the minimum occurs 
not at the geomagnetic equator, 12°S, but 15 
degrees or so to the north of this point. At the 


*R. A. Millikan, Year Book of the Carnegie Institution 
34, 343 (1935); Rk. A. Millikan and H. V. Neher, Phys. Rev. 
50, 15 (1936). 


same time it is noted that the curve is not 
symmetrical with respect to any point of latitude 
as it should be according to the theory based 
upon the simple dipole field. Both of these 
effects, attributable to local field anomalies, are 
here much accentuated because of the uni- 
directional character of the measurements, as 
compared with their appearance in the curves 
representing the ionization chamber measure- 
ments of Millikan and Neher along the same 
route, although some asymmetry with respect 
to the magnetic equator is readily noticeable in 
their results. Clay, Bruins, and Wiersma*® have 
also called attention to similar effects in both 
ionization chamber and counter measurements. 
The intensities measured in the western azimuth, 
Fig. 5, also indicate a distortion due to local 
field anomalies, although here the points are 
somewhat more scattered, partly because of 
lower counting rates. In the eastern azimuth, 
on the contrary, the intensity has a more 
symmetrical distribution, although the minimum 
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Fic. 7. Orbits defining the energy limits for rays incident 
from the east, vertical, and west in the plane of the 
magnetic equator. 


still occurs somewhat north of the magnetic 
equator. 

A dependence upon direction of the distortion 
of the intensity curves due to local field anomalies 
finds a qualitative explanation in the dipole 
theory of the orbits although the introduction 
of the real field into the theory has not yet been 
accomplished. The three critical orbits, corre- 
sponding to positive rays of lowest energy 
incident in the plane of the equator from the 
directions 45°E, the vertical, and 45°W are 
shown in Fig. 7. These have been calculated 
from Stérmer’s equation® with y= —1, which is 
exact for orbits in this plane. A circle on each 
orbit indicates where the ray begins its last 90° 
deflection before it strikes the earth, and it is 
evident that this appreciable part of the total 
deflection takes place in a much more localized 
region in the case of the western and vertical 
orbits than in the case of the eastern orbits. 
Hence any localized magnetic intensity pro- 
ducing alterations in the field of the dipole should 
have more of an influence on the energy limits 
of western and vertical orbits than of eastern 
orbits and this in turn will be reflected in the 
intensity distribution. Actually there is fairly 
close correspondence of the western and vertical 
intensity curves with the local horizontal compo- 
nent of the magnetic field as indicated in Fig. 8. 
In these directions the minimum value of the 
cosmic-ray intensity occurs at the same latitude 
as the maximum value of H7. The eastern cosmic- 
ray intensity, on the other hand, corresponds 


* Eq. (5) of reference 5. 


more closely to the field calculated from the 
dipole.'® It is clear that these effects of the local 
field must be taken into consideration in inter; 
preting the results in terms of the energy spec- 
trum and the composition of the primary radia- 
tion. 


SHOWER INTENSITIES 


In addition to the directional telescopes, the 
apparatus also contains an arrangement of 
counters for recording cosmic-ray showers. These 
counters, indicated in the lower part of Fig. 2, 
were covered with 1.6 cm of lead. The upper 
row of counters were connected together as a 
single unit, those below were alternately con- 
nected to form two units, and coincidences 
between the three units were recorded. During 
the 1936 runs this circuit developed troubles 
early in the season, and no consistent data were 
obtained. The 1935 voyage, however, resulted in 
a fairly consistent set of values for the shower 
intensities, and these are represented in Fig. 9. 
A variation of from 6 to 10 percent is indicated, 
and this is in accord with results reported by 
Johnson" in 1934 that shower intensities show 
less of a latitude effect than the vertical radia- 
tion. The same conclusion has also been reached 
by Pickering.” 
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Fic. 8. Illustration of the correspondence between 
cosmic-ray intensity, I, and horizontal magnetic field, H. 
Curve I, vertical cosmic-ray intensity, percent of high 
latitude value; II, eastern cosmic-ray intensity; III, west- 
ern cosmic-ray intensity; IV, horizontal component of 
magnetic field in gauss along the South American course 
of the Grace Line; V, horizontal component of magnetic 
field calculated from the eccentric dipole by J. Bartels. 


10 J. Bartels, Terr. Mag. and Atm. Elect. 41, 225 (1936). 
4" T. H. Johnson, Phys. Rev. 47, 318 (1935). 
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We are very much indebted to Mr. R. R. 
Adams and other officers of the Grace Line who 
have allowed us to place the apparatus on the 
ships and have cooperated wherever the need 
arose. The expenses incurred in attending the 
apparatus and keeping it in condition have been 
offset by grants from the Carnegie Institution 
of Washington, and funds from the same source 
have been used for assistance in reducing the 
data. For this, acknowledgment is due to 
President John C. Merriam and the Committee 
on Cosmic-Ray Investigations. 
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Intensity of cosmic-ray showers plotted against 
geographic latitude. 
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Effects of Certain Liquids on the OH Vibrational Band of Alcohol 


WALTER Gorpy* 
University of North Carolina, Chapel Hill, North Carolina 


(Received December 26, 1936) 


Transmission curves for several different concentrations of methy! alcohol in dioxane, methy] 
cyanide, isopropyl ether, benzene and carbon tetrachloride, and of ethy] alcohol in dioxane were 
obtained in the region 2.55 u« to 3.15 wu. Dioxane, methyl cvanide, and isopropyl! ether were found 
to shift the OH vibrational alcohol band to the shorter wave-lengths and to increase its intensity 
appreciably. There appeared to be no definite relationship between these changes and the elec- 
tric moment of the solvent molecule. It is suggested that these variations indicate some type 
of interaction of the solvent molecule with the OH group. It may be possible to interpret the 
results as indicating a linkage of the solvent molecules with those of the alcohol through the 


formation of hydrogen bonds. 


_ a recent study of the infrared absorption 
spectra of alcohol-acetone mixtures, marked 
shifts and intensity changes appeared in the 
fundamental vibrational band of the OH alcohol 
group and of the CO group in acetone, while 
other bands were apparently constant. Dipole 
interactions of these groups, possibly through 
the formation of proton bonds, were offered as 
interpretation of these variations. Investigation 
of the OH alcohol group has been continued to 
determine what effects certain other solvents, 
which differ widely in chemical structure and 
in dielectric constants, exert upon this band. 
The resolving instrument and the cell windows 
were the same as those used in the previous 
investigation.' For all the measurements the 
absorbing layer was kept at the constant thick- 


* Now at Mary Hardin-Baylor College. 
Walter Gordy, Phys. Rev. 50, 1151 (1936). 


ness of 0.002 cm. The effective slit width was 
approximately 0.03 

Measurements were made of methyl alcohol 
solutions in dioxane, methyl cyanide, isopropy! 
ether, benzene, and carbon tetrachloride, also of 
ethyl alcohol solutions in dioxane. The alcohol 
concentrations varied from 6 percent to 50 
percent (by volume) in each of the solvents. 
The percentage transmission of these solutions 
for the region 2.55 to 3.154 are shown in Fig. 1. 
The effects of dioxane on the OH alcohol band 
are given in Fig. 1 A. The bottom curve repre- 
sents the transmission of pure methyl alcohol. 
The strong broad band having maximum ab- 
sorption at 2.94 is the fundamental vibrational 
band of the OH group. From a comparison of 
this curve with the curves for the different 
mixtures it will be observed that the influence 
of the dioxane is to shift the OH band to the 
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shorter wave-lengths, and to increase its in- 
tensity. In the curves for 12 percent and 6 
percent alcohol the center of the band appears 
at about 2.75u. This corresponds to a shift of 
0.154. As the alcohol content is increased above 
12 percent, the absorption approaches that for 
pure alcohol. It may be that for these mixtures 
the curves are made up of two broad overlapping 
bands with centers at 2.754 and 2.9u. It was not 
found possible, however, to resolve the band 
into two components. 

To ascertain whether these variations in the 
OH band are general effects of dioxane, solutions 
of the same concentrations of ethyl alcohol in 
dioxane were studied. The results are given in 
Fig. 1 B. From a comparison of Fig. 1 A and 1 B, 
it appears that dioxane produces somewhat the 
same variations in the OH band regardless of 
whether the solute is methyl or ethyl alcohol. 
Likewise methyl cyanide and isopropyl ether 
were found to produce changes in the OH band 
of methyl alcohol similar to the changes observed 
for the OH band in dioxane solutions. These 
results are given in Fig. 1 C and 1 D. Although 
there was some increase in the intensity of the 
band for carbon tetrachloride and benzene, these 
solvents failed to produce the marked shifts and 
intensity changes in the OH band observed in 
the case of alcohol solutions in dioxane, methyl 
cyanide, isopropyl ether, and acetone. Curves 
for carbon tetrachloride-methy! alcohol mixtures 
are given in Fig. 1 E and those for benzene- 
methyl alcohol mixtures in Fig. 1 F. 

It is advantageous to compare the observed 
transmission of the mixtures with that which 
would be expected if the mixtures were purely 
mechanical ones in which neither component 
influences in any way the spectrum of the other. 
Theoretical curves were computed for 6 percent 
and 12 percent alcohol in each of the solvents, 
on the assumption that no interaction occurs 
between the liquids. These, with the experi- 
mentally observed curves for the same mixtures, 
are shown in Fig. 2. The method used in com- 
puting the theoretical transmission of a me- 
chanical mixture is the same as that used in the 
previous paper! and will not be described here. 

Numerical ratios of the band _ intensities, 
experimental to theoretical, with the approxi- 
mate shifts in position are listed in Table I. It is 


assumed that the intensity of the band is 
proportional to the area of the absorption curve. 

The OH vibrational band is much broader and 
appears at longer wave-lengths in liquid alcohol 
than in alcohol in the vapor state.* The broadness 
of the band in the liquid state indicates an 
association between the molecules in which the 
OH group participates. It is possible that the 
molecules in the liquid state are linked through 
the formation of hydrogen bonds.* Working with 
the harmonic of the OH band in the region of 
1.4u, Kinsey and Ellis have reported evidence 
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Fic. 1. Percentage transmission of various alcohol mixtures 
in the region 2.554 to 3.15. 
2 E. L. Kinsey and J. W. Ellis, Phys. Rev. 49, 105 (1936). 
> W. M. Latimer and W. H. Rodebush, J. Am. Chem. 
Soc. 42, 1419 (1920). 
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Fic. 2. Comparison of experimental and theoretical 
transmission of some alcohol mixtures in the region 2.55u 
to 3.15. Solid line, experimental; dashed line, theoretical. 


TABLE I. Ratio of band intensities, experimental and 
theoretical, and shifts in position. 


| | 


PERCENT | APPROXI- INTENSITY 
| OF | MATE RaTIo 
| | experimental 
SOLVENT ALCOHOL (u) | theoretical. 
"Methyl cyanide 6 0.16 | 1.8 
12 O15 | 2.0 
Isopropyl ether 6 0.12 | 2.9 
| 12 O11 | 3.1 
Acetone 6 | O15 2.7 
12 0.15 2.9 
Dioxane | 6 | 28 
| 12 | 0.15 | 2.0 
Benzene 6 0.03 | 1.1 
| 12 0.00 | 
Carbon tetrachloride 6 0.00 | 1.1 
12 0.00 | 1.3 
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for the existence of vapor-like molecules for 
dilute solutions of methyl alcohol in carbon 
tetrachloride.* In the present work, no band was 
observed in the methyl alcohol-carbon tetra- 
chloride solutions which could be attributed to 
such molecules. However, it is possible that the 
band is too sharp to be resolved by the instru- 
ment used. Two of the curves for alcohol in 
benzene show a slight depression on the shorter 
wave-length side of the band which may be 
caused by vapor-like molecules in the solution. 

There appears to be no relation between the 
electric moment or dielectric constant of the 
solvent and the variations produced in the OH 
band. Dioxane is practically nonpolar‘ and has 
a dielectric constant’ of about 2.2. These con- 
stants are very nearly the same as those for 
carbon tetrachloride and those for benzene. 
From these considerations, one would expect the 
influence of these three solvents on the OH 
band to be very nearly the same unless there is 
some type of interaction between the alcohol 
and the solvent molecule. From a comparison of 
the curves, it is evident that this is not the case. 
As has been previously mentioned, dioxane, 
methyl cyanide, acetone, and isopropyl ether, 
all of which differ widely in dielectric constant 
and electric moments, produce very similar 
changes in the OH band. Although the shift is 
to the shorter wave-lengths, the band remains 
broad, and it is not likely that these results 
indicate the production of vapor-like molecules. 
While the association between the alcohol mole- 
cules may be decreased, there appears to be some 
interaction between the alcohol and the mole- 
cules of the solvents. Dioxane, methyl cyanide, 
acetone, and isopropyl ether each contain 
strongly electronegative atoms, oxygen or nitro- 
gen, having unshared electron pairs. It is possible 
that the observed variations are produced by 
the formation of hydrogen or proton bonds 
between the solvent molecules and those of the 
alcohol. 

The writer is grateful to Dr. E. K. Plyler for 
the use of his instrument. 


*C, P. Smith and W. S. Wells, J. Am. Chem. Soc. 53, 


2115 (1931). 
5 J. W. Williams, J. Am. Chem. Soc. 52, 1831 (1930). 
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The Displacement of Principal Series Lines of Rubidium by the Addition of 
Rare Gases 


Ny Tsi-Zé anp Cu'EN SHANG-YI 
Institute of Physics, National Academy of Peiping, Peiping, China 
(Received October 24, 1936) 


The positions of the principal series lines of Rb displaced 
by spectroscopically pure helium, neon and argon are 
measured. The constant shift of the high terms reaches 
37.24 cm™ towards violet under 12.75 atmospheres (at 
570°K) for helium; 2.08 cm™ towards violet under 13.59 
atmospheres (at 590°K) for neon; and 33.4 cm™ towards 
red under 7.12 atmospheres (at 563°K) for argon. The 
effective cross sections for all the atoms in 1 cm’ at 1 mm 


pressure calculated by Fermi's equation is 15.47 cm?/cm® 
for helium, 0.24 cm?/cm® for neon, and 23.7 cm?/cm® for 
argon in good agreement with the results of Fiichthauer 
from his experiment with Na and K, and of Wahlin by 
direct electrical method. In the range of rare gas pressures 
above relative density 5 the mean shifts of the high mem- 
bers of the Rb principal series increases more rapidly than 
linearly with density. 


HE effects of foreign gases upon spectral 

lines, such as displacement, broadening and 
asymmetry, has opened up a new way for 
investigating the perturbations of the energy 
levels of the absorbing atom and consequently 
attracted the attention of many experimental 
and theoretical physicists. V. Weisskopf' and 
H. Margenau and W. \W. Watson? have given a 
general review on the subject with an adequate 
bibliography. In the present research the dis- 
placements of the principal series lines of 

1V. Weisskopf, Physik. Zeits. 34, 1 (1933). 


?H. Margenau and W. W. Watson, Rev. Mod. Phys. 8, 
22 (1936). 


rubidium by the addition of helium, neon, and 
argon are studied. The results for the high series 
members are in agreement with those of other 
experimenters® using sodium and potassium as 
the absorbing atoms, and verify again Fermi's 
theory. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


A nickel steel tube 100 cm long, 2.6 cm inner 
diameter, with two quartz windows, one right- 
3E. Amaldi and E. Segré, Nature 133, 141 (1934); 


Nuovo Cimento 11, 145 (1934); C. Fiichtbauer, P. Schulz, 
and A. F. Brandt, Zeits. f. Physik 90, 403 (1934). 


Fic. 1, Absorption tube. C, water coolings; M, pressure gauge; 7, steel constantan thermo- 
couple; W, quartz windows. 
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Fic. 2. The quartz window. 


handed and one left-handed, on its ends was 
used as an absorption tube. The absorption tube 
was heated by an electrical furnace so wound 
that the central part of the absorption tube 
about 40 cm long was at uniform temperature. 
Water cooling was used near the two ends of the 
tube to prevent the rubidium vapor from being 
deposited on the windows. The general arrange- 
ment of the absorption tube, which could stand 
high vacuum and moderate pressure, is shown 
in Fig. 1. The quartz window W was held by a 
steel frame threaded on to the main tube (Fig. 2). 
Lead rings were employed in the gaskets. The 
steel ring R was to hold the quartz in position 
when the tube was evacuated. 

Near one end of the absorption tube, a side 
tube led to the rare gas tank, the vacuum 
system, and the manometer through three 
separate valves. The manometer ./, read to 
0.01 atmosphere, was carefully calibrated both 
before and at the end of the experiment. 

The construction of the valves is shown in 
Fig. 3. The opening O could be closed or opened 
by a steel tip ¢ which was located in a small 
chamber enclosed by a strong nut N, so that no 
gas could leak out from the outer screw. When 
the valve O was closed, only the gas in the 
enclosure r could leak out very slowly through 
the inner screw; but when O was opened the 
flat part } of the tip will be in firm contact with 
the surface p, thus the leakage was stopped 
immediately. Only a fraction of a turn of the 
handle // would be sufficient to open or close the 
opening O, thus there was no appreciable leakage 
due to the closing and opening of the valve. 
V was an opening for admitting a foreign gas or 
to be connected to the manometer or the pump. 

The temperature of the absorption tube was 
measured by a steel constantan thermocouple. 
The temperature was sufficiently uniform over 
a 40 cm length at the central part of the tube 


AND 


C. SHANG-YI 


Fic. 3. The valve. 


to justify the assumption that within this column 
the concentration of the perturbing gas had a 
constant value. 

The absorption spectrum was photographed 
with a Hilger E-1 quartz spectrograph which 
gives a dispersion of 5.2A to 4.6A per millimeter 
from the 7th member to the end of the rubidium 
principal series. Carbon arc light was used as the 
background of the absorption spectrum, and the 
plate used was Kodak's type 0-IV spectroscopic 
plates. 

The absorption tube was carefully cleaned, and 
the contained gas in the tube wall was removed 
by long pumping and heating. Metallic rubidium 
was inserted from the window into the tube in a 
current of CO,. The tube was then pumped 
again immediately. The absorption spectrum of 
rubidium was taken when the tube was heated to 
a certain temperature. Then the spectroscopically 
pure foreign gas (He, Ne, or A) was introduced 
into the tube, and spectra were taken at con- 
venient steps under different pressures. For the 
determination of the displacement of the low 
members the tube was heated to a lower temper- 
ature (about 220°C) in order to avoid too much 
broadening of the lines. 

An iron arc spectrum was used as a com- 
parison spectrum. The wave-length of the ab- 
sorption lines were measured by Hilger L-1 
comparator. In determining the wave-length of 


the lower members, a Moll type A micro- 
photometer was used. 
RESULTS AND DiIscussION 
Our spectrogram contained all rubidium 


principal lines except the first member. When 
the pressure of a foreign gas is pretty high, the 
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Fic. 4a. Spectra showing the violet shifts of the high members of Rb principal series perturbed by helium. Above: 
without He; below: under 12.75 atmos. (570°K). 


Fic. 4b. Spectra showing the red shifts of the high members of Rb principal series perturbed by argon. Above: without 
A; below: under 5.79 atmos. (566°K). 


Fic. 4c, Spectra showing the slight violet shifts of the high members of the Rb principal series perturbed by neon. 
Above: under 13.59 atmos. (590°K); below: without Ne. 


low members of the absorption lines become very 
broad and the high members are much diffused, 
especially in the case of argon. By the addition 
of helium or neon all the lines shift towards the 
violet, while by the addition of argon they shift 
towards the red (Figs. 4a, b, and c). For helium, 
the shift first increases with the ordinal number 
of the lines in the series, then attains a weak 
maximum at the line 5S-13P, and finally ap- 
proaches a constant value for the lines near the 
limit of the series. For argon, the shift con- 
tinuously increases from the second member and 
approaches also a constant value for the high 
members. The shift for neon is very small and 
practically constant for the high members. The 
results for the three gases under different 
pressures are shown in Figs. 5, 6, and 7. The 
shift of a given line increases with the concen- 
tration of the perturbing gas. 


According to Fermi's theory,’ the displace- 
ment of the high terms of the series is due to the 
superposition of two effects: the polarization of 
the foreign atoms in the field of the core of the 
alkali atom, and the perturbation of the valence 
electron by the foreign atoms. The first corre- 
sponds always to a displacement towards the 
red and is given by 


V.= 200? (€9— 1) n**/8rhcN 

cm", 
where N is the Loschmidt’s number, ¢) the 
dielectric constant of the foreign gas under 
normal conditions, and NV the number of atoms 
per ce. 

The second depends upon the collision cross 
section, ¢, of slow electrons in the gas and is 


4E. Fermi, Nuovo Cimento 11, 157 (1934). 
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Fic. 5. The displacement of the Rb principal series 
lines perturbed by helium. (Violet shifts in cm™ against 
series members.) 


given by 
mc = +1.09-10- na}, 
Tables I-III summarize the constant shifts V 
of the high terms perturbed by He, Ne and A. 


In column 5 are given the calculated values V.. 
The differences between V and V, give the 


values V, in the 6th column of the tables. 

The shifts V, are plotted, in Fig. 8, against 
the relative densities d, the unit being the density 
of the same quantity of a foreign gas as that 
used in the experiment but at 0°C and 1 atmos- 
the 


phere. For relative densities less than 4, 
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Fic. 7. The displacement of the Rb principal series lines 


perturbed by argon. (Red shifts in cm™ against series 


members.) 


relation between the shift V, and the concentra- 
tion is linear which is in agreement with Fermi's 
theory. From the slopes of the straight parts of 
these curves, we deduce the cross sections for 
the single atom of perturbing gases. Multiplying 
these values by N/760, i.e., 3.553X 10", we get 
the effective cross sections for all the atoms in 
1cm* at 1mm pressure and 0°C. For comparison, 
they are listed in Table IV with the results of 
Amaldi and Segré and Fiichtbauer, Schulz, and 
Brandt from Na and K, and those of Wahlin® 
by direct electric method. 

The agreement between the results of the 
present work using Rb as an absorbing gas and 


5H. B. Wahlin, Phys. Rev. 37, 260 (1931). 
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Fic. 8. Shifts y, vs. the relative density of perturbing gas. 


those from Na and K confirms again that the 
shift is independent of the kind of absorbing 
atoms. The deviations of Amaldi and Segré’s 
results for He and A are chiefly due to the 
impurities of the perturbing gases they used. 
(Their helium contained 5 percent of Ne, and 
argon 25 percent of N» as impurity.) 

It is to be noticed that when the relative 
density of the perturbing gas is above 5 the shift 
increases more rapidly than linearly with density. 
This has also been observed by Watson and 
Margenau in their experiment on the pressure 
shifts of the first three doublets of the principal 
series of potassium produced by nitrogen.® The 
departure from linearity of the shifts is to be 
expected as a regular occurrence, if the perturb- 
ing influences are of van der Waals type as 
pointed out by Margenau.’ 
aoa W- Watson and H. Margenau, Phys. Rev. 44, 748 
? Marsenau, Phys. Rev. 48, 755 (1935). 


SPECTRUM LINES . 


TABLE I. The mean shifts of the high terms of the rubidium 
principal series produced by helium. (v, and Vz are values 
calculated from Fermi's equations.) 


PRESSURE, TEMPERA- (RELATIVI v Ve Ve SERIES 
(ATMOS.) TURE Density | (cem™) | (em!) | (cm™!) | MEMBERS 


| (violet); (red) | (violet 


2.11 | 578°K | 1.00 | 5.92] 0.20 | 6.12] 12-32 
3.97 | 569 | 1.91 | 10.21) 0.56 | 10.77, 12-28 
6.24 568 3.01 17.76) 0.89 | 18.65 14-30 
6.25 | 573 17.81 | 0.88 | 18.69 14-29 
8.99 | 574 | 4.26 | 25.14) 1.42 | 26.56! 13-26 
11.39 | 569 5.46 | 32.64| 1.97 | 34.61. 18-27 
12.75 | 570 | 6.10 | 37.24! 2.27 | 39.51| 13-23 


TABLE II. The mean shifts of the high terms of rubidium 
principal series produced by argon. 


PRESSURE| TEMPERA- |RELATIVI v ve Ve SERIES 
(ATMOS. ) TURE Density | | (em™!) | | MEMBERS 
(red) (red) (red) 

1.94 580°K 0.912 9.26, 1.41 7.85 | 15-30 
3.88 571 1.86 17.60 | 3.62 | 13.98) 15-24 
3.92 569 1.88 18.21 | 3.69 | 14.52} 15-24 
5.97 | 566 2.88 29.06 | 6.50 | 22.56) 17-21 
7.00 | 555 3.45 33.65 | 8.25 | 25.40) 15 
7.12 563 3.46 33.40} 8.29 | 25.11} 15-17 


TABLE III. The mean shifts of the high terms of rubidium 
principal series produced by neon. 


PRESSURE| TEMPERA- | RELATIVE v Ve | |. SERIES 

(ATMOS.) | TURE Density | (cm™') | (cm™) | (cm~!) | MEMBERS 


(violet)} (red) | (violet) 


| 
2.63 | 609°K | 1.18 | 0.57 | 0.51 | 1.08 | 16-32 
5.06 | 580 | 2.38 | 0.43 | 1.31 | 1.74 | 15-26 
8.07 | 570 | 3.86 | 0.43 | 2.50 | 2.93 | 15-24 
11.25 | 584 §.26 | 1.52 | 3.75 | 5.27 | 16-24 
13.59 | 590 | 6.29 | 2.08 | 4.77 | 6.85 | 16-26 


TaBLeE IV. The effective cross sections of certain gases 


for electrons of very low velocity as measured by different 


authors. 
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Temperature Shift of the Potassium Resonance Lines* 


Gorpon F. Hutt, Jr. 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 


(Received January 26, 1937) 


Effect of temperature on the shift and half-width of the pothssium resonance lines in 
absorption perturbed by nitrogen pressures from 1 to 15 atmospheres has been investigated. 
Two sets of data, one at 270°C and the other at 90°C, show that both shift and half-width 
vary with temperature in the same manner. As the “relative density”’ of nitrogen is increased, 
the effect of temperature decreases, and becomes unobservable above relative density 7. The 
results are in agreement with Margenau’s theory which predicts the presence of both the 
Lorentz type of broadening dependent on temperature and “‘statistical’’ broadening inde- 
pendent of temperature for low relative densities. The temperature effect decreases with 
increasing relative density, because of statistical broadening. 


INTRODUCTION 


is a recent communication by the author! the 
effects of various pressures of argon and 
nitrogen on the broadening of the potassium 
resonance lines in absorption were reported. It 
was shown that for low pressures the half-width 
and shift of the resonance lines varied linearly 
with the relative density (pressure in atmospheres 
reduced to 0°C) of foreign gas. Above relative 
density 14, the curves of half-width or shift versus 
relative density departed from linearity. These 
results were explained by Margenau’s theory of 
the broadening of spectral lines by foreign gases’ 
which is based on the assumption that the line 
contour is due to a composition of two distribu- 
tions: one, termed velocity distribution in the. 
recent review by Margenau and Watson,’ and 
another, the statistical distribution. The latter 
arises from a statistical arrangement of the per- 
turbers about the absorbing atom and is charac- 
teristic of the specific intermolecular forces; it is 
asymmetrical and independent of temperature. 
This distribution, whose half-width varies with 
the square of the relative density of foreign gas, 
also produces a shift which likewise depends upon 
the square of the relative density. Velocity distri- 
bution, on the other hand, takes account of the 
motion of the molecules (causing incoherence in 
the radiation). It is a symmetrical distribution, 


* Part of a dissertation to be presented to the Faculty 
of the Graduate School of Yale University in candidacy 
for the degree of Doctor of Philosophy. 

1G. F. Hull, Jr., Phys. Rev. 50, 1148 (1936). 

2H. Margenau, Phys. Rev. 48, 755 (1935). 

3H. Margenau and W. W. Watson, Rev. Mod. Phys. 8, 
22 (1936). 


producing no shift, whose half-width Av, is given 
by Lorentz’ formula :* 


Av,;=p*Vn (1) 
or by Weisskopf’s modification of this :*: ® 
Avy =2.2 (2) 


where p is the molecular collision diameter, |” the 
root mean square velocity, } a constant related to 
the intermolecular forces, and m the number of 
perturbing molecules per unit volume. One of the 
consequences of this theory is that for low relative 
densities, the shift Av of a spectral line, produced 
by the cooperation of the two distributions, 
should be approximately related to the half-width 
Av, by: 

(3) 


Since the root mean square velocity varies with 
the square root of the absolute temperature, the 
half-width should, in view of Eqs. (1) and (2), 
depend upon either the one-half or (nearly) one- 
third power of the absolute temperature when 
the relative density of foreign gas is kept con- 
stant, and by Eq. (3) the shift should show a 
similar dependence. 

Few investigators have been careful to dis- 
tinguish between broadening due to the density 
of perturbing molecules and broadening due to 
temperature. It has been customary to increase 
the pressure of the gas by raising the temperature, 
and so to measure shift and half-width as a func- 
tion of both relative density and temperature. 


Av= 


4H. A. Lorentz, Proc. Amst. Acad. 8, 591 (1906). 
5 V. Weisskopf, Zeits. f. Physik 75, 287 (1932). 
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Fic. 1. Wave-lengths of the maxima of the potassium 
resonance lines vs. relative density of nitrogen at 270°C 
(X) and 90°C (O). a, \7699, b, \7665. 


Michelson® was one of the first to establish a 
dependence of half-width on temperature only 
for very low gas pressures. This, however, arose 
from Doppler broadening, which is negligible in 
comparison to velocity and statistical broadening 
for gas pressures above atmospheric. The validity 
of Eq. (1) was demonstrated by Orthmann’ for 
the 2537A line of mercury using mercury vapor 
pressures between 0.1 and 5 mm. Shifts of spec- 
tral lines at different temperatures were investi- 
gated by Holmes* for some of the iron arc lines 
in emission. However, it is difficult to distinguish 
between a shift due to temperature and one due 
to an interatomic Stark effect which arises from 
the high ion density in the arc. Consequently it 
seems more advantageous to work with absorp- 
tion rather than emission lines. 

The present work on the potassium resonance 
lines in absorption was undertaken primarily to 
investigate the effect of temperature only on the 
shifts of the lines, as well as on their half-widths. 
To do this two sets of spectrograms were used, 
one at a high temperature and one at a low, each 
pair having the same relative density of nitrogen. 
Altogether, six spectrograms of the resonance 
lines were obtained at 270°C for various nitrogen 
pressures between 1 and 15 atmospheres. Four of 
these could be compared with the spectrograms 
used in the earlier work,! which were taken at an 
average temperature of 90°C, since the relative 
densities of nitrogen at the high temperature 


* A. A. Michelson, Astrophys. J. 2, 251 (1895). 
*W. Orthmann, Ann. d. Physik 78, 601 (1925). 
*F. T. Holmes, Phys. Rev. 35, 652 (1930). 
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Fic. 2. Half-widths of the potassium resonance lines vs. 
relative density of nitrogen at 270°C (xX) and 90°C (OQ). 
Points are averaged values for the two components of 
the doublet. 


were adjusted to equal those at the low. Ap- 
paratus for obtaining spectrograms at 270°C and 
metheds of measurement and calibration were 
identical with those previously used,' except that 
the thermometers which recorded the tempera- 
ture of the absorption tube were replaced by 
thermocouples and the amount of potassium in 
the tube had to be limited so that there was ap- 
proximately 75 percent absorption at the line 
maxima. The correct amount of potassium to use 
was determined by trial. 


RESULTS 


The results are shown graphically in Figs. 1 
and 2. In Fig. 1 the wave-lengths of the maxima 
of the potassium resonance lines are plotted 
against relative density, and in Fig. 2 the aver- 
aged half-widths of the doublet are plotted. The 
lower curves of each pair, taken from the earlier 
work,! are for 90°C and the upper ones for 270°C. 
Again the curves are straight lines for low relative 
densities, as far as can be determined from the 
position of the experimental points. The effect of 
temperature on shift and half-width decreases 
with increasing relative density, until at about 
relative density 7 it is no longer observable. 

This is entirely in agreement with Margenau's 
theory,’ which predicts the presence of both 
velocity and statistical broadening for low rela- 
tive densities. Although statistical broadening, 
which is independent of temperature, does not 
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become predominant until about relative density 
14,' nevertheless it is always present and becomes 
increasingly effective as the relative density of 
nitrogen is raised. Consequently the high tem- 
perature curve approaches the low temperature 
curve as the relative density increases, rather 
than diverging from it, as would be the case if 
velocity broadening alone were present. Further- 
more from the ratios of the slopes of the curves in 
Fig. 2 to those in Fig. 1, we obtain Ay,;’, Av’ =2.0 
for 270°C and Av; Av=2.5 for 90°C, which is in 
sufficient agreement with Eq. (3). 

Comparison of the results with Eqs. (1) and (2) 
is facilitated by reference to Tables I and IT. In 
Table I the shifts, as reckoned from the line 
origins 7664.84A and 7698.96A,' at 270°C 
(primed) and 90°C (not primed) are listed to- 
gether with the ratios of the shifts, while in Table 
II the averaged half-widths of the two compo- 
nents of the doublet at the two different tem- 


TABLE I, Shifts of the potassium resonance lines perturbed by 
nitrogen at 270°C (primed) and 90°C (not primed). 


RELATIVE DENSITY 


7665A COMPONENT 


0.80 0.360A | O.149A 


2.42 
3.50 0.716 ' 0.615 1.16 
5.89 1.061 1.052 1.01 
8.08 1.358 1.379 0.98 

76994 COMPONENT 

0.80 0.378A O.181A 2.10 
3.50 0.853 0.728 1.17 
5.89 1.366 1.363 1.00 
8.08 1.727 1.739 0.99 
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TABLE II, Averaged half-widths of the two components of 
the potassium resonance lines 7665A and 7699A perturbed 


RELATIVE DENSITY Any’ An Avy’/ Ary 
0.80 2.32cm™ | 0.85 cm™ 2.73 
3.50 4.35 3.63 1.20 
5.89 5.35 3.83 1.39 
04 6.16 0.98 


8.08 6. 


peratures, as well as their ratios, are given. If 
velocity broadening alone were present, the 
ratios AX’ AX and Av,;’/Av, should be constant 
and equal to either the square root or cube root 
of the ratio of the absolute temperatures, 1.22 
or 1.14, respectively. This is seen not to be the 
case, although AX’// AX and Av,;’,; Avy are approxi- 
mately of the right order of magnitude. These 
ratios, however, are in good agreement among 
themselves, which lends further support to 
Eq. (3). 

One may conclude then that the effect of tem- 
perature on the breadth of the potassium reso- 
nance lines is due to the dependence of velocity 
broadening on temperature. As the relative 
density of nitrogen increases, the effect of tem- 
perature becomes less because of the influence of 
statistical broadening. Furthermore the shifts of 
the resonance lines vary with temperature in the 
same manner as the half-widths for constant 
relative density of nitrogen. This is a crucial test 
of Margenau’s theory,”? which these results have 
verified. 

Again the author wishes to thank both Profes- 
sor W. W. Watson, who suggested the problem, 
and Professor H. Margenau, for their advice and 
continued interest throughout the entire in- 
vestigation. 
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On the Behavior of a Certain Supraconducting Circuit 
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The circuit considered consists of a hollow wire of one supraconducting metal and a return 
wire of another supraconducting metal inside the first wire. It is shown that under appropriate 
conditions the current is a function of the temperature. The application of thermodynamics to 
the system results in an equation equivalent to the latent heat relation of Casimir and Gorter. 


CCORDING to present day conceptions the 

current flow in a supraconducting wire 
takes place on the surface of the wire for all 
values of the current not exceeding a certain 
critical value. After this critical value is passed 
resistance returns. Silsbee has proposed the 
hypothesis that this return of resistance results 
from the secondary action of the current’s own 
magnetic field. 

Silsbee’s hypothesis suggests that a different 
behavior might be observed if the current’s 
magnetic field were compensated or shielded in 
some way. For instance, suppose the supra- 
conducting wire were a hollow cylinder and that 
the return circuit were another wire inside the 
first. In this case the configuration of the mag- 
netic field would be quite different from usual, in 
particular the field would be zero at the surface. 
We shall attempt here to give an account of the 
phenomena to be expected in such a circuit, at 
least when the supraconducting metal has certain 
ideal properties. 

The circuit to be considered is shown in Fig. 1. 
We shall imagine the outer cylinder, 8, and the 
inner cylinder, a, both to be supraconducting 
metals. More precisely we wish the inner 
cylinder to be a metal of very high transition 
temperature so that for the range of conditions 
considered below, it shall always be in the 
supraconducting state. 

Suppose the temperature of this system is 
beneath the transition temperature of either 
metal and a weak current is started by means 
of the generator. What path will the current 
take? We may answer this question by using 
the fact that for low temperatures and weak 
currents supraconductors behave as perfect con- 


*Now at Purdue University in the department of 
mathematics. 
5 


ductors. In perfect conductors the current follows 
the path of least inductive energy,'! consequently 
the current will flow along the outer surface of a 
and along the inner surface of 8. 

It would be most remarkable if 8 were to 
remain without resistance for arbitrarily large 
currents, so we shall proceed under the more 
probable alternative that there is a critical 
current strength above which 8 possesses re- 
sistance. Furthermore it is to be expected that 
this critical value is a function of the outer 
diameter of 8. To show this we make the reason- 
able supposition that the behavior of the system 
is independent of whether 8 is a solid piece of 
metal or a series of sleeves. Consider a case in 
which the current is greater than the critical 
value ; then as more sleeves are added there must 
be a stage at which resistance returns. This is 
clear, as one could choose a sleeve of sufficiently 
large diameter such that at no point of the sleeve 
would the magnetic field be comparable with the 
threshold field. Consequently this sleeve could 
carry all the current by itself and there would be 
no resistance. 

One might suppose that the current flows 
along the inner surface of 8 for all values of the 
current strength less than the critical value. 
Suppose that this picture were correct and 
consider a case in which the current is slightly 


Fic. 1. The cylindrical circuit in cross section. 


1M. von Laue, Physik. Zeits. 33, 793 (1932). 
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less than the critical. Then, since the critical 
value depends on the diameter, resistance would 
return if a sufficient number of sleeves were 
removed. By hypothesis the sleeves removed 
were not carrying any current or performing any 
ather apparent function. This would be such a 
remarkable behavior that we may reject this 
picture as improbable. We proceed under the 
alternative that the current takes up an inter- 
mediate configuration between the initial state 
of flow along the inner surface of 8 and the state 
when resistance appears. 

To determine the path of the current in this 
intermediate state we apply Gorter’s condition 
which states that there may be no magnetic 
flux through any supraconducting region. This 
is, of course, an-ideal condition but some actual 
supraconductors behave approximately in this 
fashion. In the first place it is clear that this 
condition requires that any current through the 
supraconducting phase must flow on a surface, 
for otherwise by taking the line integral of the 
magnetic field around a small tube of current 
we should find a non-zero magnetic flux. Since 
there is no dissipation, the current must not 
flow through the normal phase. A conceivable 
path would be for the current to flow on a 
cylindrical surface of radius p. Inside this radius, 
8 must be in the normal phase for the magnetic 
flux is there different from zero. Outside this 
radius the magnetic field is zero so the material 
is in the supraconducting phase. The writer sees 
no other possibility consistent with Gorter’s con- 
dition and with the symmetry of the system. 

To determine the radius, p, we apply Gorter 
and Casimir’s rule? which states that at the 
point of contact of the normal and supracon- 
ducting phases the magnetic field equals the 
threshold field, HZ; Consequently if 7 is the 
current strength 

H,=2i1/p. (1) 


In what follows we shall only be interested in 
the intermediate state when p>a. 
The work done on the system by the generator 


in the time dt is 
dW = —iedt, 
where ¢ is the counter e.m.f. Since the current 


2C, J. Gorter and H. Casimir, Physica 1, 305 (1934). 
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may always flow along a resistanceless path 
e=—dN/dt where N is the flux linked with the 
circuit. Hence 
dW=idN. (2) 
The flux may be calculated from the in- 
ductance, L, since N=Li. For a cylindrical 
inductance in the notation of Fig. 1 we have 
L=21 log p/a 
and by Eq. (1). 
N =12l log (3) 


We wish to consider now a case in which no 
work is done on or by the system but in which 
the temperature is varied. Since dN =0 Eq. (3) 
gives 

i log (4) 
But //; is a function of the temperature so we 
obtain the interesting consequence that the 
current depends on the temperature. 

The first law of thermodynamics applied to 
the system is by Eq. (2). 


dE=dQ+idN, (5) 


where dQ is heat added and dE is the increase 
of internal energy. Since we are considering an 
ideal case with complete reversibility, the second 
law may be applied. We do this by requiring that 
dQ/T be a perfect differential, dS. Thus Eq. (5) 
becomes 


dE=TdS+idN. 
This may be rewritten in the form 
d(TS+Ni—E) =SdT+ Ndi. 
The left side of this equation is obviously a 


perfect differential, hence the right side is also 
and we may equate the cross derivatives. 


or 
If we now make use of the “equation of state,” 
Eq. (3), we obtain 
(6) 
The left member of Eq. (6) is to be interpreted 
as the amount of heat which must be added to 


maintain the temperature constant when the 
current increases one unit. 


to 


he 


HARTREE AND HARTREE-FOCK METHODS 


To gain a clearer conception of the significance 
of Eq. (6) we shall transform it by a change of 
variable. The volume of a cylinder of radius p 
and height / is 

By using Eq. (1) 
Then at constant temperature 
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Inserting this expression in Eq. (6) we obtain 
(dQ/dV)r= —T(H,/ 42) dT). (7) 


It is clear that the left member of this equation 
is to be interpreted as the amount of heat which 
must be added to maintain temperature constant 
when a unit volume transforms from the supra- 
conducting to the normal phase. Eq. (7) is 
therefore identical with the latent heat equation 
of Gorter and Casimir.* 


REVIEW VOLUME 51 


Note on the Hartree and Hartree-Fock Methods 
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Certain limitations on the Hartree and Fock methods of 
handling many-body problems in quantum mechanics are 
discussed with the help of simple examples. One may treat 
the interaction in such problems as composed of two parts: 
first, a part which is correctly accounted for by the above 
methods, and secondly, a part which is entirely omitted 
in these approximations. If the latter type of interaction is 
not negligible one may be led to two types of solution of a 


1 


T is commonly stated that the Hartree method 
of obtaining an approximation to the solution 
of a many-body problem in quantum mechanics 
yields the best possible approximation subject 
to the restriction that the wave function of the 
system be expressed as a product of one-particle 
wave functions. In the extension of this method 
by Fock and Slater to include exchange effects, 
the restriction on the form of the wave function 
is that it be an antisymmetric linear combination 
of products of one-particle wave functions. Such 
statements are apt to imply a false notion as to 
the accuracy of the results obtained by these 
methods and in this note we propose to analyze 
a limitation of these methods which in certain 
cases may be of a serious nature. 

For the sake of simplicity, let us consider the 
problem of two coupled oscillators in one 
dimension. If x and y denote the coordinates of 
the two particles, m the mass, w) the uncoupled 
angular frequency of each particle and \ a 
coupling coefficient, the Hamiltonian of the 


problem; first, a complete set of orthogonal Hartree wave 
functions, and secondly, additional wave functions of 
entirely different nature from the Hartree functions. These 
two types of solution are not independent, however, and 
it is shown how to establish connections between them 
with the help of perturbation methods. The importance of 
these limitations in connection with the theory of nuclear 
structure is pointed out. 


system can be written as 


2m dx? 2 (1) 
h? 0? 

2m 2 


This problem can be solved exactly by intro- 
ducing normal coordinates (3)'(x+y) and 
(3)'(x—y). The normalized wave function for 
the ground state is 


Yo(x, ¥) = (~) (wow)! 
th 


wo | 
2 2 


and the energy of the system in this state is 
tha; (2a) 
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If we apply the Hartree method to this 
problem, setting Yo=u(x)u(y), we readily find 
that the ~’s must satisfy the equation 


+ 2mm u*(é)Eu(E)dE 


400 


with u*(y)TTeu(y)dy. 


In Eq. (3) the first integral vanishes, since 
u*()u(&) is an even function of & This integral 
arises from the term 2mAxy in the Hamiltonian 
and we thus see that the Hartree procedure fails 
to take this term into account. The Hartree 
wave function for the ground state is 


mo'\ mo’ 
Yo(x, vy) = (= ) exp ——(x?+y? 
wh 2h r (4) 


and the energy is Eo = hw’ jw" = wo? 


For purposes of comparison of the solutions 
(2) and (4) it is advisable to consider the case of 
small coupling and expand the results in a 
power series in A. The energies agree only to 
terms in the first power of \ and the wave 
functions do not agree even to the first-order 
terms. More explicitly the exact energy is 
+--+) and the Hartree 
energy is E=hwo(1+X/wo?—A*/wo'+---). For 
the wave functions of the ground state we find, 
keeping terms in the first power of \, 


th 


| mr Mao 
xX} y)? exp| 
| 2h 


May 
(Hartree) yYoo= (= ( +.) 
| mx 


| Mw o 
(x?+ y?) exp| -— 
dh 


x}1-— 
| 


The difference in the wave functions clearly 
arises from the fact which we have already 
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mentioned that the Hartree method does not 
recognize the existence of the interaction energy 
term 2mAxy. It is easy to see from a physical 
viewpoint why the Hartree method is inadequate 
to handle the above problem correctly. The 
motion of two uncoupled oscillators may be 
described in terms of the motion of their center 
of mass (essentially the normal coordinate x+y) 
and of the motion relative to the center of mass 
(essentially the other normal coordinate x—y). 
The interaction affects only the relative motion, 
the motion of the center of mass being entirely 
unaffected by the presence of the coupling forces. 
Thus we would expect that the coupling would 
alter the frequency of the relative motion, 
leaving the frequency of the motion of the center 
of mass unchanged. As we see from Eq. (2), 
this is exactly what the correct solution of the 
problem yields. On the other hand, it is char- 
acteristic of the Hartree method that the co- 
ordinates of the various particles are treated in 
an identical manner. Consequently we find both 
the motion of the center of mass and the relative 
motion equally perturbed in this method, and in 
such a manner that the perturbed frequencies 
give the correct energy to the first order of the 
coupling coefficient \. 

The example which we have cited shows that 
we must recognize the possibility that a part of 
the interaction energy in a many-body problem 
may be entirely neglected in the Hartree approx- 
imation. In general we may think of the inter- 
action as being composed of two parts, the first 
of which is correctly accounted for by the 
Hartree procedure, and the second part the effect 
of which is entirely omitted in the solution. 

In the following it will be convenient to limit 
our discussion to problems involving only the 
second type of interaction so as to avoid the 
necessity of performing a preliminary Hartree 
calculation. This imposes no restriction on the 
generality of the results, since we may always 
perform a Hartree calculation and utilize the 
wave functions and energies so obtained as the 
starting point of a new calculation. Thus in the 
case of the coupled oscillators we may equally 
well consider the problem described by the 
Hamiltonian 


(5) 


f 
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a form very often employed and entirely equiva- 
lent to (1) if we imagine the unperturbed 
frequencies wo of (1) to be replaced by w’. In the 
case of the Hamiltonian (5), the solutions of the 
problem in which we omit the coupling term 
2mdxy are the exact Hartree solutions. It is not 
difficult to see that in any problem (still re- 
stricting ourselves to the case of one dimensional 
motion of the individual particles) in which the 
potential energy of the system in the uncoupled 
state is an even function of the particle co- 
ordinates, any interaction which is an odd 
function of these coordinates cannot be handled 
by the Hartree method. This fact is exceedingly 
evident in the case of the problem under con- 
sideration. The correction to the potential field 
in which one of the particles moves is obtained 
by averaging the interaction energy over the 
motion of the other particle. Since the interaction 
is proportional to the coordinate of the second 
particle, this correction is proportional to the 
latter's average position, which is zero. Further- 
more, this correction to the potential is inti- 
mately related to the first-order perturbation 
energy when the problem is treated by standard 
perturbation theory. The latter is obtained by 
averaging the perturbation energy over the 
unperturbed wave functions, and it is just this 
average over all the partial wave functions with 
the exception of that for the particle under 
consideration which constitutes the Hartree 
correction to the potential. Thus, whenever a 
standard perturbation calculation yields no cor- 
rection to the energy in first approximation, a 
Hartree calculation is totally inadequate to 
handle this sort of interaction. This is in accord 
with the result obtained by Moller and Plesset,' 
who have shown that if the Hartree-Fock 
solutions to a problem are used as the starting 
point of a perturbation calculation, the first 
order correction to the energy vanishes. 

It should be mentioned that the Fock scheme 
does vield nontrivial results when applied to the 
problem described by (5), of course, only for 
states other than the lowest. It is interesting 
to note that here we have a case where the 
exchange corrections alone are different from 
zero, instead of being small compared to the 
Hartree terms as in the atomic case. Except for 


‘ Moller and Plesset, Phys. Rev. 46, 618 (1934). 


the first excited state, this particular example is 
not a particularly simple one to discuss from the 
standpoint of the Hartree-Fock method due to 
the presence of degeneracies apart from exchange 
degeneracy which must be removed before the 
method is applicable. Nevertheless, it is clear 
that in more complicated problems the same 
sort of restrictions which we have demonstrated 
for the ordinary Hartree method may well exist 
for the exchange terms in the Fock equations. 


2 


The limitations which we have discussed in 
the preceding section may result in a situation 
in which the Hartree or Hartree-Fock procedures 
converge to results which are widely different 
from the correct solutions of a many-body 
problem. In particular, the ground state of a 
system may possess an energy much lower than 
that yielded by the above methods and the wave 
functions may be of entirely different character 
than the Hartree-Fock functions. In this section 
we shall illustrate this by a simple example. 

Let us again focus our attention on the 
problem of two particles, each performing recti- 
linear motion, and imagine that we have con- 
structed a surface representing the potential 
energy of the system, using x and y as before as 
particle coordinates, plotting the potential energy 
vertically as ordinate with the «—y horizontal 
plane representing configuration space for the 
system. We shall suppose that for uncoupled 
particles there is a single minimum in this surface 
which we place at the origin of coordinates, and 
that the surface approaches positive infinity as 
we move away from the origin in order to 
restrict our considerations to the case of a discrete 
eigenwert spectrum. Thus for the case of a pair 
of uncoupled harmonic oscillators the potential 
energy surface is a paraboloid of revolution, and 
the wave function for the ground state is a 
two-dimensional Gauss error function with its 
maximum at x=y=0. Let us now imagine an 
interaction energy present which is an odd 
function of x and y and which so modifies the 
potential energy of the system that new minima 
appear in the surface. If these minima are far 
from the origin and deep enough so that the 
minimum potential energy is large and negative, 
we are certain that the ground state now 
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possesses a wave function with maxima at the 
positions of these new minima in the energy 
surface and that the energy of the state may 
lie much below the unperturbed energy. This 
means that the interaction has produced new 
configurations of stable equilibrium. 

As a specific example, let us consider the case 
of two harmonic oscillators with a coupling 
energy of the form —kxye~**+”), Although this 
example is of little physical interest, it will serve 
to illustrate the essentials of the problem and 
possesses the advantage of avoiding unnecessary 
mathematical complexity in the later analysis. 
The Schridinger equation for this problem is 


2m\dx* 
Mwy" 


+—(x?+-y*) 
? 


=0. 


If we use the quantity (#/ma)! as the unit of 
length and fiw/2 as the unit of energy, the 
above equation becomes 
ax? oy? 
with a=k/ maw and B=hc/ma. In this system 
of units the potential energy becomes 


V(x, y) =x? + y? — +0"), (7) 


One easily verifies that this function possesses 
minima along the line x=y at the positions 
Xo = Yo= +(20/28)!, where 2 is the solution of the 
transcendental equation 


a(1—z) 


The point x=y=0, which was a minimum for 
the surface =.°+y", now becomes a saddle 
point and no longer represents a configuration of 
stable equilibrium. The transcendental equation 
possesses real positive solutions only for a>1, 
and furthermore 0=2~=1. The value of the 
potential energy at the points x»= yo turns out 
to be 
Vo= —20°/B(1—20), 


which is negative for all values of 2) and positive 8. 
We may now obtain a good approximation to 
the solution for the ground state of the system 
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by expanding the potential energy in a power 
series about the minima and breaking off the 
expansion with quadratic terms. The result is 


1+25 
+ (a+ ~) (Lx 


-0 


—2(1—29)(x—xo)(y—yo) (8) 


and we can solve the problem of the motion of 
the system with this potential energy with the 
help of normal coordinates. Introducing the 
abbreviations and v?= 
2 (1—20), and denoting the normal coordinates by 
(3)'L(x—x0) ] and 
—(v—yo) ], the wave function for the ground 
state has the form 


¥=const. a2? (9) 
and the energy of this state is 
E=Votut». (9a) 


The complete wave function will consist of a 
symmetrical or antisymmetrical linear combina- 
tion of two terms, each of the form (9),—one 
corresponding to and the 
other to x» —(20/28)!. If the points 
lie far enough apart so that the two wave 
functions do not overlap appreciably, the energy 
will be given very nearly by (9a). If, however, 
there is considerable overlapping, the energy 
level (9a) will be split into a singlet and triplet 
due to spin degeneracy, and the splitting can be 
calculated by the ordinary Heitler-London 
method. 

In contrast with the above results, the 
Hartree-Fock wave function for the ground state 
is simply 

y) =u(x)u"(y), 


where the «°’s are the ground state solutions of 
the harmonic oscillator problem. Each of the 
u°’s satisfies an equation of the form 


@u°/dx?+ (e—x?)u°=0 


and the ground state energy is +2 units in this 
method. Thus we have a situation in which we 
are in possession of a complete orthogonal set of 


V(x, y)=-— 

B(1—2zo) 
t 
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Hartree solutions to a problem and in addition 
other solutions. Care must be taken in such cases 
not to treat the additional solutions as inde- 
pendent of the original set, or conversely if the 
latter type of solutions are utilized, one must not 
add to them all the original Hartree solutions. 
As we shall show in the next section, some of the 
original Hartree solutions go over into the second 
type when a perturbation calculation is per- 
formed which takes into account the interaction 
which is omitted in the Hartree scheme. One 
meets such a case as this in the theory of metals. 


3 


In the preceding section we have demonstrated 
the existence of wave functions for a system 
which are markedly different from the Hartree- 
Fock wave functions. It is of importance to 
investigate the relations between these solutions 
and to see how the latter solutions go over into 
the former when the interaction is taken into 
account. This we shall do by using the standard 
perturbation theory and we shall see, at least in 
a favorable case, that the result of such a 
perturbation calculation actually modifies the 
Hartree-Fock solution in such a way as to yield 
a result which bears a great resemblance to the 
wave function obtained utilizing the expression 
(8) for the potential energy of the system. It 
turns out that there is interference among the 
contributions of the excited unperturbed states 
to the perturbed wave function of such a nature 
that new maxima appear in the latter function 
and presumably if the perturbation calculation 
were continued sufficiently, we would be led to 
the correct ground state wave function. 

Using a potential energy of the form (7) and 
treating the last term in this energy as a per- 
turbation, we have as unperturbed wave func- 
tions : 
= 


where ,° and x;° are solutions of the problem 
of the harmonic oscillator. For the ground state, 
the first-order perturbed wave function is given 
by 


¥) =Yoo"(x, y) 


00 


Ex Evo 
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where the term k=/=0 is to be omitted from 
the summation and 


wom f  Yoo'dxdy. (10a) 


The matrix elements /7;:, 99 may be written in 
the form Wi, with 


Min= f 827 x)dx, (10b) 


so that the perturbed wave function (10) can 
be written as 


Yoo(x, ¥) = Wo0"(x, 


Mio: Mio 
y). (11) 
kik 


+00 


From our previous discussion we see that the 
first-order correction to the energy vanishes (it 
is equal to Hoo,0=0) and the second-order 
correction is 
kl Eo — Ex: 


Thus the calculation consists essentially of the 
evaluation of the matrix elements ./. In our 
system of units the denominators are simply 
given by 

2(k+1). 
The matrix elements may be evaluated in closed 
form by making use of the generating function 


for the Hermite polynomials. If we consider the 
two generating functions 


F(x, = ——-o* 
k! 
I1)(x) 
and G(x, r) = r! 


+-co 
and write "dx 


we can evaluate this integral and expand the 
result in a power series in o and +. Comparison 
of the terms in this power series with those 
obtained by replacing F and G by the series 
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Fic. 1. Cross sections of the potential energy surface and of the wave functions along the line 
x=y. I, potential energy; Il, wave function from parabolic potential; III, wave function from 
first-order perturbation calculation; A, unperturbed ground state energy; B, ground state energy 
from second-order perturbation calculation; C, ground state energy from parabolic potential. 


expressions in terms of the Hermite polynomials 
then leads to the values of the matrix elements 
My. The only elements Wj which are different 
from zero are those in which & is odd and these 
can be written in the form 

(—1)"7 (2n4+1)!} 


J2n+1 


1 B\" 
x--— ( ). (13) 


In order to compare the results of the perturba- 
tion calculation with those expressed by Eqs. (9) 
and (9a), it is necessary to insert numerical 
values for the parameters a and 8. Furthermore 
it is desirable to choose these values in such a 
manner that the first terms of the perturbation 
calculation as expressed by Eqs. (11) and (12) 
yield sufficiently large corrections to the unper- 
turbed functions and energy to clearly indicate 
the trend of the method. This can be accom- 
plished by utilizing such values as place the 
minima in the potential energy surface suffi- 


n! 


ciently close to the origin, so that the unper- 
turbed wave function and the perturbation 
energy are both appreciably different from zero 
in the same region of configuration space. The 
numerical values which have been used are as 
follows: 2=0.60, 8=0.20, a turns out to be 
equal to 4.55, the minima of the potential energy 
surface lie at xp=V¥o=+1.225 units and the 
minimum value is —4.5 units. The energy of 
the ground state as calculated from (9a) is 
—0.21 unit. 

From the perturbation calculation it turns out 
that the second-order correction to the energy of 
the ground state is —1.80 units, giving a per- 
turbed energy of +0.20 unit which is not very 
different from the value given by (9a). The 
results of the calculation are shown in Fig. 1, in 
which are plotted cross sections of the potential 
energy surface and of the wave functions along 
the line x=y. The curve I is a plot of the po- 
tential energy, curve II the symmetrical wave 
function of the ground state as calculated from 
Eq. (9), and curve III the perturbed wave 
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function as calculated from Eq. (11). These 
wave functions are normalized. The energy 
levels A, B, C indicate the unperturbed energy 
of the system, the second order perturbed energy 
and the ground state energy calculated from the 
parabolic potential. 

It is evident from the figure that the approxi- 
mation of a parabolic potential energy leads to 
too high an energy for the ground state, since 
the period of the motion is actually larger than 
that which we obtain from the above approxi- 
mation. The effect of this will not only be to 
lower the level C but also to alter the wave 
function II so that the maxima will lie nearer 
the origin than shown and the minimum at the 
origin will be less pronounced in the correct 
wave function. Keeping this in mind, the simi- 
larity between the curves II and III is note- 
worthy and it seems certain that a continuation 
of the perturbation calculation to higher order 
terms would yield a function which would 
approach the correct one. The Hartree-Fock 
function is not plotted but is simply a Gauss 
error curve with Hs maximum at the origin. 

It must be pointed out that in general the 
perturbation method may not be practicable as 
a means of getting at the correct solutions of a 
many body problem from the Hartree solutions. 
In fact other calculations than the one given 
above indicate that unless both the Hartree 
functions and the perturbation energy are both 
fairly large in the same region of configuration 
space, the perturbation method will converge 
very slowly indeed, and a prohibitive amount of 
labor would be necessary to obtain a reasonable 
degree of accuracy. At this point we should like 


to emphasize the need of better methods for 
dealing with such situations. More general wave 
functions than those composed of products of 
one particle functions must be employed. Per- 
haps some generalization along the lines sug- 
gested by normal coordinates, i.e., wave func- 
tions built up of products of partial wave 
functions each of which depends on_ linear 
combinations of the individual particle coordi- 
nates, may lead to the desired result. 

The limitations of the Hartree method have 
been recently emphasized in connection with the 
theory of nuclear structure. Bohr? especially has 
insisted on the inadequacy of the concept of 
configurations, i.e., one particle wave functions, 
for the purposes of describing the structure of 
the heavier nuclei. Bethe and Bacher*® have also 
discussed the weakness of the method, and 
Elsasser* has briefly discussed this question in 
connection with the calculation of the density 
of nuclear levels. It seems to the writer that not 
only the fact that the energies of different con- 
figurations overlap considerably in the heavier 
nuclei, but also the limitations on the Hartree- 
Fock method stressed in this note are of im- 
portance in connection with the inadequacy of 
the description of nuclei by means of one- 
particle wave functions. 

In conclusion, it is a distinct pleasure to 
express my sincere thanks to Professor E. Wigner 
both for suggesting the problem presented in 
this note and also for many valuable discussions. 


N. Bohr, Nature 137, 344 (1936). 
’ Bethe and Bacher, Rev. Mod. Phys. 8, 150 (1936). 
*W. M. Elsasser, Phys. Rev. 51, 55 (1937). 
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T. BAINBRIDGE, in 1929, reported a 
e “search for element 87 by analysis of posi- 
tive rays.”"' Although positive ray analysis is un- 
doubtedly the most sensitive method for the de- 
tection of the presence of a new element, Bain- 
bridge, studying likely sources (lepidolite and 
pollucite, minerals rich in the alkali metals) failed 
to detect in them the presence of element 87. 

Since then, others have searched for this new 
element. F. Allison and E. J. Murphy?’ reported 
the presence of element 87 in lepidolite and 
pollucite. Their method is open to question. 

J. Papish and E. Wainer* reported the finding 
of several of the stronger x-ray lines of element 
87, using a concentrate derived from samarskite. 
The validity of their findings will be discussed 
later in this paper. 

Recently, Hulubei,* using a curved crystal 
focusing x-ray spectrograph of the transmission 
type, reported weak lines at 1032 and 1043 X.U., 
which he assumed were the La, 2 doublet of 
element 87; from which he concluded that his 
pollucite extract contained element 87. Fig. 1 
shows a Moseley plot ((v//R)! vs. atomic number) 
of data (indicated by circles) for the Lay, » 
doublets of various elements as secured from 
Siegbahn’s Spektroskopie der Roentgenstrahlen. 
Also plotted, for atomic number 87, are values 
of (v/R)' calculated from Hulubei’s data (indi- 
cated by crosses). 

The Siegbahn data (indicated by circles) follow 
the Moseley law accurately (to less than 0.01 
X.U.), but at Z=87 the crosses calculated from 
Hulubei’s data are in error by 6 and 7 X.U. for 
La, and Lag, respectively. This may be due to a 
consistent error in the wave-length scale used by 
Hulubei. He remarks in his paper that his x-ray 
tube target is contaminated by mercury, and it 
does not seem unreasonable to suppose that these 
(so designated) lines of element 87 might actually 


1K. T. Bainbridge, Phys. Rev. 34, 752 (1929). 

? F. Allison and E. J. Murphy, Phys. Rev. 35, 285 (1930). 

3]. Papish and E. Wainer, J. Am. Chem. Soc. 53, 3818 
(1931). 

*H. Hulubei, Comptes rendus 202, 1927 (1936). 


be the mercury lines L8; and L8@; (A=1030 and 
1047 X.U., respectively). Hulubei’s measured 
wave-lengths for 87 differ from these mercury 
wave-lengths by less than his error in determining 
the wave-lengths of 87 Lay, ». 

Recently the present author undertook a 
search for x-ray lines of element 87, at the request 
of Dr. T. G. Kennard of Claremont Colleges, 
California. The x-ray target material, very kindly 
supplied by Dr. Kennard, was CsHSO, derived 
from lepidolite; the x-ray tube was operated at 
30 kv d.c. and 8-10 ma (the Z excitation potential 
for element 87 was determined to be 19 kv). A 
Moseley graph was used to determine the Bragg 
angle. The rocking device,® used to oscillate the 
crystal in the spectrograph, was set for minimum 
amplitude of oscillation. 

An eight-hour plate was taken at the proper 
angle to register the La lines of element 87, the 
target material being CsHSO, concentrated for 
element 87 (as supplied by Dr. Kennard) 
pounded into a roughened copper wedge. A faint 
narrow band or “‘line’’ appeared at the center of 
the crystal reflection zone at the position calcu- 
lated for the La doublet of element 87. Fig. 2 (a) 
shows a 10-hr. spectrogram taken with CsHSO, 
more strongly concentrated for element 87, used as 
target ; the concentration of element 87 here was 
estimated by Dr. Kennard at about six times the 
concentration for the previous sample. The band 
or “‘line’’ which appeared on the first plate seemed 
to have increased roughly sixfold in intensity. 
Such an increase in intensity with increase in 
concentration for element 87 would seem to be 
suggestive. 

However, the band or “‘line’’ under discussion 
here is broader than an L x-ray line such as Th 
La;,». Accordingly, an investigation was made 
to determine whether the Th (90) La lines (un- 
resolved by this spectrograph, and very nearly in 

5F. K. Richtmyer, Phys. Rev. 37, 472 (1931). 

® At this time, the concentrate which was assumed to 
contain element 87 in strong concentration, was tested for 
a, 8 and vy radioactivity with the result that none was 


detected within experimental error (5 percent of the total 
number of counts). 


584 


‘| 
| 


SEARCH FOR 


the same spectral region) would broaden as the 
concentration of ThOs in the CsHSO, target 
material was decreased. No such broadening 
could be detected in concentrations as low as 
(0.125 percent of ThO, in CsHSO,. A 10-hr. plate 
was made using a blank copper wedge as target. 
The band or “‘line”’ again appeared at the position 
where the La doublet of element 87 should be 
(see Fig. 2 (b)). It was then suspected that this 
band or “‘line’ might be due to a crystal imper- 
fection located at the center of the calcite crystal : 
an abnormally strong reflection of Cu continuous 
radiation occurs for a narrow vertical zone of the 
crystal surface. Moreover, with close scrutiny, in 
Figs. 2 (a) and 2 (b) the band or “‘line”’ is ob- 
served to have the same peculiarities, viz.: it is 
split (double) at the top and center. This would 
point to the band or line as undoubtedly due to 
crystal imperfections which were accidently per- 
mitted to register photographically. Play in the 
rocking device, coupled with small amplitude of 
oscillation (rocking) would keep the crystal from 
moving appreciably. It is known that crystal re- 
flection defects can be eliminated by rocking a 
crystal which is reflecting x-rays since the locus, 
on the crystal surface, of reflection of a given 
wave-length, then moves back and forth along 
the length of the crystal thus eliminating the 
crystal defects which are “ironed out.”’ With a 
view to the possibility of eliminating the band or 
line in this manner, the amplitude of oscillation 
of the calcite crystal was increased to about 45’; 
the ‘‘line’’ vanished (see Fig. 2 (c), a 10-hr. plate) 
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Fic. 1. Moseley graph for the x-ray lines Lay»: circles 
are from data of Siegbahn’s Spektroskopie der Roentgen- 
strahlen; crosses are from data given by Hulubei.! 
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Fic. 2. (a) 10-hr. spectrogram taken using lepidolite- 
derived CsHSO,, strongly concentrated for element 87. 
Note “‘line’’ at center. (b) 10-hr. spectrogram taken using 
blank Cu wedge, showing same “‘line”’ as in Fig. 2. (a) With 
identical peculiarities. (c) 10-hr. spectrogram taken with 
blank Cu wedge as target, with increased amplitude of 
crystal oscillation. Note ‘‘line” has vanished. 


showing conclusively that it was due to a crystal 
surface defect. 

An important point of this study is that the 
calcite crystal used by the present author is the same 
calcite crystal which was used by Papish and 
Wainer.’ They make no mention in their article of 
oscillating the crystal to eliminate the photo- 
graphic registration of its surface defects ;’? hence 
the possibility cannot be ruled out that their 
reported “‘lines’’ may be due to the same crystal 
defect as reported in this note. Without crystal 
oscillation, faint lines would be masked by the 
crystal defect “‘line.”” Since the zone of intense 
reflection occurs at the exact center of the crystal 
(see Fig. 2 (a)), and since the crystal is cut to 
just match the crystal holder in length, the 
“‘line,”’ which is due to this crystal defect, will 
appear at the position where the element 87 line 
should appear irrespective of the magnitude of 


7 If the crystal rocking device had been set for a small 
amplitude of oscillation, say, for example, 45’, the pla 
(~25’) in the crystal rocking device on the spectrograp 
which Papish and Wainer (reference 3), used would be 
sufficient to keep the crystal from oscillating. If the crystal 
rocker had been set for a large amplitude, for example, 
120’, the heavy background (see Fig. 2 (c), zone A—>A’), 
resulting from the ‘‘smearing”’ out on the plate of the false 
‘line’ (see Fig. 2 (a)) would very easily mask faint lines of 
element 87 as the crystal would only ‘‘reflect’’ the line of 
element 87 for half of its excursion. This can easily be seen 
from a simple geometrical construction, as the locus of 
“reflection” of the line at 9° 45’, the Bragg angle for element 
87La, passes beyond the ends of the crystal for each final 
quarter of its travel across the crystal face. (The spectro- 
graph radius is 18.36 cm; the crystal is 2.5 cm long.) In 
any event, in order to have true lines of element 87 appear 
on their plates, Papish and Wainer must have rocked their 
crystal through a large amplitude and must have had 
element 87 present in considerable concentration. An 
examination of Papish and Wainer’s plates would clear up 
the whole question; however, they do not seem to be 
available to me in Ithaca or elsewhere. 
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the Bragg angle at which the element 87 line is 
calculated to be found. 

This crystal which Papish and Wainer* used, 
was temporarily transferred to the spectrograph 
which the writer is using,* as the other specially 
cut calcite crystal was not available at the time. 

The author feels that the search for element 
87 is still open. The writer failed to detect the 


* It was not used in research leading to other publications 
of the present author. 
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presence or absence of the Lay, » lines of element 
87, due mainly to the lack of a sufficient amount 
of concentrate to permit a conclusive search. 

The author wishes to thank Dean F. Kk. Richt- 
myer for permission to use certain x-ray ap- 
paratus granted to Dean Richtmyer by the 
Heckscher Research Council of Cornell Univer- 
sity. The writer also wishes to thank Dean Richt- 
myer, Professor Kennard and Dr. Parratt for 
many helpful discussions during the progress of 
this study. 
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A Note on the Density of Eigenfunctions for an Electron Obeying Dirac’s Equation 


E. K. Brocu 


University of Oslo,* Norway 


(Received January 14, 1937) 


N dealing with the theory of an ideal gas by 

the method of statistical mechanics, one 
requires to know the number of independent 
eigenfunctions in a given energy range. As is 
well known. in the nonrelativistic mechanics, the 
required enumeration of the states can be carried 
out in two ways. We may enclose the electrons 
in a box, by making the wave functions vanish 
on a certain boundary, or we may impose 
periodicity conditions on the wave function, the 
resulting density of eigenfunctions being the 
same in each case. For an electron obeying 
Dirac’s equation the number of eigenfunctions is 
easily found by the latter method.' However, 
such periodic boundary conditions for the wave 
function can never in practice be realized. In 
the following note we shall therefore show, that 
for electrons enclosed in a spherical potential 
hole, the two methods lead to the same result. 
This calculation is also of special interest, 
because it recently has been questioned? whether 
it is justifiable to use the result proved for 
progressive waves, when the physical problem 
might require the use of “standing waves.” In 
the special case examined, namely spherical 


* Submitted from Cambridge University, England. 

' Moller and Chandrasekhar, M. N. Roy. Astr. Soc. 95, 
673 (1935). 

2? Eddington, Proc. Roy. So« 


. 258 (1935). 
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enclosure, ‘‘standing waves” does not lead to a 
new physical situation. 
The wave equation (J7—J/')y=0, written in 


terms of the components, gives the equations 


Il’—V Ove h 

c or r 

{- (1 ) 

Oy, | 
po th——j 

or r 

Putting 


II'—V W-v 
C 


(1) reduces to 


Eqs. (2) are of standard type, the solutions 
involving Bessel functions in the general sense. 

We have to choose solutions which satisfy the 
given physical conditions. Let us enclose the 
electron in a spherical box by assuming po- 
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tentials =0 for r<R and V=U for r>R. 
We shall take 


E<U<2lmée, (3) 


where E is the kinetic energy of the electron 
inside the box. The second condition eliminates 
the possibility of the electron’s jumping to a 
state of negative kinetic energy outside the box. 
For r<R the electron density must be finite and 
for r>R the density must decrease exponentially. 
It is then easy to show that the only possible 
solutions of the Eqs. (2) can be written 


r<R| 
¥2=iA (ar)! f;4;(ar) 

j=-il 


where @ are the values of respec- 
tively, when =0, and a@ has been taken to be 
—i(aya2)!. A and A’ are arbitrary constants. 
Outside the box we must use Bessel functions 
of the third kind which vanish exponentially as 
r increases, i.e., Hankel functions. We find 


r>R | 
(6) 

¥2= B(Ber) j=1 | 
and ¥i=B' (Bir) r>R | 
(7) 


where B= +(6:82)'. By the condition (3) 8 is 


real. B and B’ are arbitrary constants. 

Taking first the case 7=1, we must equate the 
values of (4) and (6) at the boundary, and 
eliminating the arbitrary constants A and B we 
obtain the eigenvalue equation 

(a1)*fja(ar) 
=0 (8) 


determining the energy spectrum. 

A rigorous solution of (8) will not be given, 
we shall only consider the case when BR>1. We 
assume that //’—U is not nearly equal to me’, 
in which case the electrons would leak out of the 
box rapidly. We then have 8B~mc/h, and the 
condition 8R>1 becomes, A. denoting the 
Compton wave-length, Rd. which is always 


true in any physical application. Now for 8R>1 
and Eq. (8) becomes 
f;_)(aR) =0, j 21. (9) 
When j= —1 a similar argument leads to the 
eigenvalue equation 
‘ae)  f_j.s(aR) + (8) Be)! f_; =0, 
jai. 


The eigenvalues of (9) and (10) are different, 
and correspond to opposite directions of the 
electron spin. 

The number of zeros of (9) is given by the 
number of intersection points between the real 
functions 


and (8; 82)'f;.;(aR). 


Let denote the positive zeros of 


fj-\(aR) arranged in ascending order of magni- 


tude. Then it has been shown that 


i.e., the positive zeros of f;\(a@R) are interlaced 
with those of f;,;(@R). Further the functions 
i(a,/ a2)! and (8; 82)! are real, finite and mono- 
tonic. Hence it follows that we can take the 
number of eigenvalues of (9) in the range RAa 
to be equal to the number of zeros of, say 


fisi(aR), in RAa. We shall use Debye's* asymp- 


totic formula for Bessel functions, so that its 
zeros are given by 


cos [j(tan y—y)—7, 4]=0, (12) 


where cos y=j/Ra, the formula being valid 
provided Ra>1 and 0<j/Ra<1. The first 
condition means that we have to exclude those 
electrons whose de Broglie wave-length is com- 
parable with R, corresponding to a few eigen- 
values of very low energy. The second is satisfied 
since a/27 represents the linear momentum, and 
hence Ra/27 the maximum angular momentum, 
regarding the electron classically. 

Let a change from a to a+Aa. The increase 


§ Jahnke-Emde, Tables of Functions, (1933), p. 204. 
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of the argument in (12) becomes 


da) (tan y —y)Aa= RAa sin 


and the number of zeros of (9) in the range Aa 
becomes (R/7)Aa@ sin y. A similar argument ap- 
plied to Eq. (10) yields the same number of 
eigenvalues in the range RAa. Now every energy 
level belongs to several independent states, 
corresponding to the various possible eigenvalues 
of a Cartesian component of the total angular 
momentum. The number of these states, for any 
value of j is 2|j|.4 We thus get for the total 
number AN of independent eigenfunctions in the 
range Ja, combining the contributions for posi- 
tive and negative values of /, 


AN j504j(R/m)Aa sin 
or AN = sin y COs y. 


Since cos y changes very little as 7 increases to 


4 Dirac, Proc. Roy. Soc. A118, 355 (1928). 
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j+1, we can convert the sum into an integral 
by taking 


dj d 
—=—(aR cos y)= — aR sin jy, 
dy dy 
4R%a®Aa fo 
ie., AN=—— f sin? y cos ydy, 


since y=7,/2 and y=0 for j7=1 and j—1=Ra, 
respectively. This gives 


AN =4R*a*Aa/3r. 


Inserting a= and =p 
the volume of the potential hole, we get 


This is the same formula! which is obtained by 
imposing periodic boundary conditions on the 
wave functions representing progressive waves. 

I wish to thank Dr. S. Chandrasekhar for his 
interest in the work. 
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Erratum: The Precise Measurement of Three Radium B Beta-Particle Energies 


F. T. RoGers, Jr. 


The Rice Institute, 


Houston, Texas 


(Phys. Rev. 50, 515, 1936) 


HE value of e/m, quoted on the last page of this paper is incorrect; it 


should be 


e/mo= 1.75762 X10’ e.m.u./g, 


with €,,=0.00026 X10". Thus €:/(e/m»o) =1.48/10,000 and the value of e is 
reduced, though by very little. The values of 1) computed from this correct 


(gauss cm) 


1406.1 


(electron volts) 


1.513 X 10° 
2.045 10° 
2.612 X 10° 


e/mp are given in the accompanying Table X, which should replace the Table 
X of the original paper; they cannot be considered to be more certainly accurate 
than 1 part in 1,000. 

In Eqs. (6) and (7), v should be replaced by V. 


! 
| 
I 
TABLE X. Energies of three intense Ra B 8-rays. 
1671.1 
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LETTERS, TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 


addressing them to this department. Closing dates 


for this department are, for the first issue of the 


month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Concerning the Probable Magnitude of the Next Sunspot 
Maximum 


Although the sharp increase in sunspot numbers during 
the past four months suggests the approach of a high peak 
in the 11-year cycle of activity, there is some evidence 
which seems to indicate that the average degree of activity 
at the coming maximum may not greatly exceed the 
moderate values of the previous five cycles. This evidence, 
described below, is based on a correlation appearing in the 
sunspot data for the past two centuries. These data cover a 
relatively short range of experience, and the evidence, 
therefore, is far from conclusive. It appears to hold some 
promise as an indicator, however, and experience during 
the coming peak should provide an interesting test of its 
significance. 

In Fig. 1 each sunspot maximum since 1761 has been 
plotted against the preceding minimum, using three-year 
averages of the Wolf-Wolfer sunspot numbers to deter- 
mine the maxima and minima. The results show some evi- 
dence of a correlation between the activity at a minimum 
of the 11-year cycle and the magnitude of the following 
peak. The average sunspot number for the 1932-1934 
minimum was 8.5; hence, it appears from Fig. 1 that the 
three-vear average for the approaching maximum may 
possibly fall somewhere between 60 and 90. 

Three-year averages of the sunspot numbers were used 
in Fig. 1 rather than single annual values because they give 
a better indication of the average degree of activity at the 
maximum and minimum epochs. The computed correlation 
coefficient for the data of Fig. 1 is +0.84. This is reduced 
to +0.65 when average sunspot numbers for single years 
of highest and lowest activity are compared. 

In Fig. 2 each sunspot maximum, determined by a three- 
year average, is plotted against the minimum immediately 
following it. A correlation much poorer than that of Fig. 1 


+ 128 2 
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Fic. 1. Fic. 2. 
_ Fig. 1. Relation between sunspot maxima and preceding minima. 
(S-vear averages.) Fig. 2. Relation between sunspot maxima and 
following minima. (3-year averages.) 


is clearly indicated. For this case the computed correlation 
coefficient is +0.37. While the magnitudes of these coetti- 
cients cannot be given much weight in view of the small 
number of data represented, the wide difference in the 
relationships of Figs. 1 and 2 does appear to be significant. 
The evidence in these two figures suggests that the magni- 
tude of each sunspot peak is controlled to a considerable 
extent by activity originating at the preceding minimum, 
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Fig. 3. Relation between sun- 3 
spot maxima and corresponding 3 . 
terrestrial magnetic activity § 
maxima. (3-year averages.) 
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and that this new activity is not closely related to the 
activity of the previous cycle. The abrupt change of the 
sunspots from low solar latitudes to high latitudes at each 
minimum of the cycle is a further indication that a new 
phase of solar activity may originate at each minimum. 

It has been observed that years of highest sunspot num- 
bers are not always years of greatest terrestrial magnetic 
activity, although both phenomena follow an approximate 
11-year cycle of variation. Magnetic acitivity is of particular 
interest in connection with short wave radio communica- 
tion because the more intense magnetic disturbances are 
almost invariably accompanied by disturbances to long- 
distance radio transmission. Data upon which to base fore- 
casts regarding the magnetic activity maxima are even more 
limited than for the sunspot maxima. A series of dependable 
data extending back to 1835, based on day-to-day varia- 
tions in the mean horizontal intensity of the earth’s mag- 
netic field, is given in a paper by J. Bartels. Three-year 
averages taken at the peaks in this series were compared 
with the corresponding (although not always coincident) 
sunspot peaks, with the results shown in Fig. 3. On the 
strength of this rather meager evidence it appears that if a 
moderately high three-year peak of sunspot activity does 
occur in the present cycle, it should be accompanied by an 
equally moderate peak of magnetic activity. 

A. L. DuRKEE 

Bell Telephone Laboratories, Inc., 


New York. N. Y. 
March 8, 1937. 


1 J. Bartels, Terr. Mag. 37, 1-52 (1932). 
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Correction Due to Motion of Center of Gravity in the 
Hartree Approximation in Nuclei 


The use of one-particle wave functions with a fixed origin 
in the quantum mechanics of nuclei gives a state in which 
there is a motion of the center of gravity (G, say), and thus 
introduces an error into the calculation of binding energies. 
Fliigge' has allowed for this by multiplying the kinetic 
energy hy a factor (1—1/N), where N is the number of 
particles, and the same (apparently empirical) procedure 
has been followed by Feenberg and Wigner? in their calcu- 
lation of binding energies of light nuclei. Bethe and Rose’ 
have recently considered the question, and have shown that 
with individual wave functions of the oscillator type, the 
whole wave function can, under certain restrictions, be 
expressed as the product of one-particle functions of co- 
ordinates relative to G, and a function of the coordinates 
of G itself; and that, under the same assumptions, the 
kinetic energy of G is just equal to the kinetic energy of a 
single s particle. Perhaps the following very simple con- 
siderations relating to the problem, which are applicable 
whatever type of one-particle wave functions is used, are 
not entirely superfluous. 

The conditions that a wave function y of a system of V 
particles shall describe a state in which G is at rest, are 

(pit---+py)y 

where p:, etc., are the momenta of the individual particles. 
These require that ¥ be a function of the relative coordinates 
only. If an approximate y be used which does not have this 
property, however, we can only determine the probability 
that G has various kinetic energies of motion. In this case, 
we may reasonably subtract from the calculated total 
energy the mean value of the kinetic energy of G, in order 
to give the binding energy. This mean kinetic energy is 
calculated in the usual way from the operator 


2 


N N 

Te=(1/2NM)( =7/N+1/NM pi-p; 
i=l 


where M is the mass of each particle, and T is the ¢otal 
kinetic energy. The second part of Tg is of the same form 
as the perturbation of atomic energy levels which arises 
from the finite mass of the nucleus,‘ and represents the cor- 
rection to Fliigge’s procedure. This correction will not in 
general vanish, but will do so if all particles are in the same 
shell or if exchange is neglected. Since Tg is precisely the 
same operator as that used by Bethe and Rose, expressed 
in a different type of coordinates, it must vield the same 
result in the special case considered by them. Indeed, the 
only advantage in separating the coordinates of G as they 
do, appears to be in facilitating the calculation of 7@ in 
this special case. 

One further point may be mentioned: it is evidently in- 
correct to thus reduce the kinetic energy before minimizing 
the energy, as was done by Feenberg and Wigner.’ In fact, 
if the usual variational principle be applied using this re- 
duced kinetic energy, with y unrestricted except for the 
usual conditions, the correct wave equation does not result. 
The correction to the kinetic energy must therefore be 
made after minimizing. The difference will, however, be 
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small except for the very light nuclei. 

For He‘, for example, if we use the model and force 
constants used by Feenberg and Wigner,? and if we assume 
with them the one-particle functions 

la? 
where a is the depth of the potential hole and o a variable 
parameter, we find the following results for the binding 
energy: 
uncorrected for motion of G: 11.8 Mev (¢=0.846) 
corrected, as here: 22.0 Mev (o =same) 
corrected as by Feenberg and Wigner: 25.1 Mev (¢=1.30), 
A. F. STEPHENSON 
Department of Applied Mathematics, 
University of Toronto, 


Toronto, Canada, 
February 23, 1937. 


1S. Fliigge, Zeits. f. Physik 96, 459 (1935). 

2? E. Feenberg and E. Wigner, Phys. Rev. 51, 95 (1937) 

3H. A. Bethe and M. E. Rose, Phys. Rev. 51, 283 (1937). 

*Condon and Shortley, Theory of Atomic Spectra (Cambridge, 
1935), p. 418. 


Inelastic Scattering of Fast Neutrons 


It has been found that fast neutrons excite soft gamma- 
rays (~0.8 Mev) in many elements.’ ? We have undertaken 
to determine whether the fast neutrons are absorbed during 
the interaction, as in the case of slow neutrons, or whether 
the gamma-rays are produced as the result of inelastic 
collisions of the fast neutrons with atomic nuclei. Previous 
measurements of the absorption of fast neutrons have been 
rendered ambiguous by scattering,’ or by the possibility 
that the energy spectrum of the neutrons may have been 
altered during their passage through the absorber.* We 
have measured the absorption of fast neutrons using 
cylindrical absorbers to minimize the loss of neutrons by 
scattering and using a method of detection which made any 
change in the energy spectrum of the neutrons relatively 
unimportant. 

A source of neutrons (Be +200 mg Ra) was mounted in 
the center of a cavity A (Fig. 1) within a large lead block. 
Measurements were made with and without the cvlindrical 
absorbers around the source at A. The wall thickness of the 
iron and lead cylinders was 4.2 cm and of the copper 
cylinder was 3.5 cm. The fast neutrons were detected by 
slowing them in the optimum thickness of paraffin (5.7 cm 
and counting by means of a Geiger-Miiller counter the 
gamma-rays produced by the absorption of these slow 
neutrons in cadmium. A layer of BgC (0.6 g¢/cm*) used asa 
shutter served to define the effect of these slow neutrons in 
the usual manner. 


| Absorber AWN AWN 
e | | YAN MIO \ 
| | Pb Block {| Par | 
ZAKS Muller 
Counter 
Pyrex Sheets, 4 BL 
30 x30cm 


Fic. 1. Arrangement forjmeasuring absorption‘of fast neutrons. 
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Using this method we find a transmission through Fe, Cu 
and Pb which seems slightly to exceed 100 percent. Hence 
we conclude that within the limits of the experimental 
method no absorption occurs. 

On the other hand, when the fast neutrons are detected 
by the soft gamma-rays produced by their interaction with 
lead,? a very different and striking result is obtained. In 
these experiments the paraffin and cadmium (Fig. 1) were 
removed. When the cylindrical absorbers were introduced 
at A this time a marked decrease in the count was observed. 
The transmissions of fast neutrons measured in this way 
are as follows: Cu, 76+2 percent; Fe, 78+2 percent; Pb, 
85242 percent. 

These results show that although fast neutrons are very 
largely transmitted through the absorbers, they lose much 
of their ability to excite gamma-rays in lead. Since the 
neutrons are not absorbed, the energy of the gamma-rays 
must come from the kinetic energy of the neutrons. In 
producing these gamma-rays, then, the neutron must lose 
energy and be slowed down. The results show, however, 
that below a certain energy, perhaps 0.5 Mev, the neutrons 
are no longer capable of exciting gamma-rays and hence no 
further slowing is to be expected. The fact that the trans- 
mission as measured by the first method actually exceeded 
100 percent (Fe, 102+1 percent; Pb, 104+2 percent; Cu, 
107+2 percent) can now be attributed to this marked 
slowing in the absorbers and to this extent our method of 
detection was evidently affected by changes in the energy 
of the fast neutrons. 

We conclude, therefore, that many of the collisions of 
fast neutrons with nuclei are inelastic, leaving the nucleus 
in an excited state. The energies of the gamma-rays 
subsequently emitted correspond to the energies of these 
low lying nuclear states. From the work of Kikuchi and 
associates on the energies of the gamma-rays from Cu, Fe, 
Cd, Zn, Sn and Ni,! it appears that these energy states are 
roughly the same for each of the elements investigated. 
D. C. GRAHAME 
G. T. SEABORG 
G. E. Gipson 


Department of Chemistry, 
University of California, 
Berkeley, California, 
March 15, 1937. 


1 Kikuchi, Aoki, Husimi, Proc. Phys.-Math. Soc. Japan 18, 115 
(1936). 

? Gibson, Seaborg, Grahame, Phys. Rev. 51, 370 (1937). 

3 Fleischmann, Zeits. f. Physik 97, 265 (1935). 

‘ Dunning, Phys. Rev. 48, 265 (1935). 


Potential and Luminescence of Insulated Willemite 
Cathode-Ray Screens | 


The results of experiments on two phosphors bring out 
the following results: 

(1) For electron energies less than about 175 volts, the 
phosphor charges up to cathode potential. The reason is of 
course that the secondary emission yield or ratio is less than 
one for electrons of this energy and therefore no stable 
potential can be maintained in the absence of gas ions or 
some other discharging agent. 


(2) For electron energies higher than Vynin 175 volts 
the phosphor takes on a potential equal to or slightly 
positive with respect to the anode if the phosphor coating is 
very thin, while if it is thick and rough the surface potential 
may be as much as 10 volts negative with respect to the 
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anode and yet be stable. For the thin coating the surface 
potential rises to as much as 40 volts positive with respect 
to the anode at about 1000 volts and then falls again with 
increasing electron energy. 

(3) For electron energies of about 5000 volts and greater 
the secondary emission yield again approaches a ratio of 
unity and the phosphor potential tends to remain at a 
fixed value relative to the cathode and therefore the anode to 
phosphor potential difference increases linearly with the 
anode potential, the latter being taken relative to the 
cathode. For the phosphors studied it is therefore impos- 
sible to bombard the screen with electrons of an energy in 
excess of about 6000 volts even though the anode is 
maintained at 10,000 volts. Temporary increases in 
secondary emission as indicated by screen potentials, 
relative to cathode, greater than 6000 volts were produced 
by evaporating a trace of thorium from the cathode. This 
brought the screen potential to about 9000 for an anode 
potential of 10,000 volts. 

(4) Luminosity measurements taken with a photoelectric 
cell showed photocurrents more or less proportional to the 
anode voltage but when the screen potential is used 
instead of the anode potential the following law accurately 
represents the observed results for the thick screen when 
values of 1( V,— Vo)? are less than 5.0. 

L=Qi(V,— Vo)". 
This result depends on the assumption that the spectral 
distribution of the light is independent of current and 
voltage. 
L=photoelectric cell current, Q=constant, 

i=primary current density in amp. per sq. cm to the 

screen, 
V,=potential difference screen to cathode in volts, and 
Vo=‘‘dead” voltage of about 350 volts. 
This dead voltage is not related to the voltage V,,jn. The 
latter is determined by secondary emission properties, while 


" 


the former is no doubt an average property of the surface 
layer on each individual crystal. For screen potentials 
greater than 1000 volts and less than 4000 volts the 
luminosity is accurately proportional to the current for 
smaller current densities. As the current density increases 
the light falls below that expected assuming true linearity 
an amount which is a function of the light intensity only. 
The current at which a given deficiency in light is observed 
decreases as the screen voltage is increased. For screen 
potentials greater than 4000 volts the failure of the linear 
relationship to current is less marked. 

(5) It has been found experimentally that the evapora- 
tion of thorium atoms onto the screen results in enhanced 
secondary emission at high voltages. Although the thorium 
lost its effectiveness in a few hours’ time, it seems that the 
best solution to the problem of high screen luminosity must 
come throdgh increasing the surface Secondary emission 
efficiency so as to permit the effective use of higher voltages 
since for a given power input to the screen the luminosity 
should increase in proportion to the voltage according to 
these observations. 

W. B. NorrincHam 

Massachusetts Institute of Technology, 


Cambridge, Massachusetts, 
February 6, 1937. 
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An Application of the Theory of Neutron Absorption in 
Aqueous Solutions 


Recently! the writer suggested that the constants describ- 
ing an element’s resonance absorption of slow neutrons— 
energy Eo and half-width I of the resonance, and maximum 
cross section ores—might be determined by measurements 
of the activation of an aqueous solution. In this way one 
could avoid geometrical complications by using a volume 
of homogeneous solution large enough for practically all the 
neutrons to be absorbed in it. The basis of the method is 
the fact that the number of atoms per second activated by 
resonance neutrons, 


Ares = Of = 208 = Eo) {1+ Nadcres} (1) 


(Q=number of neutrons/sec. leaving source; 
\=mean free path for scattering by protons), 
‘ 


is a nonlinear function of the detector concentration Ny 
(atoms/cm). This makes it possible to calculate T'/Eo and 
Nores from the values of ¢ for two Ny's. The validity of (1) 
depends on the conditions (T/Eo)<1, 

The method as originally described had the disadvantage 
of requiring the use of very large concentrations of cad- 
mium and boron in solution along with the detector. In 
order to insure against serious errors it would be necessary 
to take readings at various concentrations of cadmium 
and establish an extrapolation to infinite concentration. 
The boron concentration also required would be so high 
that it would present a considerable technical difficulty. 
In this note a method is described which avoids these 
difficulties and has besides the interesting feature that 
it does not depend on the assumption of the v~ law for boron, 
which is involved in methods hitherto used. This assump- 
tion is here replaced by the assumption that ocap for the 
detector accurately obeys the Breit-Wigner formula. 

In the presence of a concentration Ng (atoms, cm) of 
boron the total activation is 


A =ApestA therm = Nar/(Ne+Nart+s) (2) 
where 
r= Av/(oBv) Av 
7=mean lifetime for capture of thermal neutrons by 


protons. (The subscript Av indicates an average over the 
Maxwell velocity distribution.) From (2) we have 


-AG+H=A Np (3) 

with 
s+Ngr=G (4) 
NaQr+ (Qs) -s=H. (5) 


If A is measured for Ng=0 and for two other values of Nz, 
(3) yields a set of equations which can be solved for the 
three unknowns (Q¢), G, and H. (The A’s will of course be 
in arbitrary units, and (Q¢) and H/Ng will then be in the 
same units.) If this procedure is carried out for two different 
values of Ng, (4) then enables us to solve for s and r, after 
which Q can be determined from (5). Dividing the two 
values of (Q¢) by Q, we get the values of ¢ for two Na’s, 
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from which Aores and T, Eo can be determined as remarked 
after (1). Also since by the Breit-Wigner formula 


Av >= res(T Eo)? (6 
we have 
vo=(T'/ Eo) ay / Ores. 7) 


Thus by six relative activation measurements we can get 
Eo, T, and opes in terms of the two quantities \ and (og-v)4,. 

For this procedure to be practicable, and lead to results 
whose percent errors are not much greater than those of 
the activation measurements, the concentrations used 
must satisfy the following requirements: 


Np’: (v0/d) (8) 
Na'r2s (9) 
(10 


Here a prime indicates the larger of the two quantities. 
The smaller Ng should be about half of Ng’; the smaller 
Na should be half of N4’, or less, and should not satisfy 
(10). The boron concentration required is of the order 
(¢’/x) compared to that required by the method previously 
described. 

In water, the mean number of elastic collisions made by 
a thermal neutron before it is captured is quite large. It 
can still be regarded as fairly large in a solution satisfying 
(8), (9), (10), provided 


(v% ‘Vtherm) «<1. (1 1) 


If this is so, the Maxwell distribution law is obeyed rather 
accurately, and (2) correctly represents the competition 
between the different absorbers of slow neutrons, even if 
there is a marked departure from the v™ law. Thus cad- 
mium could be used instead of boron, if preferable on tech- 
nical grounds. 

The method can be applied to cases in which the detector 
substance has another low energy resonance besides the one 
being studied, provided that: (a) The resonance widths are 
small compared to the resonance energies and _ their 
difference; (b) the two absorptions give rise to different 
activities, so that that due to the one studied can be meas- 
ured separately, One has to apply to the results of the calcu- 
lations described above the correction factors f/f, 
(1+r’/r)>, and (1+r’/r), to obtain the actual values of 
Ores, (I'/ Eo), and vo, respectively. Here f is the fractional 
concentration of the isotope to which the activity is as- 
cribed, and r’ is the constant characterizing the thermal 
absorption connected with the other resonance (cf. defini- 
tion of r). The ratio r’/r cannot be obtained with any accu- 
racy from such measurements; it should, however, be easy 
to determine from measurements made with solid absorbers. 
A similar extension to cases with more than two low energy 
resonances is possible provided conditions corresponding to 
(a) and (b) above are satisfied. 

W. H. Furry 

Department of Physics, 

Harvard University, 


Cambridge, Massachusetts, 
March 3, 1937. 


1W. H. Furry, Phys. Rev. 50, 381 (1936). 
? E. Amaldi and E. Fermi, Phys. Rev. 50, 910 (1936). 
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LETTERS TO 


A Theory of Elasticity 


I write to call to the attention of any of your readers who 
may be interested, the fact that there is given in the 
April (1937) issue of the American Journal of Mathematics 
(Johns Hopkins Press, Baltimore, Md.), a theory of elas- 
ticity which does not make any assumptions as to the 
infinitesimal nature of the strain. The following summary 
of the paper will give an idea of the results to be found in it: 

Formulae are derived which enable one to calculate the 
stress in an elastic medium when the strain and the elastic 
energy density are known, no simplifying assumptions, 
such as smallness of strain, being necessary. For an iso- 
tropic elastic solid under hydrostatic pressure the following 
one constant formula gives good agreement with experi- 
mental observation (only two elastic constants A, wu being 
used in the expression for the elastic energy density) 
p=a(f+5f*); f=43{(Vo/V)I-1}; a=3\+2y, 
In the Young's modulus experiment the formula for the 
extensional stress (again using only the two constants A, 


is 
T=Ees 1- E (Young’s modulus) = ————— 
J 


where «=(2+e)e (1+e)*, e being the relative extension. 
Hence 7 has a maximum value (A+y)E£/4(2A+3,u) occur- 
ring when e=(A+yu) 2(2A+3u). For a true second-order 
approximation (the infinitesimal theory being regarded as 
a first-order approximation) five elastic constants occur 
and the corresponding formulae are either given or their 
derivation is immediate. 
F. D. MuRNAGHAN 
School of Mathematics, 
Institute for Advanced Study, 


Princeton, New Jersey, 
March 9, 1937. 


The Disintegration of Li‘ 


It was observed by Lewis, Burcham and Chang! that the 
B-ray disintegration of Li® is accompanied by a further 
disintegration of the residual Be* nucleus into two a- 
particles. In our letter to the Physical Review? we considered 
this reaction from a theoretical point of view. We tried to 
account for the discrepancy between the limit of the 8-ray 
spectrum? of Li® and the mass difference between Li*® and 
Be’, as derived by Rumbaugh and Hafstad‘ from the 
nuclear reaction. 


Li7+H?=Li8+H'. 


Selection rules for 8 disintegration indicate that the Be‘ 
nucleus is left after the emission of the electron in an 
excited state. We stated that this excited state should emit 
y-tays. However, estimates making use of penetration 
through the potential barrier for d waves show that the 
disintegration into a-particles should be more probable. 
Indeed, this probability is so large that it may be expected 
to broaden the excited Be® level considerably. Such a 
broadening has in fact been indicated in observations of 
Cockcroft and Lewis* in the reactions 
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= (Be*)*+ Het =3Het 
and especially in the work of Dee and Gilbert® on 
BU+H!' = 


The experiments seem to indicate a half-value breadth of 
roughly 1 Mev for (Be*)* which is most probably a 'D level. 
These experimentally observed widths are in fair 
agreement with theoretical estimates on the basis of 
Gamow'’s formula for a d wave so this interpretation of it 
appears to be correct. The width due to y-radiation from 
the Be®'D level should be smaller than 1 volt since it 
corresponds to 4?=2. Thus the y-radiation should be 
negligible in comparison with the alpha-particle emission. 
The long range alpha-particles (5 cm) are presumably due 
to transitions in a higher level. 
G. BReEIT 
WIGNER 
University of Wisconsin, 


Madison, Wisconsin, 
February 22, 1937. 


'W. B. Lewis, W. E. Burcham and W. Y. Chang, Nature 139, 24 
(1937). 

?G. Breit and E. Wigner, Phys. Rev. 50, 1191 (1936). 

*H. R. Crane, L. A. Delsasso, W. A. Fowler, C. C. Lauritsen, Phys. 
Rev. 47, 971 (1935). 

‘L. H. Rumbaugh and L. R. Hafstad, Phys. Rev. 50, 681 (1936) 

J. D. Cockcroft and W. B. Lewis, Proc. Roy. Soc. A154, 246 (1936). 
P. I. Dee and C. W. Gilbert, Proc. Roy. Soc. AlS4, 279 (1936), 
cf. especially §6. 


The Common Ion Effect in Some Aqueous Solutions As 
Shown by Means of the Raman Effect 


The suppression of ionization of zinc chloride in aqueous 
solution by adding the common chloride ion is demon- 
strable by means of the Raman effect. The binding between 
the zinc and chlorine atoms is homopolar in concentrated 
solutions and heteropolar in dilute solutions. Consequently 
any suppression of ionization will increase the intensity of 
the Raman line at A? 280 corresponding to the symmetrical 
vibration. 

In Fig. 1 are shown microphotometer tracings of solutions 


Fig. 1. The common ion effect in solution as demonstrated by means 
of the Raman effect. Curve 1, 1 molal ZnCle; curve 2, mola! ZnCls; 
curve 3, } molal ZnCl», 1 molal NaCl; curve 4, } molal ZnCle, 2 molal 
NaCl; curve 5, § molal ZnCle, 2 molal NaC 
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of zine chloride at different concentrations and of solutions 
consisting of mixtures of zinc chloride and sodium chloride. 
The Raman shift at A? 280 in 1 molal and } molal solutions 
decreases in intensity more rapidly than can be accounted 
for by the dilution effect alone. It is evident that there is an 
increased ionization accompanying dilution. However, in a 
; molal solution in respect to zinc chloride and 2 molal in 
respect to sodium chloride the zinc chloride line appears 
with considerable intensity. A } molal zinc chloride solution 
shows no Raman shift at all. The solution of } molal zinc 
chloride and 2 molal sodium chloride gives rise to a shift of 
greater intensity than observed from 1 molal zinc chloride 
solution. This would seem to indicate that the 1 molal 
solution was at least fifty percent ionized. 
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In addition to the above effects there is a notable altera- 
tion in intensity of the water band occurring at approxi- 
mately Av 150. This is attributed to the hindered transla- 
tional motion of the water molecules. This shift is apparently 
intensified by sodium chloride and diminished by zinc 
chloride although the maximum intensity occurs at concen- 


trations of } molal zinc chloride and 1 molal sodium 


chloride. 

No other line than the one mentioned for zinc chloride 
was observed with excitation by means of the 25364 
mercury line. 

James H. Hissen 


Geophysical Laboratory, 
Carnegie Institution of Washington, 
Washington, D. C., 
February 10, 1937. 


Proceedings of the New England Section of the American Physical Society 


MINUTES OF THE CAMBRIDGE MEETING, FEBRUARY 6, 1937 


HE tenth regular meeting of the New 

England Section was held at the Massa- 
chusetts Institute of Technology, Cambridge, 
Massachusetts, on Saturday, February sixth. 
The presiding officer was Professor F. A. 
Saunders, chairman of the section. 

The afternoon session was devoted to two 
invited papers, one on ‘Project Methods for 
College Students” by Professor S. R. Williams 
of Amherst, and the other on ‘Some Recent 
Work on Isotopes” by Professor K. T. Bain- 
bridge of Harvard. 

The meeting was attended by about 50 


members. 
The program of the section consisted of 12 
papers, the abstracts of which are given below. 
On the previous evening, a small group from 
the section were guests of the Harvard College 
Observatory at an open night. Cloudy weather 
prevented the observation of celestial objects, 
but the buildings were open for inspection, and 
Professor B. J. Bok gave an outline of the 
milky way research problems being studied at 
the Harvard Observatory. 
JosEru C. Boyce 
Secretary-Treasurer 


ABSTRACTS 


1. The Application of Paschen’s Law to the Reignition 
of an Arc. R. B. Power anv J. D. Cosine, Harvard 
Graduate School of Engineering.—An a.c. arc, in order to 
restrike after passing through current zero, requires a 
potential considerably higher than the normal burning 
potential. This reignition potential is investigated for short 
gaps in nitrogen, using pure graphite electrodes with 
spacings up to two mm and pressures up to 500 cm Hg. 
The reignition potential is found to have two character- 
istics. One of these characteristics is followed for the ares in 
which the cathode spot is maintained by field emission, and 
the other is followed for the ‘‘thermionic” arc. Relations are 
obtained between this reignition potential and both the gas 
pressure and the gap separation, with the are current as 
parameter. These’ relations may then be combined into a 
single relation between the reignition potential and the 
product of the pressure and spacing. This gives a function 


that is of the same form as the Paschen law for the initial 
sparking potential of a gap. The constants of this function 
are dependent on the arc current, the gas in the space, and 
the constants of the circuit. That the constants of the 
reignition function are not the same as those for the 
sparking potential law is to be expected from the very 
dissimilar conditions applying in the space. Extreme purity 
of electrode materials and of gas, as well as extreme care in 
obtaining the experimental results, are essential. 


2. Preliminary Observations of Absorption Spectra in 
the Schumann Region. Frep W. Pau, Massachussetts 
Institute of Technology.—An electric furnace with a tung- 
sten tube has been used in conjunction with a 2-meter 
grazing incidence vacuum spectrograph to observe the 
absorption spectra in the Schumann region. A charge of 
silver placed in the tungsten tube was heated to the 
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desired temperature by the passage of a large variable 
current through the tungsten tube. The source of con- 
tinuous light employed was a Lyman capillary discharge of 
the conventional type. Observations have been made 
between 41300 and A1930A. Spectra obtained at several 
temperatures show an increase of intensity of certain lines 
with increasing vapor pressure. It is believed that some of 
the lines of wave-length greater than \1636 belong to the 
principal series of silver, and that the remainder belong to 
series of the type described by Beutler* for other elements. 


* H. Beutler, Zeits. f. Physik 87, 176 (1933-34). 


3. Temperature Shift of the Potassium Resonance Lines. 
Gorpon F. Hutt, Jr. Yale University.—Effect of tempera- 
ture on the shift and half-width of the potassium resonance 
lines in absorption perturbed by nitrogen pressures from 1 
to 15 atmospheres has been investigated. Two sets of data, 
one at 270°C and the other at 90°C, show that both shift 
and half-width vary with temperature in the same manner. 
As the “relative density”’ of nitrogen is increased, the effect 
of temperature decreases, and becomes unobservable above 
relative density 7. The results are in agreement with 
Margenau’s theory* which predicts the presence of both 
the Lorentz type of broadening dependent upon tempera- 
ture, and ‘‘statistical” broadening independent of tempera- 
ture for low relative densities. The temperature effect 
decreases with increasing relative density, because of 
statistical broadening. 


*H. Margenau, Phys. Rev. 48, 755 (1935). 


4. Pressure Effect in Bands of Several Dipole Mole- 
cules. S. D. Cornett, Yale University—The lines of 
several photographic infrared bands of three dipole mole- 
cules, HO, HCN and NH;, have been investigated for 
broadening due to increased pressure of the vapors them- 
selves. The spectrograms, taken at high dispersion, were 
microphotometered at a magnification of 40 : 1. Below one 
atmosphere pressure, the broadening per atmosphere, in 
excess of that attributable to dispersion forces, differs 
considerably for H»O and HCN. The half-widths found for 
one atmosphere are : 0.65 cm for H,O, 1.56 cm for HCN, 
1.45 cm™ and 0.81 cm for two NH; bands. The ratio of 
excess broadening for H,O and HCN is of the order of the 
ratio of wt for the two, in agreement with theoretical 
predictions. The NH; breadths are larger than indicated by 
the » of NH, perhaps because of the greater possibility of 
resonance between the multiple levels of the symmetrical 
top. For HCN the lines at and near the maximum of the 
rotational velocity distribution are broader than lines of 
lower and higher J values, in accord with theoretical 
predictions of Watson and Margenau (Phys. Rev. 51, 48 
(1937)). The slit-width, for slits of the widths generally 
used, is found to have a marked effect on the observed rate 
of broadening with pressure. 


5. Mass-Spectrographic Measurement of the C+ Band 
Resulting from the Dissociation of CO*. K. T. BAINBRIDGE 
AND E. B. Jorpan, Harvard University—A band may 
appear on the recording plate when CO+* is dissociated into 


C* and O in the region between the electrostatic and 
magnetic deflection chambers of a mass spectrograph. The 
calculated apparent mass of the carbon fragment is 
C® x C®/C®+-0!6 = 5.145; mass units on the physical mass 
scale. Energy lost in dissociation and excitation might act 
to reduce this figure. Aston* has reported a band produced 
as above with the anomalous mass value 5.18, approxi- 
mately 0.7 percent greater than the theoretical value. An 
independent examination of the C* band has been made 
using a double focusing mass spectrograph in which the 
scattering gas was introduced into the tube connecting the 
electric and magnetic deflection chambers. 5.145+0.002 
mass units was obtained for the position of the band, using 
DH! and trebly ionized oxygen for comparison masses. 


* F. W. Aston, Mass-Spectra and Isotopes, pp. 60-62. 


6. Inelastic Collision of Deuteron and Deuteron. L. I. 
ScuiFF, Massachusetts Institute of Technology.—The theo- 
retical excitation function and angular distribution of the 
products of the reactions 


H?+H*—He?+n, 


are calculated using symmetrized wave functions and 
general exchange interactions, for comparison with the 
recent experiments of Kempton, Browne and Maasdorp.* 
The experimental curves for the angular distribution of the 
resultant neutrons or protons in a coordinate system in 
which the center of mass is at rest can be represented by: 
A(1+B cos? 6), where A and B may depend on incident 
deuteron energy. It is shown that the theory predicts an 
angular distribution of this form, but that the theoretical 
quantity B is too small and has the wrong dependence on 
incident deuteron energy. This discrepancy between theory 
and experiment seems to be of a fundamental nature and is 
not dependent on the details of the assumed interaction or 
of the calculation. Curves of excitation function and 
angular distribution will be shown. Theoretital and 
experimental yield curves. 


* Proc. Roy. Soc. A157, 386 (1936). 


7. Application of the Direct-Fusion Furnace to the 
Helium Method of Determining the Geologic Ages of 
Rocks. RosLey D. EvANs AND CLARK GOODMAN, Massa- 
chusetts Institute of Technology.—Preliminary tests indicate 
that the direct-fusion type of furnace* with only minor 
modifications can be used in the determination of the 
helium content of high melting point rocks. The helium 
contained in the gases released from the boiling rocks at 
2000°C is separated and quantitatively determined by the 
Paneth-Urry Method? using a sensitive McLeod gauge. A 
knowledge of the helium content, together with separate 
determinations of the uranium and the thorium content or 
of the total rate of alpha-emission, enables the geologic age 
of the rock to be ascertained. 

*R. D. Evans, R. S. I. 6, 99-112 (1935). 


i Paneth and W. D. Urry, Zeits. f. physik. Chemie 152, 110-149 
(1931). 


etls 
ng- 
‘ter 
the 
of 
the 
> 


594 AMERICAN PHYSICAL SOCIETY 


of zine chloride at different concentrations and of solutions 
consisting of mixtures of zinc chloride and sodium chloride. 
The Raman shift at Ad 280 in 1 molal and } molal solutions 
decreases in intensity more rapidly than can be accounted 
for by the dilution effect alone. It is evident that there is an 
increased ionization accompanying dilution. However, in a 
} molal solution in respect to zinc chloride and 2 molal in 
respect to sodium chloride the zinc chloride line appears 
with considerable intensity. A } molal zinc chloride solution 
shows no Raman shift at all. The solution of } molal zinc 
chloride and 2 molal sodium chloride gives rise to a shift of 
greater intensity than observed from 1 molal zinc chloride 
solution. This would seem to indicate that the 1 molal 


solution was at least fifty percent ionized. 


In addition to the above effects there is a notable altera- 
tion in intensity of the water band occurring at approxi- 
mately Av 150. This is attributed to the hindered transla- 
tional motion of the water molecules. This shift is apparently 
intensified by sodium chloride and diminished by zinc 
chloride although the maximum intensity occurs at concen- 
trations of } molal zinc chloride and 1 molal sodium 
chloride. 

No other line than the one mentioned for zinc chloride 
was observed with excitation by means of the 2536A 


mercury line. 
James H. Hippen 
Geophysical Laboratory, 
Carnegie Institution of Washington, 
Washington, D. C., 
February 10, 1937. 


Proceedings of the New England Section of the American Physical Society 


MINUTES OF THE CAMBRIDGE MEETING, FEBRUARY 6, 1937 


HE tenth regular meeting of the New 

England Section was held at the Massa- 
chusetts Institute of Technology, Cambridge, 
Massachusetts, on Saturday, February sixth. 
The presiding officer was Professor F. A. 
Saunders, chairman of the section. 

The afternoon session was devoted to two 
invited papers, one on ‘Project Methods for 
College Students” by Professor S. R. Williams 
of Amherst, and the other on “Some Recent 
Work on Isotopes” by Professor K. T. Bain- 
bridge of Harvard. 

The meeting was attended by about 50 


members. 
The program of the section consisted of 12 
papers, the abstracts of which are given below. 
On the previous evening, a small group from 
the section were guests of the Harvard College 
Observatory at an open night. Cloudy weather 
prevented the observation of celestial objects, 
but the buildings were open for inspection, and 
Professor B. J. Bok gave an outline of the 
milky way research problems being studied at 
the Harvard Observatory. 
Joseru C. Boyce 
Secretary-Treasurer 


ABSTRACTS 


1. The Application of Paschen’s Law to the Reignition 
of an Arc. R. B. Power anv J. D. Cosine, Harvard 
Graduate School of Engineering.—An a.c. arc, in order to 
restrike after passing through current zero, requires a 
potential considerably higher than the normal burning 
potential. This reignition potential is investigated for short 
gaps in nitrogen, using pure graphite electrodes with 
spacings up to two mm and pressures up to 500 cm Hg. 
The reignition potential is found to have two character- 
istics. One of these characteristics is followed for the arcs in 
which the cathode spot is maintained by field emission, and 
the other is followed for the ‘‘thermionic”’ arc. Relations are 
obtained between this reignition potential and both the gas 
pressure and the gap separation, with the arc current as 
parameter. These relations may then be combined into a 
single relation between the reignition potential and the 
product of the pressure and spacing. This gives a function 


that is of the same form as the Paschen law for the initial 
sparking potential of a gap. The constants of this function 
are dependent on the arc current, the gas in the space, and 
the constants of the circuit. That the constants of the 
reignition function are not the same as those for the 
sparking potential law is to be expected from the very 
dissimilar conditions applying in the space. Extreme purity 
of electrode materials and of gas, as well as extreme care in 
obtaining the experimental results, are essential. 


2. Preliminary Observations of Absorption Spectra in 
the Schumann Region. Frep W. Pavut, Massachussetts 
Institute of Technology.—An electric furnace with a tung- 
sten tube has been used in conjunction with a 2-meter 
grazing incidence vacuum spectrograph to observe the 
absorption spectra in the Schumann region. A charge of 
silver placed in the tungsten tube was heated to the 
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desired temperature by the passage of a large variable 
current through the tungsten tube. The source of con- 
tinuous light employed was a Lyman capillary discharge of 
the conventional type. Observations have been made 
between A1300 and \1930A. Spectra obtained at several 
temperatures show an increase of intensity of certain lines 
with increasing vapor pressure. It is believed that some of 
the lines of wave-length greater than 1636 belong to the 
principal series of silver, and that the remainder belong to 
series of the type described by Beutler* for other elements. 


* H. Beutler, Zeits. f. Physik 87, 176 (1933-34). 


3. Temperature Shift of the Potassium Resonance Lines. 
Gorvon F. Hutt, Jr. Yale University.—Effect of tempera- 
ture on the shift and half-width of the potassium resonance 
lines in absorption perturbed by nitrogen pressures from 1 
to 15 atmospheres has been investigated. Two sets of data, 
one at 270°C and the other at 90°C, show that both shift 
and half-width vary with temperature in the same manner. 
As the “‘relative density”’ of nitrogen is increased, the effect 
of temperature decreases, and becomes unobservable above 
relative density 7. The results are in agreement with 
Margenau’s theory* which predicts the presence of both 
the Lorentz type of broadening dependent upon tempera- 
ture, and “‘statistical’’ broadening independent of tempera- 
ture for low relative densities. The temperature effect 
decreases with increasing relative density, because of 
statistical broadening. 


* H. Margenau, Phys. Rev. 48, 755 (1935). 


4. Pressure Effect in Bands of Several Dipole Mole- 
cules. S. D. Cornett, Yale University—The lines of 
several photographic infrared bands of three dipole mole- 
cules, HO, HCN and NHs;, have been investigated for 
broadening due to increased pressure of the vapors them- 
selves. The spectrograms, taken at high dispersion, were 
microphotometered at a magnification of 40 : 1. Below one 
atmosphere pressure, the broadening per atmosphere, in 
excess of that attributable to dispersion forces, differs 
considerably for H2O and HCN. The half-widths found for 
one atmosphere are : 0.65 cm for H2O, 1.56 cm for HCN, 
1.45 cm™ and 0.81 cm for two NH; bands. The ratio of 
excess broadening for H2O and HCN is of the order of the 
ratio of w* for the two, in agreement with theoretical 
predictions. The NH; breadths are larger than indicated by 
the » of NHs;, perhaps because of the greater possibility of 
resonance between the multiple levels of the symmetrical 
top. For HCN the lines at and near the maximum of the 
rotational velocity distribution are broader than lines of 
lower and higher J values, in accord with theoretical 
predictions of Watson and Margenau (Phys. Rev. 51, 48 
(1937)). The slit-width, for slits of the widths generally 
used, is found to have a marked effect on the observed rate 
of broadening with pressure. 


5. Mass-Spectrographic Measurement of the C*+ Band 
Resulting from the Dissociation of CO*+. K. T. BAINBRIDGE 
AND E. B. Jorpan, Harvard University—A band may 
appear on the recording plate when CO? is dissociated into 


C* and O in the region between the electrostatic and 
magnetic deflection chambers of a mass spectrograph. The 
calculated apparent mass of the carbon fragment is 
C® x = 5.145; mass units on the physical mass 
scale. Energy lost in dissociation and excitation might act 
to reduce this figure. Aston* has reported a band produced 
as above with the anomalous mass value 5.18, approxi- 
mately 0.7 percent greater than the theoretical value. An 
independent examination of the C* band has been made 
using a double focusing mass spectrograph in which the 
scattering gas was introduced into the tube connecting the 
electric and magnetic deflection chambers. 5.145+0.002 
mass units was obtained for the position of the band, using 
D.*H! and trebly ionized oxygen for comparison masses. 


*F. W. Aston, Mass-Spectra and Isotopes, pp. 60-62. 


6. Inelastic Collision of Deuteron and Deuteron. L. I. 
ScuiFF, Massachusetts Institute of Technology.—The theo- 
retical excitation function and angular distribution of the 
products of the reactions 


are calculated using symmetrized wave functions and 
general exchange interactions, for comparison with the 
recent experiments of Kempton, Browne and Maasdorp.* 
The experimental curves for the angular distribution of the 
resultant neutrons or protons in a coordinate system in 
which the center of mass is at rest can be represented by: 
A(1+B cos? 6), where A and B may depend on incident 
deuteron energy. It is shown that the theory predicts an 
angular distribution of this form, but that the theoretical 
quantity B is too small and has the wrong dependence on 
incident deuteron energy. This discrepancy between theory 
and experiment seems to be of a fundamental nature and is 
not dependent on the details of the assumed interaction or 
of the calculation. Curves of excitation function and 
angular distribution will be shown. Theoretital and 
experimental yield curves. 


* Proc. Roy. Soc. A157, 386 (1936). 


7. Application of the Direct-Fusion Furnace to the 
Helium Method of Determining the Geologic Ages of 
Rocks. Rostey D. EvANs AND CLARK GoopMAN, Massa- 
chusetts Institute of Technology.—Preliminary tests indicate 
that the direct-fusion type of furnace* with only minor 
modifications can be used in the determination of the 
helium content of high melting point rocks. The helium 
contained in the gases released from the boiling rocks at 
2000°C is separated and quantitatively determined by the 
Paneth-Urry Method? using a sensitive McLeod gauge. A 
knowledge of the helium content, together with separate 
determinations of the uranium and the thorium content or 
of the total rate of alpha-emission, enables the geologic age 
of the rock to be ascertained. 

* R. D. Evans, R. S. I. 6, 99-112 (1935). 
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8. Longitudinal Effect in Rochelle Salt. WaLter G. - 


Capy, Wesleyan University —The fact that plates can be 
cut from Rochelle salt crystals in such orientation as to 
become polarized in a direction parallel to the direction of 
compression, and conversely to undergo extensional strain 
in a direction parallel to an applied electric field, seems to 
have been lost sight of, although the essential equations, 
applicable to all piezoelectric crystals, were given by Voigt. 
Experiments are here described with a plate cut with its 
normal making equal angles with the three crystalline 
axes, also with a bar whose length is in the same direction. 
The piezoelectric constant for compression along this 
direction has been measured, and both the plate and tod 
have been made to serve as resonators, in accordance with 


theory. 


9. The Elastic Hysteresis in Wires, Sheets and Yarns. 
IrvING J. SAXL, Providence, R. I.—The stiffness, resiliency 
and recovery in wires, sheets and yarns is a factor that is of 
importance for textile, paper, wire and other products. 
Quantitative methods for measuring elastic after-effects in 
a manner that is adaptable to routine laboratory investi- 
gations have been developed therefore. The instrument 
with which these determinations have been made will be 
shown, and the characteristics of elastic hysteresis, stiffness 
and after-effects will be discussed for a number of different 
materials. The method has found practical application for 
the evaluation and the development of new textile ma- 
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terials and other products in which exact means for 
determining elastic recovery, stiffness and similar factors 
were put for the first time on a precise scientific basis. 


10. Measurement of the Velocity of Light. Witmer C. 
ANDERSON, Harvard University.—A high frequency method 
for measuring the velocity of light with a relatively short 
base line has been developed. A Kerr cell is used for 
modulating a light beam at very high frequencies and a 
photoelectric cell is used for detection of the modulations. 
The modulated light beam is divided into two portions by 
means of a half-silvered mirror and the two parts sent over 
different optical paths. By noting the difference in path 
length required to produce a minimum of voltage in the 
tuned circuit of the photoelectric cell, and the frequency of 
modulation of the light, it is possible to measure the 
velocity of light with considerable precision. Modifications 
of this method will be described and the theory discussed 
briefly. The mean result of a series of measurements 
involving 651 observations is given as 299,764 km/sec. 
in vacuum. 


11. Project Methods for College Students. S. Rk. 
Witurams, Amherst College. 


12. Some Recent Work on Isotopes. K. T. BAINBRIDGE, 
Harvard University. 
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MEETING OF MARCH 12, 1937 


HE second meeting of the Metropolitan Section of the American Physical 
Society for the season 1936-1937 was held on Friday, March 12, 1937 in 
the Pupin Physics Laboratories of Columbia University. The following papers 


were presented. 
Afternoon Session 3:00 P.M. 


The Effect of Similarity of Particles on Gas Kinetic Quantities, Orro HALPERN 


Nuclear Moments, I. I. 


The Structure of Organic Crystals, L. H. GerMER 
Energy Levels of Fast Electrons in a Periodic Lattice, J. C. SLATER 


Evening Session 8:00 P.M. 


The Nature of Cosmic Rays, W. F. G. SwANN 


The following officers were elected for the season 1937-1938: 


Chairman 
Vice Chairman 
Secretary-Treasurer 


Elected Members of the Executive Committee 


W. BLEAKNEY 

C. J. Davisson 

W. S. Gorton 
A. MEssENGER 
Dana P. MITCHELL 


W. S. Gorton 
Secretary-Treasurer 


